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The Mathematics of the Physical Properties of Crystals 
By WALTER L. BOND 
SECTION 1 
INTRODUCTION 
HE use of crystals as oscillating elements and as light valves in electric 
circuits has given the mathematics of crystalline media an engineering 
importance. Soon after the first simple quartz oscillators were made it was 
noticed that some ways of cutting the block from the natural crystal gave 
lower temperature coefficients of frequency than did other ways. This led 
to studies of the change of elastic modulii with direction and temperature 
and finally to the discovery that there are directions in quartz for which the 
shear modulus does not change with temperature. 

Such computations are rather involved, and there is, in the English 
language, no general reference book on these new problems. The existing 
works were evidently not written with the idea in mind that anyone would 
ever actually do much numerical work with directional properties of crystals, 
since the methods used are not the best suited to this. The matrix algebra 
has the advantages of a symbolic algebra and is also, through the concept of 
matrix multiplication, a scheme for computing numerical results. 

As the problem of temperature coefficients of frequency involves the 
temperature coefficient of expansion, the temperature coefficient of density 
and the temperature coefficient of elastic modulii, these problems must be 
put into the language of matrix algebra so that they will fit into the general 
structure being built for more difficult problems. For this reason, after an 
introduction to the idea of linear vector functions, through consideration 
of the relation between the electric field and the induction in a crystal, and 
a hasty sketch of symmetry types found in crystals, we proceed to the 
consideration of stress and strain and their relations to each other. 

Following these, we take up piezo electricity and the converse piezo 
electric effects; these are important as they tell us the ways a crystal may be 
driven. We have not seen anywhere a general proof that the modulii of 
the converse effect are the same numbers as the modulii of the direct effect 
~to the first order of small quantities, though Lippman predicted the 
converse effect and demonstrated its magnitude to be about this; he ap- 
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parently didn’t consider the general case of six simultaneous stress 
ponents, six strain components, three electric field components and 
induction components. The fact that the mentioned relation is true 
to the first order of small quantities seems to have escaped the attenti 
some experimenters who have sought to show non-linearity of the p 
electric effect by demonstrating non-linearity in the converse effect. 
As a basis for light valve problems, we handle the propagation of 
through crystals, then the electro optic effect and the piezo optic et 


SECTION 2 
A LINEAR VECTOR FUNCTION 


For almost every physical constant of an isotropic medium a crystallin 
medium has several constants. For instance, a piece of glass has a 
efficient of thermal expansion but a crystal has many coefficients of thern 
expansion, the coefficient depending on direction. It might be thought 
that there were no necessary relations between the coefficients in different 
directions but there are necessary relations. 

As an example of the simplifying relations between the values of physica 
constants in different directions let us consider a crystal in an electric fiel 
Measurements show that the dielectric constant varies with direction in a 


crystal. If the field is not in the direction of greatest dielectric consta: 
the displacement current might veer over a little, much as a nail tries to 
follow the grain of the wood. We shall assume that for any electric field 
vector & there corresponds an electric induction vector D which may not 
coincide with £. Also we assume that the magnitude of D is proportiona 
to E, that is, if E results in D, then nE results in nD. Lastly, we assume 
that if F, results in D,, Fy in D....and £,, results in Dm, then FE; + E>. 
E,, results in D; + Dz + ...Dm. If these assumptions hold, then as an) 
arbitrary field E can be expressed as the vector sum of its three components 
E,, E2, E3 along three arbitrary unit vectors i, 7, k, the induction vector 
resulting from E can be computed from the induction vectors resulting from 
FE, Fo,and £3. For, let FE; result in Dyt + Daj + Dsk, Eo result in Dyot + 
Dj + Dsk and E; result in Dyi + Dj + Dyk, then Ey + EE, + EF; = fF 
results in the induction vector: 


D= (E\Diy + EsDy. + E3D\3)1 
+ (E,Da + E2D2 + E3Do3)j ..... eT 
+ (E\Ds, + E2D32 + EgDs3)k 


bo 


It is seen then, that not more than 9 constants are needed to describe 
The energy required to establish the 


the dielectric properties of a crystal. 
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tric field is half the product of the component of the induction in the 
tion of the field and the electric field. This is, therefore: 


a EiDu + E:Dx + E3Dx + F,E3(Dx3 + Dy) + E3Ei(Dy + Dis) + 
FE, E2( Dye + Day) 


Considering then a condenser made from a unit cube of crystal, the charge 


; D and the energy content is W’. If there is no leakage loss, the charge 
, : ow ou 

it can be drawn from the condenser is D = oE Whence D, = ‘E = 
of Of, 


DuE,y + 4(Diz + Da) E2 + 3 (Dis + Du) F3. If, therefore, the induction 
is derivable from a potential, Dy, = 4 (Dy + Da) or Dig = Dy. Similarly 
Dy; = Ds, and Dey = Dy. By a proper choice of axes the remaining six 
D’s can be reduced to three. In the case of isotropic dielectrics Dj, = De: = 
Ds; and 4rD\, corresponds to k, the dielectric constant. 


SECTION 3 
THE SYMMETRY OF CRYSTALS 


If a crystal has certain sorts of symmetry the number of constants re- 
quired to describe each property is materially reduced. For this reason 
we now turn our attention to a study of symmetry. 

In general, plotting a vector property of the medium for a crystal gives a 
complicated surface which we shall call a property surface. Each property 
surface of a homogeneous isotropic medium is a sphere. 

Because of the orderly arrangement of matter in a crystal, the property 
surfaces of crystalline media are commonly symmetrical. If a casting of a 
property surface were made it might fit into its mold in several positions. 
A property surface for quartz for example, if lifted from its mold and rotated 
through a third of a turn about the proper axis, would fit back into the mold. 
That is, quartz has a three fold axis. The natural requirement that mole- 
cules be laid down in a way economical of space limits the kinds of symmetry 
possible for crystals to axes of two fold (binary) symmetry, of three fold 
(trigonal), of four fold and of six fold symmetry, planes of reflection symme- 
try and combinations of axis-reflection symmetry, besides a simple sym- 
metry through a center. From these elements it is possible to divide all 
possible property surfaces into 32 classes. No other classes built from these 
elements could be self-consistent. 

A diagram study will prove this point. Ona sphere let us mark axes of 
two fold symmetry by means of a solid boat shaped figure, three fold with a 
solid triangle, four fold with a square, six fold with a hexagon, planes of 
symmetry with a solid line (great circle) and combination axis reflection, 


by means of similar hollow figures. Finally, we shall project the sphere 
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and markings onto a plane through the center. Figures 1 to 32 is a 

such diagrams. Fig. 23 for instance showsa six fold axis. Fig. 1 repre 
a medium with no symmetry whatever. The cross represents a ty 
vector property, the vector piercing the sphere above the projection s 
If the vector pierced below the sheet it would be marked with a circle. 

dashed circle of Fig. 23 indicates the boundary of the sphere without 
plying it to be a plane of symmetry. The presence of six fold symn 
requires the typical vector to be shown in six places. If an axis of two 


J 


symmetry is added at right angles to the six fold axis, it must appea: 
times and the typical vector must now appear twelve times, six times al 
and six times below the projection sheet. Continuing in this way we s 
find the self-consistent classes of symmetry to be the 32 shown in the 
grams. Often the symmetry of a crystal class is expressed by means o' 
formula. A center of symmetry is symbolized by the letter C, a binar 
axis by Ag, a trigonal axis by A3, a ternary axis by Aq, a six fold axis by 
a plane of reflection by P, and a combination rotation reflection by th: 
combination symbols &s or As. In this way the symmetry formula o/ 
quartz for example, is 342° A3. 
SECTION 4 
MATRIX ALGEBRA 

In the solution of problems of crystal physics we are involved in thi 
handling of many sets of linear simultaneous equations. As the matrix 
algebra lessens the work involved in handling sets of linear simultaneous 
equations we turn now to a study of matrix algebra. 

Several independent variables x, x2... 2, are linearly related to several 
other independent variables y,, yo. . . ¥m as 


VW = Ay} + Q12X2 + oe + AjnXn 


Vo = AX) a. oun 


CODE CHEHOKD EH 4EC O66 Ee D-8 


Vm = Gmi%1 + one AmnXn 
or briefly 
n 
ym = D>, 04%; Pee eS SES cite OR) 
j=] 


In most all such equations as (4.1) the variable to be summed over appears 
twice in the subscripts of one side. As a convention we agree to omit thi 
summadion sign and sum wherever subscripts are repeated. 

Thus: y; = a;;x; is to be summed over j 
again, if x; = bj,%, the z’s being a third set of variables we have: 


Yi = @;;bj2, to be summed over j and &. 
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e can think of this as a special multiplication of hyperquantities a, } 
If we define 


Ck = ai jb 5 ‘ 4.2 
may go from the y’s to the z’s directly thru y; = cus. We can now 
nsider the ‘“‘table’’ 

/QyjQ\2... 7 

Ghee 

\Gmi - - - Amn / 


as being the quantity a, and the table 


bi dy2 eeeee bin 


as the quantity 6. 

These ‘‘tables”’ are called matrices. 

Going to eq. (4.2) we see that the quantity c is to be a “table,” the typical 
element ¢,; of which is to be gotten by multiplying the ith row of a by the 
ith column of b, term by term thus: 


Cy = aind;; + dibs ; + tee 


After a little practice it becomes almost automatic to form the ith term 
of the product of two matrices by letting the index finger of the left hand 
follow across the ith row of the left matrix while the right index finger 
follows down the jth column of the right matrix. The fingers step along in 
synchronism and at each pause the quantities under the two fingers are 
multiplied and the product added algebraically to the accumulated sum. 

The algebra of these special multiplications is not commutable, i.e. 
ab # ba. 

Eq. (4.1) can be considered as a special case of eq. (4.2), in which the 
matrices x and y have one column only. In _ this manner a vector with 


1 
components x; x2 “3 can be considered as the matrix | x2 


X3 


If eq. (4.1) has the same number of «’s as y’s we may solve (by means of 


determinants) for the x’s in terms of the y’s. We would then get a new set 
of equations 


al 
Xi = Aij Yj. 
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owe ° fe P —l]. . 1: 
The significance of the a is that the matrix product of aanda isa: 


with ones on the major diagonal and all other terms zero. Whenever { I 


product of two square matrices gives such a matrix (known as the ide: 
factor, J) they are said to be reciprocal. Only square matrices hay, 
reciprocals. Multiplying any matrix by the idemfactor leaves the matr 
unchanged. We might consider, as part of our mathematical short han 
that eq. (4.1) was solved for x by multiplying through by a’, as 


-_ —] 
@y=e ec=2tlz= «x. 


We must remember that the order must not be disturbed as the quantities 
are not commutable, and that only square matrices have reciprocals 

The major diagonal of a square matrix is the set of terms running diagon 
ally from the upper left to the lower right. 


A symmetrical matrix has any term M;; = M;; 
An anti-symmetric or skew symmetric matrix has any term M;; = — MV 
for 1 # }. 


Rotation Theory 


The matrix algebra can be used to express a vector as a function of another 
vector, that is to handle such relations as exist between £ and P of section 2 

There is another important aspect of matrix multiplication, that of trans- 
forming a function from one set of axes to another. Let us assume that th 
new set of unit axes, x Xo and v3 are merely the old ones rotated throug! 
angle about some axis A which is a unit vector passing through the origin 
From Fig. 33 we see that in the expression: 


, , , 
X= 21% + dg\X2 oe a31%X%3 


the a;;’s are the cosines of the angles between x; and the three quantities 
x;. Conversely they are the cosines of the angles between the x,’s and } 
Consequently, if the primed unit vectors are given in terms of the unprimed 
ones by the three equations 


’ 
XL, = AgjX;j 


then the unprimed x’s are given in terms of the primed ones by the ex- 
pression: 


, 
Xp = AziX; 


This reversible relationship is well depicted by the table: 


xy Xg X3 
, 
x1 ay a2 33 
/ 
X2 Q21 a2 A23 


a32 
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[n this direction cosine table we can “‘look up” the components of any unit 
vector in terms of the other system. 
The matrix a;; is merely the matrix a,;; with rows and columns inter- 
ianged. aj; is called the conjugate of a,;. We shall denote the conjugate 
of any matrix M by M.. 
Obviously V’ is the vector sum of the 3 components (on the new system) 
of each of its 3 components on the old system. 
—ly, —ly, —ly; 
4 Vi + diz Ve + a3 V3 
‘ i os 
(V)new = (OI V fb vere reeeeees 
x5 
Fig. 33—The direction cosines of X; on X’; X’z X’; 
If the expression giving the components of V on the new system is de- 
noted by V’ we may write 
V'’=aV 
, =} o- 
conversely V=a V' 
Since x; is of unit length, the sum of the squares of its three components 
(on the primed system) is unity. 
That is ai, + aio + a3 = 1 
wer 2 2 2 
similarly a2; + do + a23 = 1 (4.3) 
and a31 + azo + a33 = 1 


Now a, can be considered as a rotation similar to the rotation a. Con- 
sequently their product aa, is a similar rotation. Let us consider this 
product. 
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The squares of its terms must sum to zero, row by row as in (4.3) 


411412033 23142143) 1 
221099093 Q}2092032 | = son 
431432033 213023033 ae 


Because of the relations aj; + ai2 + ai3 = 1, etc., we see that the ter: 
the third matrix are zero for all terms not on the major diagonal. T! 
fore, aa, is an idemfactor and the reciprocal matrix of a is the same a 
conjugate matrix. 


—] 
ag =a 


Se. ° . r ee 
Also x; is of unit length, and the sum of the squares of its component: 
the unprimed system is unity. Thus we find: 


ai + dy + ay, = 1 
aj2 +. a + 32 =]. ' } 5 
@33 + d23 + a33 = 1 


We now introduce from vector analysis the concept of the scalar prod 
The scalar product of two vectors u and 2 is u,v. It is the product of + 
lengths of the two vectors and the cosine of the angle between them. 

If we take the scalar product of x, and x, as expressed in the prime: 
system we have, since they are mutually perpendicular: 

a21 
(dit, G12, 13) dz | = Aude + di2d22 + dy3d23 = 0 
a3 


Similarly multiplying «2 and x3 scalarly, and x3 and a; we find: 


44142, + Ay2d22 + 413023 = 0 
21031 + 990392 + 493033 = . (4.5 
43101, + 32012 + d33d13 = O 


» ° U , e 
If we multiply x; and x2 etc. as expressed on the unprimed system we get 


the relations: 
41032 + 2129 + 431432 = O 
212013 + A22d23 + d32033 = 0 4.7 
213011 + de3de1 + 3303, = O 
The vector product of two vectors u and v requires the defining oi 
speciai matrix, the cross matrix. 
. 0 = Ue 
ux = 




















al 
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\Ve note that this is formed by writing zeros on the major diagonal, then 


ng back from the lower right corner writing 4, #2 and “3 around the edges. 
then make the lower left term negative, then operate on the opposite 

ie of the major diagonal so as to make the matrix skew symmetric. 

lhe reciprocal of any matrix m is 


M;(—1)**3 


- 
—— m 


aceires aro) 


ere M;; is the ji minor of | m 


he cross matrix has no reciprocal as for it (4.9) becomes indeterminate. 
Since the vector product of two vectors u and v is another vector per- 
pendicular to both w and 2 and of a length uv sin (uv) we may write, on the 
primed system 

v, X x = 2X3 in the form 


UD dy a2 a21 — 013022 + dy2023 a31 
33 0 — 211 22 = @)302, — 03 = 132 
— a\2 ayy 0 dog — Q)2d9} + 1022 33 


Matrices including vectors are equal only when their corresponding terms 
are equal. Hence, we get the relations 


431 = 2023 — 413022 
Gag = Giada — Girdag .. .- sec ccecee 105 


Q33 = 31022 — 12021 


Similarly we get the relations: 


Qi, = 22033 — 23432 
dy = e343, — 21033 
a13 = 1032 — 22031 
opie (4.11) 
Qa, = 39013 — 32033 
Qaq = 433011 — 31013 


423 = 431:0;2 — 43201 


The 21 relations between the a;;’s allow us to complete the matrix given 
four terms. 


Several Useful Matrix Relations 
0 ‘0 xy 
The del operator is the pseudo vector\Y = | 0/0 x» 


0/0 X3 
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It transforms upon a rotation of axes as does an ordinary vector: 


grad u = u,, a matrix 4.14 
div u =. 4u,a scalar 1.15 
curl «w =V X u, a matrix 11 
grad radius vector = Vp = J, the idemfactor 1 
Gke.. J) a 1.60% C 118 
(abc...) = . 5. 6Dae ; 1.19 

len @ @O5.45° t/a, 0 O 
O tte O46 cs = QO 1/az2 O _ 
oS me ciscx © © tf... - 


, ° —] 
: ; (matrix) 
~ l c / 
(Scalar times matrix) © = ; cealveieye ROE 
scalar 


SECTION 5 
THE GEOMETRY OF ROTATIONS 
As a first application of the matrix algebra let us compute the a matrix 
for a few general rotations. Although we can consider a general rotatio. 


as one of angle ¢ about the unit vector* A, it is easier to consider a genera 
rotation as three successive rotations about coordinate axes. } 





A study of Fig. 34 shows that for a counterclockwise rotation ¢ about . 
the new components of a vector V are: 


w=; 
V> = Vo cos @ + V3 sin 07) 
V3 = —Ve2sin ¢+ V3cos¢ ] 
whence V’ = aV where 
‘1 0 0 
a={0 cos@ sing PR eae ee oe 
0 -—sing cos@¢ \ 


* A general rotation of amount ¢ about the unit axis A is given by 
a = AA, + (I — AA.) Cos¢ + Sing A 


See Vector Analysis (Gibbs Wilson, Yale Press) pp. 338. 
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Similarly, for a counterclockwise rotation ¢ about x2, we have 


cos@ 0 —sing 
a=1{0 1 0 
sing 0 cos @ 


and for a counterclockwise rotation ¢ about x3: 


cos@ sing 
a= |-—sing cos@ 
0 O 





0 
0 
1 





+t 
w 


t Xo 





Fig. 34—The relationship between the components of a vector on one coordinate 


system and on another. 


In the appendix we give the special transformations corresponding to the 


symmetry operations of the 32 crystal classes. 


rotations: 
e = 46 
o" = a’ x! 
Z'" = gx" 
the resultant rotation is 
a!" = e"a'ax 
or o’" =e Re 





where R= a"'a'a 


If we have three successive 
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The I.R.E. Orientation Angles and the I.R.E. Matrix 


The Institute of Radio Engineers has proposed that, for quartz crys 
all orientations be given in terms of three rotations ¢, 9, y about 
and x3. respectively, starting with the plate length along x; width ; 
x2 and thickness along x3. (Here x3 is the z or optic axis, x; is the el 
axis. ) 

Whence, here: 


cosy siny 0 cos@é 0 —sin@ cos@ sing 
R=1|{-siny cosy 0 0 1 0 —sing cos¢ 0) 
0 0 1 sin@ 0 cos 6 0 0 1/ 


and carrying out the two matrix multiplications: 


xy Xe x3 
( . ° ur 
| cos¢@cos@ cosy sing@cos@ cosy —sin 6 cosy] x 
—sin ¢ sin y +cos ¢siny 
—cos¢@cos@ siny —sin@cos#@ siny (5. 
| ° ° ° "rt 
| —sin @ cosy +cos ¢ cosy sin sin | x2 
. . . “7 
cos ¢ sin 8 sin @ sin 8 cos @ |} x3 








If we denote the unit vectors along the length, width and thickness as PP, 
and P; respectively we have as a matrix defining the plate: 


| dr Pete an cn Lind cots 


The I.R.E. orientation system is useful to the designer of crystal plates 
because his problem is to choose such values of ¢, 0, y as to give the plate 
certain physical properties along its length, width and thickness. The man 
who cuts the plate has a different problem, that of moving the crystal (and 
hence the x; x2 3 axes) about a fixed saw so that the plate cut parallel to the 
saw blade is what the designer ordered. 

Let us consider such a system as shown in Figs. 37, 38 and 39. In Fig 
38 the crystal stands with its optic axis along P3, its + electric axis (for 
right hand quartz) along P. Since the shop man considers clockwise rota 
tion as positive we now rotate the crystal through angle U’; about P3 clock 
wise, we then turn the crystal through angle U2 clockwise about P;, and 
finally, after cutting out a slab of required thickness, we turn it clockwis« 
through angle U; about P; to cut its length and width. 

On the plate axes P; the crystallographic axes +; x2 x3 are now given by 
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x3 
1/ 
X2 
x) 
{ Fig. 35—The initial position 0, 0, O) for the I.R.E. direction angles. 
= 
P; 
es 
te 
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xy 


The final position (#, 0, ¥) for the I.R.E. direction angles. 





lig. 36 
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sin U0 1 6 6 cos U; sin U;0 
cos U, 0 0 cos Ugsin U2 —sin U3 cos U;0 
0 1 0 —sin U2 cos U2, 0 O 1 


{ , , , ° y . , . > \ 
| cos U,; cos U; cos U, sin LU’; sin U; sin U2 
—sin U,; cos Up sin Us; +sin U; cos U2 cos Us 


—sin U; cos U; —sin U, sin U3 cos U; sin Us! (5.9 
—cos U, cos Uz sin Us; +cos U; cos U2 cos Us 


sin U2 sin Us —sin U. cos U3 cos Us 












































Fig. 37 Fig. 38 


Fig. 37—The (0, 0, 0) position of a shop system of direction angles. 
Fig. 38—The second position of a shop system. 


From (5.8) we see that P = rx and hence, if this is to be the same plate 
the designer specified by P = Rx we must have R = r whence we may 
equate the terms of (5.6) and (5.9) to get the relations 
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cos U. = cos @or Uz = + 8 
tanU, = coty. 
tan Us; = cot¢ 


U,=y—-9 + nr 
Us = 6 
U3; = ¢@—-90+ 

















Fig. 39—Cutting the slab and trimming it to the piezoid boundaries 
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Where » is any integer positive or negative, including zero. If we ta 


U,=y—90 
U.= 06 


Us=¢+% 


The matrices are consistent term by term. 


SECTION 6 
CRYSTALLINE DIELECTRICS 


As a first application of the matrix algebra considered as a linear vector 
function let us reconsider the problem of the crystal in an electric field. 

The relations of chapter II, equation (1) can be written in the abbreviated 
form: 


D = DE where D,, = Dey 


in accordance with the system of abbreviations adopted in the appendix 
If we put 
4aDre = Rre 
equation (1) can be written 


1 
Vis eG Ak ava ive bok wee MOD) 
An oi 


In order to investigate the effects of crystal symmetry in determining 
the least number of dielectric constants that are required for a given class of 
symmetry it is desirable to find the electric induction D for any system oi 
axes, Suppose that we choose a system for axes x; , x2, Xj related to x), 
x2 , x3 through the relations: 


, , , 
= A%X1 + Ay2X2 + a43X3 
’ ’ ’ ’ 
= A9%] + d22X2 + Ce he dae ks ‘ . (6.2 


, , , 
= dgiX1 + AgXq + 33% 


. . , . . e 
where a); is the cosine of the angle between x; and 2; , diz is the cosine of the 


/ 
angle between x; and x2 etc. 
Equation (6.2) can be abbreviated to 


x’ = ax 


where a is the matrix 
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Qy1 G2 213 
a = | 421 G22 doz 
31 A32 233 


V; 

It is shown in the sec. 4 that any vector V = | V2 }can be written on the 
V3 

new system of axes as V’ where V’ = aV, conversely V = a ‘V’; a ‘ is the 


matrix reciprocal to a. Since the induction D and the electric field E are 
simple vector functions they transform as the vector V, that is: 


D mal .... (6.3) 
E’'=aE....... (6.4) 
But by (6.1) 
1 
D = — kE 
4a 
whence: 
aD = — aka.aE 
or 
D’ i k'E’ (6.5 
= 4 (O 
4r 
if 
k’ = aka, (6.6) 


We see that the form of (6.5) is the same as that of (6.1) for any set of 
axes if (6.6) is used to define the new dielectric matrix k. 

To apply this relation (6.6) to a particular crystal let us consider a tetrag- 
onal crystal (which has its properties unchanged by a rotation of 90° 
about a four fold axis). Let us choose the four fold axis as x3 and then rotate 
the axis 90° about x3. In this case 


uv =-10 010 
a={[1 00 )and the reciprocal matrix a’ = {| —100 
0 O01 001 


whence equation (6.6) becomes: 
010 ky Rio k3y 0-10 Roo — kp hog 
k’ =(-100] | Ri kokn (1 OO) =| —ke kun —ka 
00 1 k3y hog ka 0 01 hog — ks; 33 
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But because of the symmetry & is unchanged by this transforma 
hence 


k=k' 
Two matrices can be equal only if corresponding terms are equal, henc: 
ky = ko, ky = —kip = 0, hes = —kg and kgy = hog 
whence 
hog = key = 0. 


We are left then, with the dielectric constant matrix for the tetragona| 
bisphenoidal class: 


kn 0 0 
k = 0 ky 0 
0 0 kgs 


Applying other transformations possible for tetragonal crystals gives no 
further simplification. 

If we go through all the symmetry transformations possible for the 32 
classes we find that cubic crystals require but one dielectric constant, hex- 
agonal, trigonal and tetragonal crystals require two constants, orthorhombi 
monoclinic and triclinic crystals require 3. 

As the triclinic class has no fixed axes or planes of symmetry the reduction 
of its 6 constants to 3 is not so obvious. It may be seen by expanding into 
ordinary xyz coordinates, that p.kp = 1 is the equation of an ellipsoid, (p 
is the radius vector) where the six k’s are the coefficients of x”, y’, 2°, ys, 
zx and xy respectively. If we choose the coordinate axes along the axes of 
the ellipsoid the yz, x and xy terms drop out and only three ’s are needed 
With triclinic crystals then, if we determine the axes of the ellipsoid, then 
choose the coordinate axes along them, only three dielectric constants are 
needed to completely specify the polarization in terms of the electric field. 
The determination of the ellipsoid axes must be made experimentally as 
there are no symmetry elements to guide us. It is possible to compute the 
positions of the axes from the 6 k’s by solving a cubic equation. 

The values of the k’s depend on the frequency of the applied field. In 
crystals of low symmetry the ellipsoid axes for different frequencies do not 
necessarily coincide. 


; : R = - 
Another vector quantity of interest is the polarization, P = D — i E. 
v 


vite ee 
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Using (1) this becomes 


P= bs (k —I)E ... (6.7) 
4a 


SECTION 7 
QUADRATIC FoRMS 


Often the elements of a matrix are themselves functions of other quanti- 
ties. In order to relate the elements of one matrix with those of another by 
means of a matrix multiplication, we may make a single column matrix of 
each of them. We then wish to know how a transformation of axes changes 
the elements of this single column matrix. Consider a symmetrical matrix 
b that relates two vectors u and 0: 


u = br. 


A transformation of axes, a, changes wu andv tow’ andv’. Multiplying « = 
bv through by the prefactor a we have 


au = abv. 
» . . . ls 
We now replace au by its equivalent a v’ whence: 
, “59 
u’ = aba v 
so that 
u’ = b'r’ 
if we define 6’ as 
—] - 
b’ = aba (7.1) 


To be in accord with common usage we now rearrange 6 according to the 
arbitrary scheme: 





bi By 
bo» Be 
bs3 = B; 
bog 

bai ; 
bie Be 











We wish to know what operation to perform on B to get B’ corresponding 
too’. If we expand b’ = aba ’ it is easily seen that b’ = aB where 
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( ait aie ai3 2a12 diz 243411 2an a12 | 
an Sm a33 2 doe d23 2423 G21 2a) a2 | 
| a3 ase a33 2432 G33 233 @31 2a31 G32 | 
Lili 22432 23 A33 a cee re | 
a= + 423432 +2343: 9+ 421432} 2 

} 
431411 3232 33. Q13 oe —s men 

| + 413432 +and33 +4145: | 
| ; Q12 223 13 G21 yd | 
om Se ae +1302 + 411d23 Stealted| 


Because we shall often need to form the a matrix from the @ matrix we 
need an easily remembered mechanism for doing so. We notice that thos 
are four kinds of terms in the @ matrix and that the four kinds can be sepa 
rated from each other by two center lines, one horizontal, one vertical. 
This gives us four squares of nine terms each and we can correlate each term 
of any square to a term of the a matrix by means of its position in the squar 
The terms of the upper left square are the squares of the corresponding terms 
of the a matrix. To form any term of the lower left square we cover the 
corresponding term of the a matrix with our finger and multiply the visibli 
terms of that column. To form any term of the upper right square we cover 
the corresponding a term and write down double the product of the visible 
terms of that row. To form any term of the lower right square we find the 
corresponding a term, strike out that row and column and write down thi 
sum of the remaining cross products. A study of the following diagram 
will help to remember these rules. 





Terms are squares 
of corresponding 
a terms 


Omission products 
doubled 


Sum of omission 


Omission products pane 
cross products 











Fig. 40 


SECTION 8 
CRYSTAL ELASTICITY 
Stress 


Consider a point P in a medium acted on by forces. If a small area is 
chosen about P the medium on one side of the area exerts a force on th 
medium on the other side. The force will depend on the size of the area and 
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the direction of its normal n. We shall choose a triangular area ds such 
.t an arbitrarily chosen set of mutually perpendicular unit axes x , x, 
pass through the vertices of the triangle. Let us consider the conditions 
equilibrium of the tetrahedral element of volume so formed. The areas 
rmal to x1, %2, %3 are ds, , ds2 , ds3 , respectively, and the forces per unit 


rea acting through these faces are: 


fu far fou 
F, =| fr], F, = | fo |, F; = | foe 
his fos fsa, 


\ny body forces (such as gravity) depend on a higher order of smallness 
that is on the volume rather than on the area) and hence are negligible. 


Whence for equilibrium: 
Fds = F,ds, a Fidsy a F3dsz 


But 
ds, = mds, dso = nods and dss; = n3ds 


ny 
where | #2 | is the normal to the area ds. Whence we may write: F = fn 
Ns 


where f is the matrix 
‘fu fie fis 
far foo fos 
fan fae fas 


For the body to be in rotational equilibrium the tangential forces must 


balance, hence fiz = fai , fis = far and fos = fee . 
Transformation of Axes 
A change of axes that transforms vectors through F’ = aF changes F = 
" 1 il ae ‘ ‘ 
fn to aF = afa ‘an so that if f’ = afa ‘ then F’ = f'n’. 

In order to relate the stress to other quantities through a matrix we wish 
to convert it into a single column matrix. We put fu = X;, fe = Xo, 
fs = X3, fos = fro = Xa, for = fis = Xs and fiz = for = Xe. 

Changing to the X representation we find 


xz’ = aX’ (8.1) 


where @ is the matrix eq. (7.2). 








24 BELL SYSTEM TECHNICAL JOURNAL 


Strain Theory 


pi 
If the dimensions of a body change, a point p = | p2 | is moved to p + ¢, 
3 
o1 
where ¢, = | o2}. A neighboring point p + wu is moved by an amo 
93 


Tp+u given by gpiu = (Yor) +o,. The movement of p + wu relative to 
p is CG = O pt u ~ op = (Voce) cu. 
The 9 components of (7a). describe the sort of movement in the neighbor- 


hood of a point; they are the strain coefficients. If the strain matrix js 


«= (Vo-.).,a transformation x” = ax causes this to become aea, = (aVo.u 
‘ we 5 , a 
and if a7 =’ and ao = o’ so that o.a, = o, we have & = (V’a¢)¢ if 
= aea,.... ie kta eens (8.2 


When we arrange € as a single column matrix e we shall, following custom, 


Ooo , Oas 


take e, = ,@; = etc. This has the effect of moving the 2’s of the a 


Ox3 Ox3 
matrix to the conjugate position so that, while x transforms as x’ = ax, 
e transforms as e’ = a, €. 

We shall take tensions as positive stress elements, and elongations as 
positive strain elements. The shear strain, e, = (0, 0, 0, 0, 0 e) becomes 
& 


° Py. é rr. ° 
upon rotating through 45° about 1; , e. = (3, — = 10,0, 0,0}. This shows 


that to be consistent, a positive shear strain about «3; must mean an expan- 
sion along the line x; = x2 and an equal contraction along the line x, = —.» 
A positive shear stress is one that tends to produce a positive shear 
strain. 
By superposing such strain elements we see that the € matrix (useful in 
displacement problems) may be formed from the e matrix (which is useful in 
stress strain relation) as 


C1, 306 365 
€e= 4&6 €2 he, . (8.3 


5e5 34 €3 


This slightly awkward relation is used solely to make the ‘‘work done in 
straining” expressible as 


2W = Xe =e,.X patina tamed ae 


If the e’s were taken as equal to the e’s the work would be: 2W = Xye, + 
Noes + Xse3 + 2X4eq + 2Xs5e5 + 2Xee5. This would be awkward in some 
later problems. 


Cutie eae aaa 
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oA he lth Ca api aa lela 
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If the scalar W is to be unaffected by a transformation @ we must have 
i)’ = e-X unaffected. If we write 


. -] - Fos 
W =ea aX =ea X’ 


we have 
aa , 4 oye 
W=W=eX 
if 
’ 2 
Cc = ca 
when 
wl fo” 
e’ = a. €; : (8.5) 


This substantiates our previous statement. 


Relation Between Stress and Strain 
If the strain in an elastic body is proportioned to the stress we may write: 
€y = SyX1 + SwXo+ +--+ SX 
éo = SuXi+ °°: awe 
Where the S’s are elastic modulii. In matrix notation: 
e= SX (8.6) 
Conversely X = S ‘e or if S*=C 
X = Ce (8.7) 
The C’s are called elastic constants to distinguish them from the modulii S. 


vu —1 —l io 1 , . | , ° 
Ase = SX,a, e =a, Sa aX, and since a, e = e’, (the representation 


of e on a new axis system related to the old one through the matrix a) and 
re . . G =] —-ln ly / r 
aX is X’, then we may write (a, e) = (a, Sa )(a@X) as: 


e’ = S'X’ where S’ =a,'Sa’ (8.8) 
Similarly operating on X = Ce we find 

X’=C'’e’ where C’ =aCa, (8.9) 

The energy required to cause the strain e is 
; . lie eae ae P 
W = | X,de, = 5 Xrep = 5 Sre X,X, (8.10) 

whence, if W is a perfect differential, 

aw aw 

Se = So (8.11) 


ai OX, OX, 7 OX, OX, 7% 
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Similarly 


Cre = 678 
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This reduces the constants and modulii to 21 of each. 
If a transformation is performed that is permitted by the symme: 


the medium the elastic modulus matrix is unaltered. 
system has a binary axis. 


The monoclin 


If we choose this as «3; and rotate the axes | 


-1 0 0 
about this by means of the matrix a = | 0 —1 OJ] we have s 
0 oO 1 
a. 'Sa’ = S. 
(1 00 0 00) (SuSpSpSuSisSe) (100 0 00} 
010 0 00) | SiSoSo5SoSe55/ 1010 0 00 
001 0 00] | SiSeSsSuSsSx} 1001 0 OO] 
(000-1 00) |SySuSSuSisSe| |000-1 00| 
1000 0-1 0} | SuSrsSsSiSisSie] |000 0-10] 
(000 0 O01) [SieSxSaSeSsoSee} (000 0 01) 


Sun Sie 
Si 9 Soo 
| Sy Ses 


|—Su —Sr 
| —Sis —Sos5 
| Sie Soe 


Siz 
S23 
S33 
S34 
S35 
S36 


—Su —Sis Sis 
—Soq —So5 Sag 
—Sx —S35 Soe 
Su Sas — St 
Sis S55 —Ss6 
—Sig —Sse Seg 





Equating terms, those whose signs differ in S and S’ must vanish. 
Proceeding in this way through the 32 crystal classes we arrive at the 
ten following matrices that cover the elastic behaviour of all 32 classes 


Triclinic System 
21 modulii 


Monoclinic System 
x3 axis binary 
13 modulii 





Su Su Su 


| Si2 
| Sis 
| Sis 
| Sis 
| Sie 


Su 


} Sie n 
| Siz Se: 


S =— | 
10 


lO 


So 
So3 
So 


S 25 


0 


2 S23 


Sog 
S33 
Ss 
Sas 
Sx 
Sis 
S33 
0 

0 


Sx Su Su 0 


Sis 
So4 


Sy 
So 
Sax | 
Sis| 
Sse 
See 


S) 5 
So5 
Ssi S35 
Sia S45 
S45 a 


Sa : 





S45 


The C matrix is 
entirely analogous 


(8.13 


The C matrix 


entirely analogous 


(8.14 















SSeS, 


ombis System 
binary 
modulii 


letragonal System 
r; a fourfold axis 
Classes 9, 10, 13) 


7 modulii 


letragonal System 
v3 a fourfold axis 

x, a twofold axis 
Classes 11, 12, 
14, 15) 

6 modulii 


Trigonal System 
x3 trigonal axis 
(Classes 16, 17) 


7 modulii ; 
> = 


Trigonal System 
x3 trigonal axis 

x, binary (Classes 
18, 20, 21) Ss 
6 modulii 

(alpha quartz) 


Hexagonal System 
x3 a sixfold axis 

x, a twofold axis 
(Classes 19, 22, 
23, 24, 25, 26, 27) 
5 modulii 
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Su Sz S30 0 0 ) 
Siz Sx Sx 0 0 0 | 
= | Sas Soz S330 O 0 The C matrix is 
' 10 0 O Syd O entirely analogous 
YY ¢ & 8 Sa O | 
0000 0 Se| (8.15) 
(Su SeSs0 0 Sy) 
Ss SnSnO0 0 —Sy 
1Sis S33 S33 O O 0 The C matrix is 
10 YY 0 Se 0 0 | entirely analogous 
0 00 0 Su 0 | 
[Si —Se0 0 0 Ses} (8.16) 
(Sn Sze 530 0 0 ) 
Su Sa Sa O 0 O | 
om Sis Si3 S330 O O The C matrix is 
0 0 0 SyO 0 entirely analogous 
0 0 0 0 SuO | 
(0 0 0 0 O Se) (8.17) 
) The C matrix 
Su Si Sis Sig — St 0 |is analogous 

Sy Syn Sis —Sig Sos O | except that 
Siz Sig S33, O 0 0 | Cae = Crs 

| Su—SuO0 Su O 2Se \Cos = Crs 

1—Sey Ss0 O Su 2Su 1Ces = 3 

1 O 0 0 25% 2S 2(Su —Siz)| (Cu —Crz) 

( (8.18) 
(Sir Se Si3 Su 0 O ) The C matrix 
1Sio «=3Sy Sis —Suy O O | isanalogous 
\Su- Su Ss 0 O 0 | except that 

= |Su —Su 0 Su 0 0 | Ce = Cu 
0 0 0 O Su 2Su Ce = } 
| 0 0 O QO 2S 20S —Si2)| (Cu —Ce) 

( (8.19) 


Si Si Sis 0 0 0 
|S Su S30 0 O 
1 Sis Sis Sn 0 0 O 


i0 O O 


| 


Sy O 


10 00 0 Sy 


10 0 0 0 0 2(Su —Sy)} 








The C matrix 
is analogous 


except that 


Cos = 3 
(Cu — Ci) 
(8.20) 
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Cubic System (Su Sx S20 0 0 
X1, %2 and x3 1Se Sy Sep O O O The C mat: 
fourfold axes —_ |Si2 Se Su 0 O O entirely anal 


) 
| 
| 
3 modulii ~ 10 0 0 Sud 0 | 
10 000 Sy 
(0 0 0 0 0 Sul 


Q 1 
«1 


Isotropic bodies Su Se Se 0 0 0 1 The C matri 
2 modulii 1 Sie Su Sp OO O O | analogous except 
S= |Si2 Se Su 0 0 O that 
0 0 0 So a | C2 = A (Cu Tt 2 
10 000 & 0 
ee a ee | 


Se = 2 (Su — Si) 


Several Elastic Ratios in common use are given here for reference: 
Young’s Modulus: A tension stress X divided by the component of 


strain in the direction of X , ¥; = —. If the coordinate axes are chosen so 
ej 
that the stress lies along X, , Y; = Sn" To find the value of Y in an arbi- 
ll 


trary direction, (6, ¢) find S’ for a transformation that puts X’ in the di- 
rection (6, ¢) 


. aah gd 
S’ =a, Sa 


Where a is taken as form (21.4). Whence we obtain: 


1 4 4o 44a 4 222 222 22 4a 
(5 — } = €y8aSi + $182S22 + €2S33 + 515202844 + 616252555 + C151 52S 66 
9.9 
2 30 . * 22 3.4 
+- 2¢1510259S 56 -b 26181 CoS2S 46 a 201510259556 4- 2615159596 
i. de 3.3 4 ¥% 2220 
- 2¢151595 16 +t 2¢1520)S 15 + 2015189625 14 + 26162525 13 
22 he 222 3 22 
+ 2cisisoSiz + 2510252523 + 251S2C2S34 + 2¢1516252S45 
3 « a YY ae 12) 
+ 2610252535 ot 261512595 05 + 25362595 04 Me A a a ee (8.23) 
Rigidity Modulus: The shearing stress divided by the component of shear 


about the axis of shearing stress. For shear about x; , 


| re woe} 
1= > ( ) 


Its value in another directions can be found as V4, was above. 





and 


ture 


anc 


we 


wl 
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e bulk modulus: The change in volume per unit volume for unit hydro- 


static pressure is the bulk modulus, 7. Fora stress X, = (1, 1, 1, 0, 0, 0) 
é¢ = (Su + Sie + S31 ’ Siz + S33 + So3 ; S31 + S23 + S33 ia ) 
H etete=Sut Set S33 + 2S + 2531 + 2523) (8.25) 


This is obviously independent of the choice of axes. 


The Temperature Coefficient of the Elastic Modulii and Constants 


C= C+th+fh +f + --- (8.26)* 


Se S°+th+fH' + fH +.--- (8.27) 


C° and S° denote the values of the C’s and S’s for some standard tempera- 
ture ¢ = 0) then as the transformations are 

C=aCa, and SJ = as Sa ' or 

C=alCo+th + thn + th" «+a 
and 


S’ =a, (S° + tH + fH' + fH'---)a’ 


we see that 


C=C + th + Ph ---.. (8,28) 
S’= S°4+tH' +H" ---. (8.29) 
where 
h! = aha, etc. (8.30) 
H’ = az'Ha ’ etc. (8.31) 


That is, the /’s transform as the C’s do, and the H’s transform as the S’s 
do. Consequently we may copy their respective forms from the C and S$ 
matrices for any particular crystal class. 

When the temperature coefficients of the constants or modulii are known 
in the form: 


Ci; = Ci; (1 + tT ¢;;) (8.32) 
Si = Si AUT 


* The n of é” denotes the nth power of the scalar ¢; the » of hk” is merely another matrix, 
it does not mean a power. 


) (8.33) 


Si; 
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we may write: 
ee & oe ee 
Hi; = Si; Ts;; 
Multiplying (5) by (6) we get: 
SC =I = S°C° = SCC? + U(S°h + HC?) + F (SH + Hh + WC?) 4 


whence, for this relation to hold for all values of ¢: 


gt ee 8 3 
whence 

h = —C°HC 

H = —S°hS° ws 
also 

Sth + Hh+H'C°’=0............. (8.38 
so that 
hi = hS°h — C°H'C° 


H' ee. ae et 


From these we can compute the h’s given the H’s and vice versa. 


SECTION 9 
TEMPERATURE EXPANSION 


The change in the dimensions of a crystal caused by a temperature change 
can be considered as a strain. The shift of the terminus of a vector | 
relative to its origin is given from the strain matrix e by the equation A/ = 
tel 

Since ¢ is symmetric a proper choice of axes makes it possible to make the 








strain per degree a diagonal matrix, 


A, 0 0 
Alw=j4l whee Ami U de. Oeeii.ccrs.- (9.1 
0 0 Ag 


As | and Al both transform as vectors, a transformation @ causes A to 
transform as 





0 | 
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The strain can easily be extended to a function of ¢ and t’ as follows: 
Al = tAl + #Bl 


Applying the prefactor a to both sides and putting the idemfactor in 
between A and / and between B and / in the form J = a ‘a we have: 


aAl = t(aAa 


Al’ 


A’ 


TEMPERATURE VARIATION OF THE ISOTHERMAL ELAstTic MopuLu 
STRESS VARIATION OF THE TEMPERATURE EXPANSION COEFFICIENTS 


We can write the isothermal elastic modulus matrix at temperature 


6+ tas 


and the coefficient of temperature expansion at constant stress Y as 


Let us take a unit cube of crystal about the cycle indicated in the table; 
starting with the cube in the unstressed unstrained state at absolute tem- 


perature @: 


Operation 


Heat at zero stress 
Apply X isothermally 
Cool at Const. X 

Apply —X isothermally 


If we sum the strain changes in this cycle to zero we have 


: 
i e elongation per unit length per degree in the direction (6, ®) is 
cos @ sin ¢, sin 8 sin yg, cos ¢ A, 0 O 
0 A. O 
0 O As 
whence 


a2 9 ah SP 2 
= A,cos 0sin y+ Azgsin 6sin ¢+ Ascos ¢ 


Nal + ((aBa')al or 


where 


t Al + PB 


1 > 
aAa 3’ = aBa 


SECTION 10 


S'= S° + tH 


A = A°+ LX 


Dg nn a Change in Strain 
0 tA° 
X (S*° + tH)X 
0 —t(A° + LX) 
—X —S°X 


H=L 
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so that we may write 
S=S +4H... vi 10.1 
PS we | eee 10,3 


This tells us that we may determine the temperature coefficients 
elastic modulii by measuring the effect of stress on the temperature , 
pansion coefficients. 

In a similar way we find that if the isothermal elastic constant matrix »: 
temp. 6 + ¢ is: 


oe elec viene: 10.4 
then the relation between temperature and stress at constant strain ¢ j; 
2) | ee ae eee an 10.3 

where 


B= B’+he...... 


The Difference between the Specific Heats at Constant Stress and Constani 


Strain 
Writing for the specific heats at constant stress and at constant strai! 
a” and o’, respectively, we can perform the following cycle: 


Change Change 


Operation in Siseen| in Strata Temperature Work In Heat Out Entropy 
° — pla! 
Heat at zero stress. — tA @tod+t 0 —pta? Peay, 
+ t/2 
Restore zero strain 
A r area ) 
isothermally ... _ —tA° 6+ 5 At C A® Q P v , 
. . . plo* 
Cool at zero strain . _ 0 6+ ttod 0 pla” ; 
ae ae 
| 


Equating the sum of the entropy changes to zero: 


ae "oe 
C+ oa pio” — o°) 


Equating the work in to the heat out: 


: 


(o? —o") = 


.) 





(nt! ll nae . (10.7) 





Su 


V 
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hermal and Adiabatic Elastic Modulii 


pint Sie acs, oneal 


10.1} % Let us take a unit crystal cube at temperature @, apply any stress X 


I 


Bad ibatically, heat it to bring the temperature back to @ at constant stress 


10.3 . ra ; : 
en release it isothermally. The cycle is analyzed in the table: 
re ¢ Ps Z ee ee : Work Heat Ou Entropy 
eration | $83 Change in Strain emperature ork In eat Out Chinkee 
at mee 
adiabati- 
cally X SeX 6to@—t 4X SX 0 0 
104 ; 
Heat to @at ' , —tpa? 
5 const. X 0 t(A + HX) 6— tto@ tX (A° + HX) —tpa” 6 1/2 
11N ¢ sat S 
Remove X 
7 isother- 
10,5 mally —X — SeX 6 —4X.S°X Q 0/6 
Summing the strains to zero: 
(10.6 (S*° — S*)X = t(A° + HX) 
omstant If we equate the total entropy change to zero we obtain an expression 
for Q that can be substituted in the relation ‘work in = Heat out.” This 
gives us: 
Strat Epo? 
—32 X(S° — S*)X + tX(A° + HX) =3- “. 
== and from these two expressions we derive, writing ¢ for S'° —S": 
/ 6 47° , Zo 4} , ( 
pla’ = (A° + HX)(A° + HX)... (10.8) 
- 1/2 po? 
Q which is, to the first order of the small quantities Y: 
+t 
7. 6 Zz z > 
; o-— A°A? + 2HXA?. (10.9) 
| po 


and since X = Ce we have also 


e -- ee 
(A° + 2HCe)A° .. e+.» (10.10) 
poP 
Whence we see that as the stress approaches zero as a limit ¢ approaches 
° 6 z z « 2 to a ° ° 
¢ 2. A° A?. If we write similarly C” — C* = y we have multiplying 
po 


S° = $+ oby C° = C’ + y and dropping higher orders of small quan- 
10.7 tities: 
v= -—C°’¢C” (10.11) 
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For example we have for quartz at 20° Centigrade 


(14.4 
14.4 
: 293 x10” 1.8 


5° = 0 |(14.4 14.4 7.8 0 —0 0) 
o 2.65 X 7.37 X 10° ( 





DEE CN 








3.13 3.13 1.69 000 
5:13: 3.43 1.209 000 
Pe 1.69 1.69 .907000 1078 
0 0 o 800 
0 0 0 000 


0 0 0 000 


and as 








1298 —166 —152 —431 0 0 
(—166 1298 -152 431 0 0 
—152-152 90 0 0 0 
—431 431 0 20088 0 OO 
0 0 0 0 2005 862 fo 
0 0 0 0 862 2928 


gre me 10 15 x 





whence $n = Sif (1 — .00241) 
dr = Si3 (1 — .0189) 
du = Sig (1 — 0111) 








st 
os3 = S33 (1 — .000917) 
For Rochelle Salt we have: 
(59.9 
. 38.1 : 
= 293 X 10 44.8] (59.9 38.1 44.8 0 0 0) 
1.79 * 15.5 10° | 4 
0 
(0) . 
‘ 
38.0 24.2 28.5000 
24.2 15.4 18.1000 
me §& 2 ? - 
_ |28.5 18.1 21.3000 x 107% ! 


0 0 0 000 
0 0 GO -G:00 
0 0 0 000 

















~~ 
MERE TEIN 
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( 4690 —795 —2180 0 O 0 ) 


—795 3205 1691 0O 0 0 
—2180 1691 2815 0 0 0 |x 107 


2 vto = } 
As 5 0 0 0 6060 O 0 | 
| Oo 0 0 0 3000 0 | 
| 0 0 0 0 0 8020) 
» So that 
éu = Sii (1 — .0080) dv = Si2 (1 — .0305) 
o2 = So9 (1 <— 0050) 013 = S13 (1 —_ .0103) 
o33 = 533 (1 — .0076) go3 = S23 (1 — .0107) 
ou = Sis 
dss = Sts 
; dos = Sée 





At the temperature of maximum piezo activity the components of @ 
for Rochelle are smaller by about 33%. 
SECTION 11 
Tue Priez0-ELECTRIC EFFECT 


Some crystals develop an electric charge when subjected to mechanical 
stresses. As far as the effect is linear it may be expressed by: 
dD, = dyX, + dyXo mae 3 digX¢ 
De a dX, + we's , (11.1) 
Ds = d3X1 + ites d3gX¢ 
or in matrix notation 
D= dX (11.2) 
where the 18 constants d;; are called piezo-electric constants, and D is the 
electric induction. 
On rotating the axes by means of a transformation a, the vector D be- 


comes D’ where D’ = aD. Thestress transforms as X’ = aX whence D = 
dX becomes D’ = ada‘ X’ or D’ = d’X’ where: 


d’ = ada -:- (11.3) 
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If transformations permitted by symmetry are performed, the d 
is unchanged. Class 3 has a binary axis only, if we choose this as 
perform the transformation 





-1 00 ; 
a= 0 —1 0] we find: ‘ 
; ( 
0 01 , 3 
100 0 00 
) Te 2 60 ; 
—1] 0 0 di, dye diz dig dis dye } : 
1001 0 0 0  ( 
d= 0-10 do do dog doy dos dx f 
0 01 ‘eager are ‘(O00 -1 OO 
M31 G32 133 d34 G35 436 10 00 0-10 
(000 0 01) } 
—dy — dig — dig dis dis — dig ’ C 
d = — dy — dy — dog do4 dos — dy 
ds d32 d33 — d34 — dss 36 
For this to be consistent with the original d matrix the terms with conflicting | 
signs must vanish. ( 
Applying similar analyses to each of the 32 classes we arrive at the set ( 
matrices: 
Class 1 (asymmetric) dy, dy dig dia dis dig 
No symmetry d = do dog dog dog dos dog (11.01 
a3 d39 d3g d34 d35 6 ( 
Class 2 (triclinic pinacoidal), center of symmetry d = 0 (11.02 ( 
Class 3 (monoclinic sphenoidal 0 0 O dy dis 0 \ (sucros 
x3 is binary d={0 0 O dau dos 0 (11.03 
dz; dx d33 0 O doe 
Class 4 (monoclinic domatic) dy dye dy30 0 di 
X3 plane is plane of d= dy dx» dog 0 0 dog (11.04 
symmetry 0 O O ds dss 0 
Class 5 (monoclinic prismatic) center of symmetry, d = 0. «(11.05 
{ 
Class 6 (Orthorhombic o 8 6. d&0-6 (Rochelle 
bisphenoidal) d={00 00 dO (11.06 
%1, X2, X3 binary o.0 & 0.0 te 
Class 7 (Orthorhombic 
Pyramidal) 00 0 0 ds 0 
x3 binary, 4 and x2 d={0 0 O dy O OF (11.07 
planes of symmetry ds; dz. d33 0 O O 











» + 
sel 


1.04 


11.05 


helle 


11.06 


1.07 





i 
: 


ss 9 (Tetragonal 0 0 0 dy ds 0 
bisphenoidal) d={0 0O 0 —dy dy 0 (11.09) 
, is quaternary alternating dx —d, 0 O O ds 
Class 10 (Tetragonal 00 0 dy ds O 
pyramidal) d={0 0 O dis —dy 0 (11.10) 
3 is quaternary ds, dg, d33 0 0 0 
Class 11 (Tetragonal 00 0 dyO O 
scalenohedral) d={0 000 d&O (11.11 
’3 quaternary, x; and x» 00 0 0 O dg 
binary 
Class 12 (Tetragonal 000 dy O 0 
trapezohedral) d={|0 0 0 0 —dy0 (11.12 
x3 quaternary, 2; and x2 S60 6 6 0 
binary 
Class 13 (Tetragonal bipyramidal) center of symmetry, d = 0 (11.13 
Class 14 (Ditetragonal 000 0 ds 0 
pyramidal) d=i0 00 ds90 O 11.14) 
v3 quaternary dx dx dx 0 O O 
x, and x2 planes of symmetry 
Class 15 (Ditetragonal bipyramidal) center of symmetry, d = 0 (11.15) 


wien 
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ss 8 (Orthorhombic bipyramidal), center of symmetry, d = 0 (11.08) 


Class 16 (Trinonal dy —dyO dy dis —2dx 
pyramidal) d= — doo doe 0 dis —dy ian 2d); 
X3 trigonal d31 d3 d33 0 0 0 


(11.16) 


Class 17 (Trigonal rhombohedral) center of symmetry, d = 0 (11.17) 


Class 18 (Trigonal) dy —dy 0 dy O 0 ) (Quartz) 


trapezohedral) d={0 0 00 dy —2dy (11.18) 
x3 trigonal, x, binary 0 0 00 0 0 
Class 19 (Trigonal bipyramidal) dy —dy 00 0 —2dx 
xs trigonal, «3 plane of d={-—dxy dx» I™0 0 —2d, (11.19) 
symmetry 0 0 000 O 
Class 20 (Ditrigonal pyramidal) 0 0 0 O diy —2dx\ (tourma- 
x3 trigonal, x2 plane of d =| —dx dw O dis 0 0 line) 
symmetry d3y d3y d33 0 0 0 (11.20) 


Class 21 (Ditrigonal scalenohedral) center of symmetry, d = 0 (11.21) 
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Class 22 (Ditrigonal 


bipyramidal) dy, —dy, 0 0 0 
x3 trigonal, x3 plane of d={0 0 000 
symmetry 0 0 00 0 
%2 plane of symmetry 
Class 23 (Hexagonal 0 0 O dy dis 0 
pyramidal) d={0 0 O ds -—dy O (11.2 
x3 Hexagonal dx dy dx O 0 O 
Class 24 (Hexagonal 0.0 0 dy 0 O 
trapezohedral) Z=if © 0 0 —dy 0 (41.3 
x3 hexagonal, x, binary 0 0 0 0 0 O 
Class 25 (Hexagonal bipyramidal) center of symmetry, d = 0 (11.25 
Class 26 (Dihexagonal 0 O 0 © hes 
pyramidal) d={0 0 O ds 0 O (11.26 
x3 hexagonal, x2 plane dx, dg, d33 O 0 O 
Class 27 (Dihexagonal bipyramidal) center of symmetry, d = 0 (11.27 
Class 28 (Cubic tetrahedral- o 0 0 dy 6 0 
pentagonal-dedoca- d={0 0 O O dy O (11.28 
hedral) 0 0 ©. 9 0 dy, 
X1, %2, x3 binary 
Class 29 (Cubic pentagonal-icositetrahedral) d = 0 (11.29 
Class 30 (Cubic, dyakisdodecahedral) center of symmetry, d = 0 (11.30 
Class 31 (Cubic, hexakis- o 0 0 a& 0 0 
tetrahedral) d=10 0 0 0O dy O (11.31 
X , X2, X3 quaternary . 8. 6 @ 0 dys 
alternating 
Class 32 (Cubic, hexakis-octahedral) center of symmetry, d = 0 (11.32 


Whenever a center of symmetry exists the piezo-electric property vanishes 
since a center of symmetry requires d’ = (—J) dJ = —d = —d’. Also 
d = 0 for the pentagonal icositetrahedral class. 

Classes 6, 11, 12, 24, 28 and 31 polarize only by shear. 

Classes 1, 3, 4, 7, 10, 14, 16, 20, 23, 26 can be polarized by hydrostatic 
pressure. As an example of this let us consider tourmaline (which is ditri- 


gonal pyramidal). For hydrostatic pressure, X; = X. = X3, X, = \X;5 


= X»s = 0, whence from the polarization stress matrices we find, D, = 0, 
D,z = 0, Ds = (2d3: + dos) X pressure. As ds, = 0.75 X 10° and dss = 5.8 
X 10° for tourmaline, we get 7.3 abcoulomlos per cm’ per dyne per cm’. 











seabed 








ishes 


Also 


tatic 
litri- 


oz () 


ad 
ne A 
x 


mM’. 











FEE it ad iS Sa eg ARABS 5G 





MATHEMATICS OF PHYSICAL PROPERTIES OF CRYSTALS 39 
SECTION 12 
THE CONVERSE PIEZO-ELECTRIC EFFECT 
A stress X causes an electric induction 
D= dX : (11.2) 
and a strain 
e= SX (8.7) 


If the charge is allowed to leak away a further strain occurs, at constant 
stress. This is the strain that would be gotten if the stress were originally 
applied with surfaces rendered conducting: 


ewe SA (12.1) 


In the first sort of stress, the work per unit volume done on the crystal by 
establishing the stress X is: 


W = 3X.e = $X.SX ... (8.4) 

The energy stored electrically in the medium is: 
We = 2nDk'D ...... (12.2) 

while the work done on a conducting crystal is: 
W° = 4X.S°X (12.3) 


If a crystal be stressed in its insulated state by expenditure of energy W, 
the charges then absorbed by an external circuit taking up energy Wy, 
the strain changes from ¢ to e° at constant stress so that the stresses perform 
additional work 


W, = X.(e° — e) = X(S° — S)X 
and the crystal is left containing energy W°. Whence 
W°=W-W.t+W., (12.4) 
or: 
1X.S°X = 4X,SX — 2nD.k'D + X(S° — S)X 
so that: 
X.(S° — S)X = 4eD.k"D 
If we substitute D = dX we find 
X(S° — S)X = 4nX dk dX 
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so that: 
os CS) Se 
The change in strain caused by rendering the surfaces conducting is 
e° — e = (S° — S)X = 4nd.k “dX 


If the crystal be now insulated and the stress removed, an induction 
opposite sign will occur and because of the assumed linear dependence oj / 
on X the new induction will be equal to the negative of the previous on 
The induction D = —dX indicates an electric field: 


E = 4nk"D = 4ek"dX ........... ey © 


Also, the strain will alter by an amount —e’’, where, since the action takes 
place with non-conducting surfaces: 


This leaves a strain on the crystal, of amount: 
= —e’ = (S°—S)X. Bapeen thas 

From (12.6), (12.7) and (12.8) it follows that: 
eé=d.E.. PO lon eee (FR 


As the medium is in just the condition that an electric field E would put th: 
unstressed medium, (12.9) is the equation of the converse piezo-electri 
effect. It is to be noted that the set of constants that relates polarizatio: 
and stress is the conjugate of the set that relates electric field and strain 
For convenience in notation the converse effect will be written as 


e= gE ... xa Aan (12.10 
where 


g=d.. acne (12.11 


Rewriting (13) as ase = (az'ga “aE we see that 


where 


g@ =a, g6..... er Ppa ree RFy 














mn. takes 


(12.9 
yut the 


lectri 


iza tior 


12.11 


12.12) 
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SECTION 13 
(He CONVERSE PIEZO-ELECTRIC EFFECT AS A NON-LINEAR FUNCTION 


if the strain of a crystal is not strictly a linear function of the electric field 
Sond power as well as to first power terms. That is, the equation e = gE 
ich gives the strain e in terms of the electric field E through the 18 
constants g) must be modified to include terms E,E;. All such terms are 


E sing it we must relate the components of strain to field terms of the sec- 


included in the symmetric matrix (EF E,). 
° + . . _/ 
\ transformation a that replaces E by aE’ also replaces E. by E,a so 
3 " . fy! . (mm ° ° 
Sthat (E E,) is replaced by (a. E,a), that is (E EF.) being self-conjugate, 


* transforms similarly to the stress matrix. We may rearrange this as a one 








> column matrix similar to the stress matrix YX, as follows: 


Ej 
E} 
Ga) 9 | = 
Fak; 
E\E: 


les] 


(13.1) 








We may now relate the strain to E and E through the two matrices g 
and G: 
e= g9E+GE . (13.2) 


If transformations permitted by the symmetry of the crystal are per- 
formed, g’ must equal g and G’ must equal G, this allows us to simplify the 
matrices; g is no different than before and hence vanishes for all types having 
centers of symmetry (and for the pentagonal icositetrahedral class). 

Rewriting (1) as a,'e = (az ga aE + a-'Ga ‘aE we see that 


eé=ZE'+G’E 
where 
| a. 
§ = a ge 


G’ =a,'Ga....... oun 


The matrix G transforms as the elastic modulii matrix does but G;; + 
G;;. Applying G’ = a.Ga we arrive at the set of matrices that follow 
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x 
i 
Triclinic (36 consts) Monoclinic (20 consts) : 
(Gu Gy Gig Gu Gis Gi Gu Gy G43 0 0 Gis) : 3 
| Gar G3 G33 Ga G35 G6 | (1 3 4) G31 G32 G33 0 0 Gag is al : ri 
1Ga G42 Ga Gus G45 Cag | <a 0 0 0 Gu G45 0 | 
Gs Gs Gs3 Gs G55 Gs | 0 0 0 Gs Gs 0 | ’ ( 
(Ger Geo Ges Ger Gos Ges) Ga Ge Gs O Ge) ' 
Orthorhombic (12 consts) Tetragonal Classes 9, 10, 13) \' 
(10 consts) 
Gu Gry Gig 0 0) 0] Gu Gy Gis 0 0 G6 
Go Gr Go3 0 0) 0 Gi Gu Gi3 0 QO —Gi su 
731 G3 G3 0 0 O | G3, G 0 O 0 g 
iC 31 G32 G33 ) | (13.6) Ga 31 G33 é | 437 ge 
| 0 0 0 Gas 0 0 0 0 0 Gu G45 0 la 
6 6 8 6 tad o © ¢+Gikh © 6. 
0 8 6. 2 Gee Ge — Gey 0 0 0 Gee h: 
Tetragonal Classes 11, Trigonal (Classes 16, 17) R 
12, 14, 15) (7 consts) (10 consts) ¥ 
( . . W 
GuGuGs0 0 0 } Gu GuGs Gu-Grs 0 
in GuGunO0 0 0 | Gu GuGa 0 0 | 
Ga Gai Gas eo. 4 | (13.8) | Ce Gutin 0 0 0 | (13.9 
0 0 0 Gud 0 | | Ga-Gu0 Gu Gau2Ge 7 
0 0) 0 0 Gu oe 4 — - 0 —G45 G42 2G | 
00 0 0 0 Ge) 0 2G 2Gu2(Gu—Gz) 
Trigonal (Classes 18, 20, 21) V 
(8 constants) i 
Gu Gi G3 Gis 0 0 
Ge GuGs—Gu 0 0 ) 
G 131 G 0 0 0 
| oa Ga G33 ‘ (13.10 
» |Ga —Gy 0 Gy O 0 
10 0 0 O Gu 2Ga 
l0 0 0 O 2G 2Gun —Gy) 
Trigonal (Classes 19, 22) 
(6 constants) 
Also Hexagonal (Classes 23, 24, 25, 26, 27) 
Gu Ge G30 0 0 
Gy Gu Gi3 0 0 0 
G3 Gu G30 0 0 
age (13.11 


00 0 Gud 0 
00 00 Gud 
(0 0 0 0 0 2Gu— Gy) 








/ 





10 








; 
i 
' 


i 
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Cubic (Classes 28, 29, 30, Isotropic Bodies 
31 and 32) (5 Constants) (2 Constants) 
Gu Gey GeO O O Gu Gy Ge 0 0 0 ) 
Gyo GuGre0 0 0 Gy Gu Ge 0 0 0 | 
. . ’ | lp . + } 
Gy Gi Gu 0 0 0 (13.12) | Guz Gy» Gy ° 0 0 | (13.13) 
0 0 0 Gad O Ce 86 8 60 0 | 
0 0 0 0 Gy 0 0 0 0 0G0) 
0 0 0 0 0 Ga 0 0 0 00 G| 


G= 2(Gu —Gy) 


According to this analysis, all bodies suffer a change in dimensions when 
subjected to an electric field. These strains resulting from a field are 
generally much smaller than those strains e = gE present only in crystals 
lacking a center of symmetry. For example, quartz has a strain of about 
6.5 X 10°* cms/cm/ab volt. Glass in a field of 1000 practical volts per cm 
has a strain of about 4 X 10°”, in a 100,000 volt field it has 4 X 10° 
Rubber in the 1000 volt field strains by about 7 X 10° and in the 100,000 
volt field by about 7 X 10™*. The 1st order quartz strain in these fields 
would be about 2.2 X 107’ and 2.2  10°° respectively. 


THE SECOND ORDER PIEZO-ELECTRIC EFFECT 


If the induction stress relation is not strictly linear one can assume the 
induction to depend also on second order terms of the stress: 


D = dX + p (XX,) 


where (X.X,) is a single column matrix formed from the 21 elements of 
XX, and p is a matrix of the 63 elements pu: . . . p33,3. 

Since X transforms as X’ = aX’, (XX.) transforms as X’X. = aXX.a. 
In the same way that a was formed from a we can form a matrix & that 
transforms the single column matrix (XX,) through (XX,)' = a(XX,)’. 

aD = ada X + ap(&) &(XX,) or 
D! = d'X' + p'(XX,)’ 
where 


d'= ada and p' = ap(a)” 


The first order effect is the same as before. With the relation p’ = ap 
(a) * we could perform the operations of symmetry permitted by the 32 
crystal classes and obtain the reduced matrices. However since & has 484 
elements we shall limit ourselves to crystals with centers of symmetry. 

As X is unchanged by an inversion through the origin, a is the idemfactor 
for this transformation and a is —J, also (a) = J. Therefore D’ = —D 
= D so that D vanishes. 
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It is to be noted that although there is a sort of reciprocity betwe: 
first order piezo effect and the converse effect, in that the matrices {, 
are the conjugates of the matrices of the other, there is no such recip 
in the second order effects: if a center of symmetry exists no polariz 
can be brought about by stress either as a first order effect or as a si 
order effect; if a center of symmetry exists an electric field can caus: 
strain through the second order effect but not through the first order e| 


Dielectric Constants at Constant Stress and at Constant Strain 


Let us consider a unit crystal cube, initially unstressed, unstrained anc j; 
zero electric field. We write k” and k’ for the dielectric constant matrices 
at constant stress and constant strain, respectively, C” as the elastic con- 
stant matrix at constant electric field E, C° as the same for zero field. We 
study a cycle consisting of a strain caused by application of an electric field 
E at zero stress followed by a stress applied at constant E to reduce th 
strain to zero and completed by conducting away the electric charges at 
zero strain so that the body is left in its original state. The cycle is de- 
scribed by the table: 


, Change | Change Change i Chang " ' 

Operation in ease in Strain Diplecsinent Cuneat in Field Energy Put In 
Apply E..... 0 dE A pg | gt. | E.kPE 

4a 8r aver 
Apply -«....| —C®d.E | -d.E| — i (ke — RE 05 E.d C¥d.E. - 
1 | / 
Apply -E...| CH | 0 - x RE oi x a | EE 
vT rT 
’ BE 
whence C = C’ ; also 
eo FS re (13.14) 


SECTION 14 
PyRro-ELECTRICITY 


If the electric polarization brought about by heating some kinds of crystals 
is simply a function of the uniform temperature change, that is if this 
polarization can be produced by taking the whole body quickly from the 
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rm temperature f) to the uniform temperature é + ¢ the pyro-electric 
effect could be described by the equation: 


dD, pi 
D.) = ¢{ pe). (14.1) 
Ds 


where p is the pyro-electric matrix. 
This can be approached in another way by considering the polarization 





as due to the uniform strain. We may hence write, since Y = Ce: (i.e. 
stress matrix = elastic constant matrix times the strain matrix) 
P = dX = dCe 
where e is the strain brought about by the temperature change ¢. If 
A, 0 O 
4={0 Az O ] is the temperature expansion matrix we have: 
0 QO As 
A;) 
Ay 
. |A 
P=td | | 
|0 | 
0 | 


lo ) 
Now since d has 3 rows and the A matrix has but one column the product 
dCA has 3 rows and one column so that we may define p as dCA. 
As D of D = ftp transforms by D’ = aD, so does p: 


U 


p’ = ap 
When a center of symmetry exists a permitted transformation isa = —/, 
whence p = —p’ = —p so that p= 0. No pyro-electric effect (on this 


theory) could exist for a crystal with a center of symmetry. 
If a binary axis exists and is chosen as x3 we have 


pi —-1 OO pi “=o 
po} = 0 -—10 pe} =| —pe 
ps 0 O1 ps p3 


whence for this case 
p=10 ; a (14.2) 


If another binary axis exists at right angles to this one we find p = 0. 
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This is seen seriously to limit the number of classes showing this kind ¢j 
pyro-electric effect. In fact we find p = 0 for classes 2, 5, 6, 8, 9, 11, 1) 
13, 15, 17, 18, 19, 21, 22, 24, 25, 27, 28, 29, 30, 31 and 32. The expressio 


0 
p = | 0 ] describes the pyro-electric effect in the classes, 3, 7, 10, 14, 16, 2 
3, 
P1 
23 and 26, while the expression p = | pz} describes class 4, and only clas 
0 
1 is described by 
Pl 
p=|pol.. Pee Gs, 


It is to be noted then that this theory excludes many classes ordinari 
described as pyro-electric, such crystals as quartz in Class 18 for example 
Consequently it would seem that whether or not this effect exists we must 
seek elsewhere for the explanation of the effect in quartz. 

The effect can easily be explained as due to non-uniform temperature 
which causes stress which in turn give rise to electric phenomena in piez 
active crystals. For example a suddenly chilled crystal has its outer layers 
in a state of tension. This would produce just the pattern of positive and 
negative charges that one actually observes. As to whether the first effect 
exists, much argument between Lord Kelvin and others seems to have left 
the question still uncertain. 

In pyro-electric crystals we would expect to find a difference in the piez 
constants measured isothermally or adiabatically. If a temperature 
change ¢ causes an electric displacement D = ptf the application of an electri: 
field E should cause a temperature change ¢ given by a relation such as: 


Joos  EIPTOTTOT ET CRT eee (14.4 


Also the temperature coefficient of expansion, A, (for a crystal with faces 
rendered conducting) would differ from the coefficient A; (for a crystal wit! 
an insulated surface). 

If a crystal at temperature ¢o has suddenly applied to it a field E the 
temperature rises to fo + gE and the crystal strains, because of the converse 
piezo-effect, by amount ¢« = ga where gq is the adiabatic converse piezo 
matrix. If the field is now removed isothermally a further strain ¢\f 
takes place. If the faces are short-circuited and the temperature restored 
to fo a further strain A,t = A,gE takes place and the crystal is then in its 
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tial state. Equating the sum of the strains to zero we find (ga — gi) 

Ek = AwE or 
3 = Aw (14.5) 
Let the initial state of a crystal be, temperature = /o, stress, strain and 


feld = 0. If the (electrically insulated) crystal is heated by amount /, 
. strain A f is caused and also an electric displacement D = pt. There now 
exists an electric field E = rk ‘pt. Let this field be discharged at con- 


| stant temperature, giving a further strain of g:E = 4agik pt. The crystal 
f 


is now short-circuited and if the initial temperature is restored a strain 
A,t follows. The crystal is now in its initial state. If we equate the 
sum of the strains to zero we find: 


A, —A = 4ngik''p (14.6) 
SECTION 15 
THE THERMO-ELECTRIC EFFECT IN CRYSTALS 


It should be possible for an electric field to be set up by a temperature 
gradient. Let us assume that the vector 7 is the temperature gradient and 
is related to the vector field & through the matrix II by means of the equa- 
tion: 

Ih The Ths 
E = TIT where II = [{ II Ie Ik (15.1) 
IIs; Ize Ws, 


Examination shows that II transforms through 
Il’ = alla. (15.2) 


For Class 1 the II matrix has the 9 terms of (15.1). Class 2 has a center of 
symmetry. For a center of symmetry a = —/ but a = —/ causes no 
change in (15.2) so that class 2 has 9 constants. The thermo electro effect 
is not killed by the presence of a center of symmetry. The ordinary thermo- — 
electric effect of metals is a case in point. 


—1| 0 0 
If x3 is a binary axis a = 0 —1 0] and II reduces to 
9 OF 


In The 0 2 
II = | Ia Ie O ? (15.3) 
0 O  Isg 


Examination shows this form to answer for classes 3 and 4 and 5. 
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If x; and x3 are binary (15.3) reduces to 


re Y- 
1={0 IT» 0 ao $4 
O O IIx 
which described classes 6, 7 and 8. 
For x3 quarternary alternating: 
IIy; Ty O 
Il = | —IIy Ty, O ; ‘ a aaaA clin, os 15.5 
0 O- Ix 


This is found to handle classes 9 and 10. If x3 is quarternary alternating 
and x is binary: 


Ty, 0 O 
Ii = 0 Thy, 0 eA ‘ oad (15.6 
O O IIgs 


which is found to cover classes 11, 12, 13, 14 and 15. For classes 16, 17 
19, 23 and 25 II reduces to the form (15.5). 

If x3 is trigonal and x; is binary the matrix is (15.5) which then handles 
cases 18, 20, 21, 22, 24, 26 and 27. 


For cubic crystals, not only are xx. and x3 binary as for matrix (15.4 


1 010 
but the vector [1] is an As, for which a = {| 0 O 1 ] whence, for classes 
1 100 
28, 29, 30 and 31 we find the matrix: 
Ty, 0 O 
I={0 IIy 0 Bis ere ay (15.6 
a ee 


Reports of a pyro electric effect in quartz should probably be attributed 
to nonuniform heating exciting the piezo electric effect. Reports of a pyro 
electric effect in such crystals as topaz and colemanite which have a center 
of symmetry and hence cannot be piezo electric should probably be at 
tributed to this thermo electric effect. 


SECTION 16 
THE PROPAGATION OF LIGHT IN CRYSTALLINE MEDIA 


Maxwell’s equations are: 
CV X B= 4nj 
CV XE=-B 
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n C is the velocity of light in free space, E is the (vector) electric field, 

: the induction current and B is the magnetic induction. In a crystalline 

lium the current is given by 4x7 = RE where k is the dielectric constant 
matrix), whence: 


CY XB 


Il 


kE (16.1) 


CVXE=-B (16.2) 
As the divergence of the Curl is always zero: 


VRE = V.j = 0 


a (16.3) 


re, ee P 
applying ry to (16.1) and substituting (16.2) in the result: 


—-CVXVXE=RE or 
C(V.V — VV) E=kE (16.4) 
We shall try as a solution: 
E=- £,"™.... .... (16.5) 


where Ey is the vector amplitude of the electric field, i is ~/—1, r is the 

radius vector from the origin to any point, g is a constant, » is the unit 

normal (at r) of surfaces of equal phase, and w is 27 times the frequency of E. 
Substituting (16.5) in (16.4) we find: 


9 


w 


E — nE.n = —. 
gc 


kE (16.6) 


“ . . ¢ - W. ° = 
Examination of (16.5) shows that — is the phase velocity along n. Writ- 
q 2 


9 


F —j. +2 Ww 
ing 4rk 7 for E and V’ for — we have: 
@ 


; . alt 
k"j —njek'n = =j (16.7) 


c 


This equation is independent of the absolute value of j so let us restrict j 
to being a unit vector. 


VRE = 0 = VekEe"*" = j.nige" 
whence jn = 0 oe on 


That is, the current is always normal to the direction of propagation. 
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Multiplying (16.7) thru by the prefactor 7. and cancelling the ter 
jen we have: 
r2 
oe ee: 
jck oli Soe (1 


This tells us that the velocity is a single valued function of the direct 
of the current. 


With the idemfactor /, (16.7) may be written: 
k a j= n(k'j 
— GI) j = men 


: , : . we. %@ 
If we multiply this thru by (: — 5 1) we get 


Ce 


: - Vy? —|l - 
j=l(k - al Ms INe.... sono 496.10 


Multiplying this thru by , and dropping the scalar factor (k ‘i)n: 


72 —1 
n (et - £1) n=0Q.. er rT Ee) | ty 


If the axes are so chosen that & is a diagonal matrix (16.9) and (16.11 
become: 


l 2 ii 2 i 
2 3 / 
ee ae ee = ee ee 16.12) 
c? ky Roe k33 _ 
2 2 3 
ny No nN3 
=(Q.. 16.13 
1 v2 + i “hs 1 7% (16 


ky CC ho» Cc? k33 Cc? 


Examination of (16.13) shows that (16.11) must have two values of 
V’ for each value of the vector normal n. As V’ is a single valued function 
of 7 there must be two distinct values of 7 (j’ and j”’ say) for any particular 
n; and given n, only waves having their current vectors in the directions 
of j’ and j” can be propagated. A ray in the direction N but not having 
its j in one of the directions j’ or j’’ will be broken up into two components 
having their current vectors along j’ and j” respectively. 

If the velocity V; corresponds to 7’ and V; to 7” we have by means of 
(16.10) since n’ = n”: 


a 1 Vi -1 V3 ty ly “1 spy 
Jef =M(k — al k ee n(k 7)-n(k 7"")en 


(The quantities in the braces are scalar) 





y 





10 








J , iy Sc 4 
j = D where D is 
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By means of the identity 





(w' —v') = —u'(u — vv" 
since 
Br (et -B)- Bo", 
c c c 


the idemfactor can be multiplied into an adjacent matrix giving 


72 \-! 72 \-1 
(eB) (8's 
c c* 
2 ™ v2 1 , v2 a 
= @ Ne k eae I —_ k — ns 
V2 — Vi ( ce ) ( c ') . 


0-0=0 


tor 


so that j’ and 7” are mutually perpendicular. 
SECTION 17 
Tue E ectro-Optic EFrect 

The velocity of light in a crystalline medium is a single valued function 
of the unit current vector 7 
ae 
a =Sek 7 re (16.11) 
where c is the velocity of light in vacuo and & is the dielectric matrix, also 
1D 
dt © 

We developed the induction as a linear function of the electric field, 
deriving the relation: 


4rD = kE ere (6.1) 


If the induction is not a linear function of the electric field we can improve 
on eq. (6.1) by adding second order terms: 
4nD; = kaE: + ki2Es + kigEs + knEi + kak + kak; 
a } hwE.Es a 1 hiwE3E> + << } hk E> 


or 


4rD; 


Il 


(Ris + hak, + } hiek2 + hisE3) Ey 
+ (hie + 3 higEs + hinE, + + hE) E2 
+ (Ris + 4 hisEy a } hisE> + hi3E3) E; : (17.1) 
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Examination of (17.1) suggests that we might consider the k’s as | 
lineraly modified by the field. Writing k as a single column matrix: 














ky Ky 
Rep Re 
Rss 

K = - 
Rey 
ki: ° 
kin) (Ke 

we may write 
K=K°+hE.. ee ae (17.2 


where K° is the dielectric matrix for vanishingly small fields. 
We can develop the modified reciprocal matrix in the same manner: 














Ky kr 

= koe 

kit kia is | on at 

k= | Rig kee kes | Kt=/" } = 14 

lkis kos kas | we “ 

. ki3 

Ky" kis 

where 

K'=K "42k ares | B 


It is to be noted that K‘ is not the reciprocal of K but merely a symbol 
for the single column matrix formed from k* in the usual way. Taking 
reciprocals of both sides of k’ = aka. we find (k')’ = ak ‘a.. That is, 
k* transforms exactly as did k. Whence, K ‘ transforms exactly as A 
does, i. e. 


wg? 


K* =aK" 
We can rewrite (17.3) as 


aK'=aK" + (aza,) aE 


or 
K" = (K")'+7E' 
3’ = a 20... PET re eg ot (17.4 
In case a = ac the z’s transform as do the conjugates of the piezo 


electric constants, d. Of the transformations permitted by the symmetry) 
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e 32 crystal classes only those of the trigonal and hexagonal systems 


1 
to have @ 


These 12 classes must be examined individually but 


other classes may have their z matrices copied from the corresponding 


matrices. 
\pplying 2’ 


aza,. for a rotation of 120° about x3 we find for class 16 


211 — 220 233] 

— 21 222 21s | 

0 0 2g3| 

has 241 31 O 
Zs1 —2a O 
—l2 —2n 9 








The remaining 11 classes may be derived from class 16 by operations for 
which either a ' = a, or a center of symmetry exists. Consequently, we 
may form our z matrices from the d,’s in all cased if we leave out the 2’s. 

lhe electro-optic effect can be put in terms of the polarization instead of 
the field by substituting in (17.3). 


E = 4 (k—1)' P_ whence 


K = K+ »P (17.5) 
where 
n = 4nz(k — I) 
Conversely .. (17.6) 


1 
z= —n(k —T) 
a" 


The » matrices transform exactly as did the z’s and hence may be formed 
from the d.’s but omitting the 2’s. 


SECTION 18 
THe Prezo-OpTicaL EFFECT 
If the dielectric constants of a crystal are changed by the application of 


stress, this may be represented by: 


K~™ = (K")° + #X (18.1) 


where the 36 constants my °-: 
We may then form k* as 


kit + miX,, 
ko = kin + 16iX i, 
kis + WsiX i, 


mes are stress-optical constants. 


Ro + W5iX i 
Ne + waiX: (18.2) 
ks + 13iX ; 


kin + Weir i, 
kee + 92iX;, 
kag + maiXi, 
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As the velocity of a light ray of unit current vector j is given by 
y2 


a ae 
jek J = eo 


We can, by (18.2) and (16.9), compute the change in the velocity caused ; 
by the stress, if we know the constants 7. 
Altering (18.1) to aK’ = aK + ama ‘aX we see that: Te 
KY = K+ 'X’ where r= ama’... 18.3 


[he alteration of K ~ can be expressed as a function of the strain by 


substituting ce for XY in (18.1). 








1 -—1° p—]° 
K' =K + ae = K + me me (18.4 
. Tr 
m= TC, r= ms . pian (18.5 
Operating in (18.4) as we did on (18.1) we find m transforms as 
m' = ama, (18.6 
Applying the crystal symmetry operation to these matrices shows that 
they reduce to the following T 
Triclinic system (mri Ty. ™3 Ti ™M5 m6) 
36 constants 21 122 123 To, Te5 Toe | ren 
The m matrix | 
| 31 732 133 34 135 36 | ; , | 
is entirely 
| 41 Wao W43 144 Way 146 | 
analogous 
| 751 W52 753 Ws WHS Woe | 
| wet We2 Weg Wea Wes re) eee Ee ee Lecck eed 
Monoclinic system T1172 7m30 O mig I 
x3 is binary |aro1 Woe 23 O O a6 : 
‘ The m matrix 
20 constants (131 Te 133 0 O m6 | . . 
is entirely 
jo 0 0 ta.%08 analogous 
10 -O O weg wes 0 — 
beta ie ne) AE BE OS wicepecus we ROS (18.8 
Orthorhombic (mi to 730 0 0 ( 
system x3 | mra1 ™ 730 O O The m matrix 
is binary 731 732 1733 o. @ ® is entirely 
12 constants 0 0 0 m0 O analogous 
0 0 0 O ws 0 (Rochelle | salt 
°-6.0 60 BO wal ..... ee eee 











trix 


alt 
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fetragonal system 
rs isa four-fold 


axis 
(Classes 9, 10 
& 13) 
9 Constants 


(mu Te 3 0 O 


| Miz 1 ™3 0 O 
| 731 TW) 733 O O 
0 0 O my 0 
10 0 0 0 mm 


lr —7 0 0 O 


Tetragonal system 











x3 is a four-fold axis 
x, isa binary axis 
(Classes 11, 12,14 & ( 





x3 is a trigonal 
. | 712 
axis 


% , . | 31 
x, isa binary axis | 


(Classes 18, 20 & Py si 
21) lo 
8 Constants | 


Hexagonal system 
3 is a Sixfold axis 
x, is a binary axis 
(Classes 19, 22, 23, |0 


(mu 
| 


| 112 
| 
731 


24, 25,26 & 27) [0 
6 Constants 0 
3 Constants 


Cubic system 
| 
| 


) 0 O 


7713 —T40 O 
731 733 0 0 O 
14; 0 T1440 O 
0 0 0 
0 0 0 


712 713 Gg Go Gg 


— 116 


(aru tm 730 0 O 
Te Ty M3 0 VO O 
‘331 Ts 7330 0 0 
Ti 0 O 


— 12) | 


15) (0 0 0 0 mu 0 
7 Constants 10 0 0 0 O me 
Trigonal system Wi 712 113° 14 — Ts O 
vz isa trigonal M2 71. 13 —%14 = Tas O 
axis 31 731 733 «(OO 0 0 
(Classes 16 & 14) 14, O 144 «45 250 
17) —%52 7520 —m45 a4 2H 
11 Constants 0 0 0 T2514 (ty 
Trigonal system 
Tm 8 T2™m3 «THO O 
| 


744 241 


T7i4 (411-712) | 


m6) 

0 The m matrix 

0 | is entirely 
analogous 

0 : 


66) (18.10) 


\ 
} 
Che m matrix is entirely 


analogous 





(18.11) 


) The m matrix is analo- 
gous except that 


Ms = M52 
Msg = M4 


m1, — My 
? 


| Me = 


} (18.12) 


The m matrix is analo- 


gous except that 


Ms = M4 


} mM, — Mie 
Me = 

| 2 

| (quartz) 

} 8.13) 

} (18.13 


} 
| 


The m matrix is ana- 


Wit 713 0 0 O 
37330 O O | logous except that 
QO O m0 O My, — My 
00 0 ry 0 a 2 
0 0 0 O (an - mu) (18.14) 
Tu Te Te O OO O ) 
ltm2 Tu 7T2 O O O | 
T2 T2710 0 O The m matrix is en- 
0 0 0 ry 0 0 tirely analogous 
0 0 O O waO | 
000 0 0 se) (18.15) 
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For isotropic bodices, the + matrix is formed by setting mu = (9,— 7, 
in the r matrix of the cubic system; the m matrix is similarly forme 


i Mm, — My 
putting m4 = ——, , 
Isotropic (mu To ™2 O 0 0 The m matrix 
bodies [m2 71M. «O 0 0 

| ar12 Tem O 0 0 | myetc = = 2 
10 0 0 m—m 0 0 | . 
10 00 O m-m OO | 
(0 0 0 0 0 m-—m,).... (18.16 


SECTION 19 
APPLICATION OF THE ELECTRO AND Piezo OpTIcAL EFFECT 


In the equations K’ = K-" + zE and K' = K+ mE, etc. the Ks 
are to be used in forming k ’ for the equation giving the velocity of the light 


> 


used namely Ga j-k'j. Obviously then K~* should be formed from the 


squares of the reciprocals of the refractive indices, the lower three members 
being zero. After applying the electric field or strain a transformation oi 
coordinates may be necessary to rediagonalize, i.e. make = = K;" 
Kz. = 0. From the rediagonalized K-' we may write the new principle 
refractive indices by taking the reciprocals of the square roots of Ay. 
Kz" and K;". It should be noted that if x” = K;' + A; then 


, x 
wit Ow = wi — 5 TT UR Le 


For a given direction of the wave normal there are two velocities, a wave 
splitting inte two components traveling with different velocities. By defini 
tion the refractive indices, ya and yw» in a given direction are the norma! 
velocities in that direction divided into the velocity of light in free space 


; : l lu 
Whence in a path length / there are = waves in one component and “ 


waves in the other, where A is the wave-length in vacuum. Consequently ii 
1 is the thickness of the crystal along that path the two components can 
recombine after passing through the medium but they are out of phase by 


l : 
Tr. (ua — ue) whole waves so that the light which entered as plane pola 


ized will leave elliptically or circularly polarized, except when 2y is an intege? 








| 





3 
4 

















MATHEMATICS OF PHYSICAL PROPERTIES OF CRYSTALS 57 


(he quantity B = pa — me is known as the birefringence. 


1B 
= — | (19.2 
x ) 
If a phase difference of y'5 wave can be just detected, using a wave length of 
6000 A and a path length 7 = 1 cm the just detectable birefringence would 
be B = yA = 2X 10°: if the path were 10 cms the detectable B would be 


l 
0.2 X 10°. Obviously this detectable difference between refractive indices 
is much smaller than could be detected by measuring each refractiv¢ index 
and subtracting. ‘ 

It is customary to choose the coordinate system so that looking along x2 
the very lowest refractive index is for polarization in the plane of x; and the 
very highest for polarization in the plane of x3. That is, the x2 axis is the 
axis along which light should be passed to get the greatest birefringence. 


Birefringence in any Direction 


If the axes are so chosen that K is diagonal and K3; > Ke > A, then, 
somewhere in the plane perpendicular to x2 are two directions, the optic 
axes, along which there is a single normal velocity. These directions make* 
equal angles V with the x3 axis where 


a f= — Ky’ 
sin =-+ a Ky 


or P , eee: (19.3) 
sf =, Kx 
cos V = + kK; — Ky 


Also the two refractive indices ua and yw» for a wave normal making angles g, 
and ge with these optic axes satisfy the equation: 


2 
“3 = (Ki) + K3') + (Ki’ — Kz") cos (gi — ge) 
Ma 


jbo 


orl 


= (Ki + K3') + (Ki — Kz’) cos (gi + g:) 


= 


whence 
) ( 5) - 2 (us — Ma) (Ma + bp) 
; ba Mo 
= (Ky' — K;') cos (g1 — ge) — cos (g1 + ge) 


* Theory of Optics, P. Drude, pg. 320. 
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as B = ws» — ua is the birefringence we have: 


Ma Mb >—1 7>—1 ° . 
B= : (K, — Kj ) sin gi sin go 
Ma + Me 
By spherical trigonometry: 
cos g, = cos V cos + sin V sin @ cos @ 
cos go = cos V cos 6 — sin V sin 6 cos @ 
where @ is the angle the wave normal makes with x3 and @¢ is the angie th; 
plane containing the normal and x3 makes with 2}. 
From (19.5) it follows that: 
sin g; Sin go = 


Yi — cos? V cos* @ — sin? V sin? 6 cos? ¢)? 
° 2, Sep » 3 2 2 
—4sin°' V cos V sin @ cos @cos 


Hence if the rediagonalized K 


+A 
+ 
+ ds 
0 
0 
0 








then 


ee 

Be Ma Bb 
Ma — Md 
(Ki' — Kz’ +4: — 43) — (Kz' — Ky" + de — As) cos’ 0 
— (Ky — Kz" + A, — ds) sin’ 6 cos’ ¢ — 4(K7" — Kr" 
+ A, — de)(K3' — Kz" + de — Aj) sin’ 6 cos’ 6 cos’ 
3 
B= B+ 5 


For most practical purposes we may take 


(A; As) sin gi Sin go........ 


3.3 
‘ Ma bb 
Ma + Me 


3 
=f 
2 


where u is some intermediate value of the refractive index.* 


2 
* Note: It might seem that as Kj! = Kj + A; gives us wy = wo — 3 Ait guia; 
we could form the 3 principal birefringences directly from the y’s instead of using (6a - - « ¢ 
From the u expressions we would get Bs = ws — wi — 4(u3d2 — ids) + § which differs 
from 6 if uw. # w,. Equation 66 is correct; the one from the u expression is an approxi 
mation. : 
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In a few special cases (18.7) may be simplified. If (¢, ) falls along an 
ptic axis g. = V — V~ and gs = V + V” whence 


Opt 
i 


B, = ~[(d: — ds) cos V° — (de — Ag) sin? V°] (19.61) 


B = (us — w1)(1,— cos’ ¢ sin” V°) 


3 
+ ¥ (Ai — As — (A: — As) cos” ¢). . (19.63) 


if = 0 
3 
(us — wx)(sin” V° — sin’ 6) + aa (Ai — dg — (A; — As) sin @). . (19.64) 
if ¢ = 90° 
B = (us — w1)(1 — cos’ @ cos’ V°) 


3 
ob ‘ (A; me As _ (As — As) cos 6) 


The Electro Optics of Quartz 


° P — 
For quartz, in the equation K 


is in practical volts* 
0 | Zu = 0.47 X 10°” 
Za1 Z = 0.20 « 10°” 
0 —za 0 
0 —zy 0 








Obviously the E; component produces no effect so we shall examine the 
1 
effects due to the components £; and EF separately. If E=|0] 2, 
0 
Ky" + 2*nFy 
x;* = tuk, 
K;" 
tu Ey 
0 
( 0 


* Computed from F. Pockels data, see his Lehrbuch der Kristall-Optik, (B. G. Tuebner). 


which can be diagonalized 
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. —24 E, iis 
by a small rotation about x, of amount @ = Ko Kz giving 
1 — 3 


(Ki + mn FE; rm ” 10 
| p—1° _ — — X 0.47 X 10“! 
| Ky? aaa 24) ” 2 

| i , 
K L and B=mt " xX 0.47 * 10 a: 

| 2 

0 Ms 

0 y 





The greatest ‘“‘added birefringence” is gotten by viewing along x3, whe 
AB = 1.544° X 0.47 X10" E,. If E; = 10°AB = 1.73 X 10° a quantity 
detectable if the path length is about 1 cm. Viewing along 2; (the opt 

axis) is complicated by the rotation of the plane of polarization in quartz 
Homogeneous strains have never been found to alter this rotation, but th: 
rotation complicates and partly masks the birefringence phenomena. |; 








0 
E=1|1 1] we find 
0) 
pan} ® 
1 
K;" 
- Ky" 
vdcate 0 
—2 Fe 
| —21 Ee 


Rotating the coordinate axes through 45° about «x3; then applying the 
transformation 


: © a 
a= 01 —4 
— © 1 
where 
ae 
sia 2 “41 £2 
ile -¥ Ke 
; —Bes 
we find: 
Ky" + su Ea) aBo ne 
Ky" — tnE ™ 2 
wv K;"° Mi 
K = 0 B= mt z 211 Ey 
0 3 
0 
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h is identical to the yu for field along x; , but the final axes in this case 
do not coincide with the final axes for E = E, , but again, the greatest added 
birefringence is utilized by viewing along x3. In the second case the Nicols 
vould be best set along x; and 2», i.e., at 45° to x1 and x; whereas in the 
first case they would be best set at 45° to x, and x». 


The strain Optics of Quarts 


K"'= K" + me 


(mu My. M43 mys O 0 ) my = 138 
| my My, M3 —M14 0 0 My = .250 
lms, ma ms «0 0 0 | m3 =  .259 
m = | mas my, O mas O 0 | where* my= _ .029 
\0 0 O QO mas may | m3 = .258 
My, — My | m33 = O98 
lo 0 0 0 ms — | 
{ 2 ) my = —.042 
mya = — 0685 
If the strain is a simple tension along 2 , 
~—1° -—1° 
(Kj + mx €1| (Kj + mye) 
punh® = | naa? 
es + mye which diagonalizes, Ki + mye! 
, r—1° — ep 
K' = |}Ks + mae! through a small KY = ‘ee + ms ex! 
Mar €1 ' transformation to: 0 


0 | 0 


( 0 ) 0 


applying 18.63 or 18.64 we find the birefringence along x, to be: 
yf 
Bu = ws — fi + > (my - m31) €1 = .0091 — .0148 €1 g 


Similarly the birefringence along x2 is By = .0091 — .225 4& and By = 
0 — .207 &. With a strain of 10 * which is about a tenth of the breaking 
strain, By would be 20.7 X 10°, a quantity detectable in a thickness of one 
millimeter. 

The values of By, --- Biss corresponding to birefringence along x¢ for a 
tension along x2 etc., can be computed in just the same way. But By --: 
Beg require rotations of 45° about x3 to diagonalize, so the birefringences 
can be computed by setting @ = 45° in equations (19.6). 


* Lehrbuch der Kristalloptik—F. Pockels. 
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The following table summarizes these simple strain birefringence efi 
the rows indicating the strain and the columns the direction of light pass 


Bo — .0148 &1 Bo — .222 1 .207 1 
Bo + .222 & Bo — .0148 & .207 é2 
Bo — .298 é3 Bo — .298 €3 0 
Bo + .0536¢: Bo — .0536¢  —.107 e 
(By + 0 Bo + 0 —.107 és) 
(Bo + 0 Bo + 0 .208 €6) 


(the parentheses indicate the 45° transformations). 
A similar table for the electro-optic effect in quartz is 


By — 87 X 10 °F, By + 87 X 10°F, 1.74 X 10°F, 
Bo + 0 By + 0 1.74 X 10 “E, 
Bo + 0 Bo + 0 0) 


Since a driving voltage of E,; = 100 volts may, due to the building up of 
oscillations, cause a periodic strain of e = 10 ‘ina quartz plate, it would 
seem from the foregoing that 99.99°% of any birefringence change must b¢ 
due to the mechanical effect. 


The 18° Cut Crystal 


A crystal, the thickness of which is along the electric axis, x; , the width 
making an angle 6’ = 18° with the optic axis, x3 , can be caused to oscillate 
with a siraple motion alyng its length. (If 6’ is not about 18° or 72° the 
oscillation is not a simple extension along the length, as is shown by the node 
which then lies diagonally across the crystal.) Ona set of axes defined by 
the edges of the crystal block, «1 being in the direction of the thickness or 
x1, x3 in the direction of the width and makes an angle @’ with x3, x3 is the 
length and makes an angle 6’ with x2; on these block axes the strain is ¢: 
Rotating the axes about x; through an angle 6’ we find the strain expressed 
on the crystal axes to be: 


0 

cos’ 6’ 
ezae= sin’ 6’ e2 

. —sin 6’ cos 

0 


| 0 




















node 


} 


Ga DY 


Ss OF 
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whence 


K'=K +e gives us 
7 + (my cos’ 6’ + my sin” 6’ — my sin 6’ cos 6’ )e: ) 
| Ky) + (my cos” 6’ + my; sin? 6’ + my sin 6’ cos 6’)e: 
_ | Kz" + (ms, cos” 6’ + my sin’ 6’)e: 
it — (my cos’ 6’ + my sin 6’ cos ¢’ Jes 


| 


A small transformation removes the 4th term without altering the others to 


Ko 


the first power of small quantities. 
ry . . . e , 
lo obtain the birefringence along the width x3, we set 6 = @ 


(18.65): 


/ 


in equation 


3 
” ae Ms) 2 4 ie 
By = (us — wa) sin’ 6’ + 5 M2 COS 6’ — my cos 6’ + (m3 — m3) sin’ 6’ 


a1 + cos aii ‘ ; 
— mysin 6’ : — m3, sin’ 6’ cos 6’ > e2 
which, for 6’ = 18° is Bz = .00087 + .20 e2 


. . . , 
For the birefringence along the length x2 we set 


6 = 90° + 6’ in (6e) giving: 
uw 
B = (us — wi) cos 6’ + {my cos 0’ — ms, cos* 6’ + (mig — my, — ma;) 


sin’ 6’ cos’ 0’ — my sin 6’ cos 6’(1 + sin’ 6’) e: 
which, for 6’ = 18° is Bz = .00824 + .049 e:. 
SECTION 20 
TRANSVERSE ISOTROPY 
A material that has identical properties in all directions normal to a given 


line is called transversely isotropic. Any line parallel to this line may be 
considered as an axis of transverse isotropy. 


Dielectric Properties, Optical Properties, Thermal Expansion 


With respect to these, a transversely isotropic material behaves as does a 
uniaxial crystal, only two constants being needed to describe each. For 
example, the displacement current in terms of the electric field and the dielec- 
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tric constant matrix is D = = kE where, if x3 is the axis of tran 
tg 


Dy 0 0 
0 Dy 0 


\o 0 Ds 


isotropy: 


D 


Elasticity 


We must find the forms of S and C, (the elastic modulus and elastic ¢«; 


stant matrices) that are not changed by rotations about the axis of trans. 


verse isotropy. We can simplify the work by starting with the crysta 
class that has hexagonal symmetry only. On applying the transformatio: 
S’ = a.Sa = S for arbitrary rotations about X3 we find no further simp)i 


cation follows. Hence the S and C matrices can be copied from those for 


the Hexagonal Pyramidal Class. 
The Piezo-Electric Effect 


Again choosing x3 as the axis of transverse isotropy and starting wit! 
hexagonal symmetry about X3 we find that in order to be invariant to a 
rotations about X3 the matrix must simplify to: 


0 0 0 000 
d={0 0 0 000 
ds, day d3g 0 0 0 


A pitch solidified in an electric field would probably exhibit this kind 0 
piezo electric behaviour. It might also be expected to show an electro opt 
effect governed by a matrix like the conjugate of the above matrix. 

SECTION 21 
APPENDIX 
Transformations 


A counterclockwise rotation of the axes through an angle ¢ about the 


axis is represented by the matrices a and a as follows (where c is written for 


cos @ and s for sin @¢): 


1 00 00 0 
> 3 >) 
wee oe 
a=(0 c s} a= a (21.1 


0-sescc—-s0 0 
0 00 0 ¢ =s 
lo 00 Os e¢ 











elle dioeh ican pit hag io 


abi namie 
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\ counterclockwise rotation about x2 is given by: 


\ counterclockwise rotation about .3 is given by: 


In case one wants only the value of a tensor property in a given directi 
not all the elements of a and a need be used, but only a row or column. A 
special case is that of computing such a property in the direction (@, @ 
of polar coordinates. The Xu axis is chosen in this direction; te and x3 are 
not determined. Writing c; for cos 0, co for cos ¢, s; for sin 6 and Ss» for sin @ 
the required matrices are 


\ 
C; $2 | 
\ | $1 S2 
CiS2°** } 2 
< — =e C2 e 4 
ae = $1 $2 9 a = éi. 
$3 $902 °° ° 
Co 
€1C2S2°°° 


lans?-+<] 


From these the (11) term can be computed for any tensor. 
A few special transformations needed constantly are: A rotation of 180 
ibout X33 ’ 


(1 oo oO od) 
i010 7) 00 
—| 00 
00 1 7) 00 
000 —1 00 
10 0 0 0 —1 0 


10 00 Oo .0 4) 
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A rotation of 90° about xs: 


0100 0 8 
11000 0 O 


QO ft) j 

; M (0010 0 0 sy 
ros F ; @™ Deee-nat or" ~* 

001 


16001 0° 0 
(0000 0 —1 





A reflection in the plane perpendicular to 13: 


(1 00 0 00 
ood 0. 80 
(o> 6 
leq 6 OO 1.0. OO a. 
a= > a= ly 00 a 0 0}? a = a 


00 - 
: oe 0 —1 0] 





000 0 041] 


/ 


A cyclic interchange where x2 replaces x; , etc.: 


The line making equal angles with x, 2, and x3 is a three-fold axis 


(001000) 


100000) 
a 4 } 


oe 010000) . 
a = = , =a 
a <- ° eoueeir ” ~* 
: 1000100 
lo00010) 
A cyclic interchange, where —.3 replaces ., etc.: 


(The line making equal angles with x, x2, and x3 is a six-fold axis 
second sort) 

(010000) 

001000} 

100000 e 
00001 0)’ 

1000001) 

100010 o| 


A rotation of 90° about «3; combined with a reflection in the x3 plane: 


(0 10 00 0) 
100 00 0O| 
001 00 0| [oe 
000 01 QO)’ - 
000-10 0 
(000 00-1 


0 1 0 
a= —1 0 Q |, a= 
00 —!1 





of Lie 


(21.9 


(21.1! 








\ rot 
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station of 60° about 3: 


) 
1 , /3 0| 
tee alk 
V3 1 
—_ 5) 2 * Q@ 
0 01) 


To form a * substitute 
V3 for V3 ina. 


| 


Vv; 


He) Go Hel 


0 


Q 


wa 


\ rotation of 60° about x3 combined with a reflection 


\ ro 


To form a@ 


173 al 
+. 

< v3 1 0| a 
2 2 | 
0 0-1) 


lo form a’ substitute 
~ \/3 for 7/3 in a. 


tation of 120° about 23: 





1 +73.) 

—.- — 0 

2 a 4 

: | 

al R) 3 = 1 0 9 a 

2 2 | 
| 
\ 0 0 1) 


* substitute 
— V3 for V3 ina. 





we! Gn | 


— 


Q 


de! Go oe 


~, 


0 


te et oe 


O 


in. 


bol ke Dol GW 


0 - 


No) GW 
na 


bol GW!) bol G 
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(21.11 


21.12) | 
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Inversion through the origin (a center of symmetry) 


(100000) 
1010000) 
001000; 
a= Q —!1 Q ss | (= (144 

“~~ 1000100| -, tis 
1000010 


loo 0001) 


A rotation of 180° about xy: 


(1000 0 0) 
10100 0 0O| 
0010 0 oO] ., 
0001 0 O| ; 
(0000-1 O} 
(0000 0-1) 


ad = Q —1 


\< 


A reflection in x; plane: 


1000 0 0 


100 0100 0° O| 
0010 0 Of 


a= 010 a= ia =a (21.16 
0001 0 0 
OO 1 H 
0000-1 0| 
0000 0-1) 
A reflection in the x2 plane: 
(100 00. 0) 
1 0O (010 00 0O| 
a 0 ~—10 a= |001 00 Ola'=a=a..(21.17 
OF 000-10 OF; 
000 01 O| 
looo0 00-1) 


In computing the electro-optic and mechanico-optic effects we need 


transformation that will restore to diagonality a matrix that has very sma 


but symmetrical off diagonal terms. 
transformation. 
from an idemfactor. 


. ky + An A Asi 
If k = Ai hog + Aso Avg 
As) Aog k33 + Az; 


This transformation we call a sma 
Such a transformation has its matrix differing but slight; 
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e assume that it can be rediagonalized to the matrix k’ by means of the 
1 2) 01 
transformation: k’ = 6ké. where 6 & t & 
; a 
21.14 
Since 61; 60. + 412 de2 + 413 63 = O we find that, to the first order of small 
quantities 6;; = —6;;. 
Expanding 66, to the first order of small quantities and equating the not 
diagonal terms to zero we find that: 
Aw 2 Avs 1 3 Asi 
a. ; 03 = an 031 . 
Ry — hoo k . ka k ~ k 
1.15 
herefore, to the first order of small quantities A 
ky + An Ap Az) Ry + Ay 0 0 
6 Ay Ro» + Avo Avg 6 0 k T A. () 
Asi Avs kag +. Agz3 0) 0 ksz + A 
| A —As | 
kiy — Roo R33 — kun 
1.16 : 3 — Aj , Avs is the transformation 
where 6 = ; 
Ry = Roo Ro» — k33| x v 21.18 
As) — Ao3 1 
\ Rss ant ky ko» “— sg 
The electro and piezo-optic effects of biaxial crystals can be handled by 
these infinitesimal transformations, but uniaxial crystals and cubic crystals 
" may require finite rotations to re-diagonalize the k © matrix. In the 
ROO 
case of cubic cry stals we note that [0 k A may be diagonalized by 
OAR 
d rotation of 45° about 2), giving 
nN 
\ 
* k 0 0 
OR+A 0 21.19 
0 0 k—A 
k O As k 0 V/A? + A? ; 
upon rotation 
0 k Ay} becomes 0 k 0 a ; 
= 4 through angie 
A, Ay k V Ai + Aj YU k 
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Aj ‘ . ‘ . 
tan © = about x; and diaonalizes by then rotating through 45° at 


giving: 
k+ V/A? + a3 0 0 
k” = 0 k 0 
0 k—~VJ/d?+ a3 


The work can be handled with single column matrices A instead of us 
the square matrices k. If & is a diagonal single column matrix (i.e 
single column matrix of a diagonal matrix), than the almost diag 





matrix A + A is diagonalized by the transformation: 
1 0 0 0 shy 2d | lk +A : 
K,— K,X,-—K\||\" ce 
2A, —2Asz . 
Kz — Kz Ki, — Ko 
—2A, —2A; > 
0 0 1 0 K As 
K “fae K3 Ky — Kz ' + 
PY lg —As As —de — —As |. 
K, — Kz Kz — Ks K, — K. Ki — Ks : 
| —As 0 As i” As 1 Ay nm 
|K, — Kz K, — K3 Ki — Ky Ki — Kz 
|} —Ag Ag As Ay 
—" ' 0 . — . 1 | As 
\K, — Ky Ki — Kz Ki — Ks K2 — Kz JI 
Ki + dy : 
1 ie | 
Ke + Ae 
ee 
Ks; + As 
(21.2 
0 
0 
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transforms vectors as «” = 6x where 


( 1 As As ) 
| n.- kh E-k 
ja i oe 21.22 
K, — Ke Kz — Kz 
Wifes — Ag 1 


(Ay = Ks Ke aa Ks 


If Ai, Ke and Az are not all different the preceding analysis falls down as 
some terms become infinite; a finite transformation is needed in this case 
lhe difficulty can generally be doged by applying a 45° rotation about one 
of the axes. Sometimes the easiest solution is to rotate through an angle 
6 about a coordinate axis then solve for the value of @ that will vanis! 


certain terms. As examples of these devices we give the following: 


(Ki | Ay 
Ky A, + Ag| 
K,;| rotated 45° about x, becomes | A; — Ay (21.23 
| As | | 0 | 
0 | 10 | 
0 (0 } 
K,) ( Ki (K+ Vai + a 
K,, rotated through Ky and ro- |K 
ms tan = about m tated 45" K—-~vVA?+A 21.24 
Ay As 0 about 22 0 
As | x3 is VA? + A?| is 0 
9) : Soe \0 
D = electric induction 
E = electric field 
k = dielectric constant matrix (square) 
K = dielectric constant matrix (single column 
= transformation matrix for vectors 
a = transformation matrix for tensors (of stress tensor sort 
\ = stress matrix (single column) 
= strain matrix (square) 
€ = strain matrix (single column) 
S = elastic modulus matrix 
C = elastic constant matrix 


1 = temperature change of elastic modulus matrix 











~I 
Nm 


1. 
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temperature change of elastic constant matrix 

temperature coefficient of elastic modulus terms 

temperature coefficient of elastic constant terms 

temperature expansion matrix 

piezo-electric constant matrix 

inverse piezo effect matrix 

electro-striction matrix 

electro-optic matrix 

stress optic matrix 

strain optic matrix 

polarization elec. mom. per unit vol. = surface charge per 
area 


susceptibility 


Table of Equations rransformations 

= The form ry 
D= 4x kE lm my 
ex SX y’ ls me Ne 
X = Ce s’ ly m3 ng 
. a l lo l 
es ae H¢ is the transformation a = m, mz m 
H = —S°AS m nen 
Al = tAl + @Bl . 
) = dX A" = ada” 
e= gE D = aD 
g=d, E’ = ak 
K-1 = K-1° +. ZE k' = aka 
Kt me 1? ¥ K’ = ak 
An! = K-77! + pre & = aX 

bk — I] e’ = az'e 
I _ 4 I = 7 ty = al a 

b <" S’ = at Sam! 
1 14 = aha 

4x H = az'Ha™ 


hi; = Ci iT ij 
Hy; = S;;Ts; 
d’ = ada“ 


? 
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The Metallurgy of Fillet Wiped Soldered Joints’* 
By E. E. SCHUMACHER, G. M. BOUTON, G. S. PHIPPS 


HE seriousness of the present tin scarcity has stimulated large con 

sumers of this vital metal to develop drastic conservation measures in 
order to extend the available supplies to cover the emergency period. 
By devising new soldering methods and alloys the Bell System has contrib 
uted a substantial share in the tin conservation effort. Fortunately, the 
changes, as far as can now be determined, have not introduced weakness 
into the soldered joints. Some of the new procedures now used were already 
in the process of development at the onset of the emergency, while others 
were devised under its stress. In some instances, the newly developed 
solders were found to be more difficult to use than the alloys previously 
available, and would not have been introduced under normal conditions 
One major change made that previously had been under consideration will 
result in large tin savings. Unless service difficulties are encountered, this 
modification gives promise to remain after the emergency has passed. The 
change involves a reduction in the amount of solder placed on a wiped joint 
between the cable sheath and the sleeve. Instead of the customary full 
size wiped joint a wipe of fillet preportions is formed. Through this change, 
a solder saving of over 60% per joint can be realized. 

Plumbers and cable splicers have for many years joined lead pipes and 
cable sheath by a soldering process called “wiping.”’ The name is an apt 
description of the operation. In wiping a joint the sections to be united 
are heated by pouring molten solder over their surfaces and manipulating 
the resulting semi-liquid mass by wiping with cloth pads to a well rounded 
symmetrical form such as is shown in Fig. 1. The operation requires con 
siderable skill on the part of the splicer and close control of the solder com 
position. At first consideration, the problem of tightness in such joints 
seems simple but experience shows that even under the best conditions the 
fissures frequently found in the solder occasionally link to form a path that 
allows leakage to occur. In the case of telephone cables not maintained 
under gas pressure, such leaks permit the entrance of water to wet the paper 
covered conductors, thereby impairing the insulation value and causing 
service interruptions. By going to an extreme and wiping off all the solder 
in excess of a fillet, it has been found that many causes of porosity are elim 
inated. Figures 2 and 3 show cross sections of joints wiped the old and new 


*Reprinted from Metals Technology, A.I.M.E., 1943 


i) 














Fig. 1 
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A conventionally wiped joint between telephone cable and sleeve (third siz 





Fig. 3 


A section taken from a joint wiped using the fillet technique (magnification 1} 
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vs. The saving in solder and consequently in strategic tin is evident 
he field splicing forces find that joints are easier to make by the new 
nethod and are less apt to be porous. 

Several interesting metallurgical considerations which are responsible for 
he success of the fillet wipe will now be discussed briefly. Much has been 
written about the wiping process of soldering cable joints and the many re 
quirements of a good wiping solder have been frequently listed. The suc- 
cess of the procedure here described is dependent upon a few fundamental 
characteristics of lead-tin alloys in the process of freezing which have sound 
metallurgical explanations. 

For an understanding of the defects possible in a soldered joint wiped in 
the usual manner, the simple solidification phenomena of metals may be 
considered. As is well known, molten metal in a crucible when allowed to 
cool with free circulation of air will begin freezing near the walls of the vessel 
and with a few exceptions, will end with a concave surface due to solidifica- 
tion shrinkage. Restricting the discussion to a simple lead-tin wiping solder, 
solidification progresses as follows: a lead-tin solid solution commences to 
freeze and forms a rather porous cylinder touching the crucible walls and 
extending to a height corresponding to the volume of the melt at that time: 
on further cooling, dendrites of lead-tin solid solution grow inward toward 
the center of the crucible and at the same time many tiny new crystals form 
throughout the liquid. There are thus taking place simultaneously, shrink- 
age of metal as it becomes solid, shrinkage of previously frozen solid as it 
cools, and shrinkage of the remaining liquid as the temperature drops. The 
originally solidified outer cylinder, adhering to the crucible walls remains 
essentially at its original height. The level of the semi-liquid portion nearer 
the center of the crucible continuously falls until the precipitated crystallites 
formed in the body of the melt make a loosely piled mass extending from 
the upper surface to the bottom of the crucible. Further shrinkage of the 
liquid then leaves these primary crystallites at approximately this level 
while the liquid recedes, leaving fissures between them. This can be beauti- 
fully observed by means of a binocular microscope focussed on the surface 
of a soldifying crucible of wiping solder or, on the top surface of a solidifying 
wiped joint. 

Further insight into the mechanism of wiping solder solidification may 
be gained by another simple illustration. If two solder strips are cast by 
pouring small quantities of molten solder, one on a cold iron surface, and 
the other on a cloth-covered board and both are then bent cold to produce 

specimens as shown in Fig. 4, the chill cast sample will exhibit fewer cracks 
resulting from shrinkage than the slowly cooled one. In the slowly cooled 
sample primary crystallites form throughout the solidifying mass and pack 
at a level above that which the final volume of completely solid solder 














76 BELL SYSTEM TECHNICAL JOURNAL t 

warrants. The sample cast on the cold plate starts to freeze at its 
in contact with the iron plate and, because of the rapid extraction + = 
by the cold iron, it continues to freeze in a rapidly advancing smooth; {1 asl 

until the last liquid at the top is solid. Because of the steep temp : 
gradient there is little opportunity for nucleation and dendrite format ny 
the upper liquid. The surface of this melt is therefore smooth an ni 
from the fissures that are caused by the shrinkage of the eutectic away | ais 

j 

ia 
Bx 

Fig. 4—Bent strips illustrating the effect of variations in cooling rate on the struct 5 

of wiping solders are here shown. The upper strip was chill cast and shows a s j 

ductile surface. The lower strip of the same solder was slowly cooled and upon bending } ; 
exposes fissures between the crystallites at the surface (magnification 3 X). on 
fa 
the dendrites in the slowly cooled sample. Recession of the liquid in th: 4 i] 
slowly cooled sample leaves a multitude of shrinkage channels which, 11 : p 
they occurred at the critical portion of a wiped joint, would cause leaks 3 p 
Another illustration may be useful in demonstrating the processes taking et 
place in connection with joint wiping. Solder may be allowed to solidit: : ( 
in a crucible until its surface is quite firm to a probe. If, at this stage, the | . t 

i 

i 
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cible is tilted sideways to a position shown in Fig. 5 a portion of the r¢ 
naining eutectic may be poured out leaving spongy regions. This loss of 
eutectic is observed frequently during the formation of the old massive type 
oints which may lose several drops by drainage after the splicer has com 
pleted his shaping operations. It is also shown by the greater number of 


pores in the top half of a joint compared to the bottom and the somewhat 


sraver surface appearance of the top. 





Fig. 5—An ingot of wiping solder which had been tilted while in the crucible before 
completely solidified. The lower lip represents eutectic drainage from the partiall 
solidified mass (magnification 1} X). 


Although a solidification range in which quantities of liquid and solid 
metal may exist at equilibrium is an essential feature of a wiping solder, 
another factor of major importance is the nucleation rate of the alloy. Wip 
ing solders having high nucleation rates will develop quickly a myriad of 
points or nuclei throughout the melt from which further crystallization will 
proceed, while an alloy of low nucleation rate will develop relatively few of 
these points in the same time and consequently grow fewer and larger 
crystals. The former alloy will have a texture similar to fine clay while 
the latter will behave like coarse sand and water when subjected to wiping 
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tests. In the fine clay-like texture there are more solid particles pr: 
than in the water-sand type of texture and therefore there is more suriac; 
available for the retention of the liquid in the former type of semi-solj 
mass. Drainage is thus greatly retarded with the result that porosity is 
materially lessened. The type of texture determines in a large measure thy 
ease of shaping and potential porosity of a wiping solder. 

Having examined elementary forms of solidification, attention may nov 
be focussed on the setting up of the wiped joint itself. In practice, the parts 
to be joined are cleaned and fluxed. Circumferentia! paper pasters ar 
then applied to the sheath and sleeve to restrict the spread of the solder 
The splicer then pours hot solder from a ladle over the prepared parts 
catches the excess in a cloth held in contact with the bottom of the joint 
The caught solder is repeatedly pushed back around the cable with a wiping 
motion to aid “tinning”’ or alloying and to distribute the heat. After a {i 
such operations the prepared surfaces can.be seen to be thoroughly wetted 
by the solder. At this stage a portion of the caught solder is mixed in the 
ladle with more hot solder and the mass which now has a clay-like consis! 
ency is poured on the joint and molded into place using cloth pads. When 
solidification has proceeded to a condition where the solder can support it 
self in position, manipulation is stopped. From this point on, loss of heat 
takes place by conduction away from the joint by the sheath and sleeve, }) 
radiation, and by air convection currents at the surface of the solder. As 
a result of this combination of heat losses final solidification takes place in 
the interior of the solder mass near the important sheath-sleeve junctio: 
The action that causes pipes to form in castings draws the eutectic from thi 
critical area between the sheath and the end of the sleeve. If the solder 
has the proper characteristics there will be a shell of solder which does not 
have interconnecting shrinkage cavities, drainage cavities or fissures due to 
the wiping operation and the finished joints will be gas tight. If the solder 
is unduly coarse or has insufficient liquid eutectic at the time the mass is 
too rigid to manipulate further, the resulting joint may leak. 

The new fillet wiping technique is similar to the old up to the step where 
the splicer molds the mass to shape. At this point the new technique co: 
sists in wiping the solder to a small fillet similar to that shown in Fig. 3. 
The resulting joint has much less solder and therefore much less total shrink- 
age and tendency to draw eutectic from the space between the sheath and 
sleeve., Also, at the temperature where wiping is discontinued there is i: 
sufficient volume of solder left by the fillet wiping technique to permit drain 
age drops to accumulate and fall from the bottom of the joint. Therma! 
conduction along the sheath and sleeve cause rapid solidification of th: 
solder at the joint, eliminating the possibility of drainage. Experience has 
shown a consistently high percentage of sound joints when fillet wiping '- 
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rigidly practiced. During the development period of the fillet wiping tech- 
nique examination of the few fillet type wiped joints that were found to leak 
showed quantities of solder present much in excess of that required. Under 
the microscope such joints showed the tell-tale sponginess where the eutectic 
had been drawn away from the junction in the course of final solidification. 

Physical tests on joints made using the fillet wipe between sections of 
telephone cable and sleeving, have demonstrated that fillet joints similar 
in size to that shown in Fig. 2 made with 38°; tin, 0.1% arsenic, balance 
ead wiping solder are stronger in tensile strength, creep and fatigue than 
the cable itself. 

The application of the new technique has gone much further toward 
saving tin than any known permissible change in the composition of solder. 
Using the old technique, a reduction of only one per cent in the nominal 
tin content of a lead-tin wiping solder resulted in widespread occurrence of 
eaky joints, indicating that little tin could be saved by a simple change in 
solder specification. This observation was to be expected since many 
studies had been conducted over the years to reduce the tin content in 
wiping solders to the minimum consistent with the production of satisfac- 
tory joints. Tin has always been much more expensive than lead and for 
arge users of solder a reduction of one per cent in the tin content might 
result in savings of many thousands of dollars annually. 

While the use of the fillet wipe results in large savings in tin other avenues 
for conserving this strategic metal are available such as the substitution of 
ternary and quaternary alloys containing less tin than that required by the 
binary lead-tin wiping solders. A satisfactory alloy of this type was de- 


veloped which contains 13°¢ tin, 23 


> bismuth, 0.1°; arsenic, balance lead. 
Though readily available a short time ago, bismuth now has become too 
restricted to be used extensively in solders. A wiping solder is now being 
introduced into service in which, through the inclusion of a small quantity 
of antimony, it has been possible to reduce the tin content. This material 
appears suitable for fillet wiping although it requires more skill to use than 
the 38% tin, 0.1% arsenic, balance lead wiping solder. Other compositions 
may be usable that contain less than normal tin, but on the whole, the 
savings accomplished by composition modifications will be small compared 
to those produced by the new wiping technique that has been described. 


In SUMMATION 


By virtue of its small solder volume the fillet wipe reduces tin consump- 
tion and produces joints less liable to leakage than the conventional wiped 


joints. The reasons for the success of this type of joint are based on the 


sound metallurgical principles herein described. The use of the fillet wipe 
promises to survive the period of restricted tin consumption. 











A Mathematical Theory of Linear Arrays 
By 
S. A. SCHELKUNOFF 


MATHEMATICAL theory, suitable for appraising and contro 
directive properties of linear antenna arrays, can be based upo: 


simple modification of the usual expression for the radiation intensity oj 
system of radiating sources. The first step in this modification is clos 
analogous to the passage from the representation of instantaneous values 


r 


of harmonically varying quantities by real numbers to a symbolic re; 
sentation of these quantities by complex numbers. The second step co: 
sists in a substitution which identifies the radiation intensity with 
norm’ of a polynomial in a complex variable. The complex variable itsel: 
represents a typical direction in space. This mathematical device permits 
tapping the resources of algebra and leads to a pictorial representation oj 
the radiation intensity. 

An antenna array is a spatial distribution of antennas in which the in- 
dividual antennas are geometrically identical, similarly oriented, an 
energized at similarly situated points. The first and the last properties 
insure that the form of the current distribution is the same in all the ek 
ments of the array and that consequently the array is composed of antennas 
with the same radiation patterns. The difference between individual el: 
ments consists merely in the relative phases and intensities of their radiatio 
fields. The second property means that the radiation patterns of 
individual elements are similarly oriented and that consequently the radic 
tion pattern of the array is the product of the radiation patterns of its typica 
element and the “space factor”. The space factor of an array is detined as 
the radiation pattern of a similar array of non-directive elements. Hence i 
studying the effect of spatial arrangement of antennas, we may contin 
ourselves to non-directive elements and thus materially simplify the analy 
sis. 

An array is /inear if points, similarly situated on the elements, are colinear 
In this paper we are concerned mostly with linear arrays of eguispaced 
sources although in conclusion we shall have an occasion to say a few words 


about more general types. 


1 The norm of a complex number is the square of its absolute value 
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RADIATION INTENSITY AND FIELD STRENGTH 


Consider a linear array of m equispaced nondirective sources (Fig. 1). 
{part from the inverse distance factor, the instantaneous field stre sath of 
the array in the direction making an angle @ with the line of sources may be 


expressed as follows 


\ ©; = Ao cos(wt + do) + A: cos(wt + Y + di) + Ao cos(wt + 2P + Be) 


+ ees + Ano cos(wt + n—2 YP + Bp-2) + coscat + n—-1y), (1) 


y = Blcos 0 — 8, B= —. 








Fig. 1—A linear array of equispaced non-directive sources. If two sources are of equal 
intensity and in phase, their fields at a distant point are substantially equal in 
intensity but differ in phase by 8é cos 6 where f cos @ is the projection of the distance 
between the sources upon the particular spatial direction under consideration. If 
the sources are unequal, an allowance must be made for the relative field inten 
sities in proportion to magnitudes of the sources and the phases must be adjusted 
for the phase difference between the sources. 


In this equation: Ag, A; , --* An = 1 are the relative amplitudes of the 
elements of the array; ? is a progressive phase delay, from left to right, be- 
tween the successive elements of the array; 3), #2, ++ > Jas, Pn. = 0 repre- 
sent the phase deviations from the above progressive phase delay; 8 = 21/X 
is the phase constant, where \ is the wavelength. The radiation intensity, 
that is the power radiated per unit solid angie, is proportional to the square 
of the amplitude of +/ 4; 

Forming another expression similar to (1) but with sines in the place of 
cosines, multiplying the result by i = ~/ —1 and adding it to (1), we have 


1+ Ase? + 
+ Ac.sd” 2y +i n—s 4. e” ad PP 


V ; = [Aye”® a Aievt” 


)) 


lhe true instantaneous value of the field strength is the real part of (2). 
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Hence the amplitude +/¢ of the field strength’ is the absolute va 
))- .3 
(2); thus 


2 -2 -1 
Vb = | ao + ays + aon” + ++ + apo” + 2” |, 


3 
z= e¥ y = Blcos0 — 8, dm = Ame”. 


In this equation: do , a), @2 , *** @n—2, @n-1 = 1 are complex numbers r 
senting the relative amplitudes of the elements of the array and the plas 
deviations of these elements from a given progressive phasing. Th 

all the coefficients are real and positive, they represent the relative a: 
tudes of the elements of the array. If the algebraic sign of a particular 
coefficient is reversed, the phase of the corresponding element is chai 
by 180°; if some coefficient is multiplied by i or —7, the phase of the 
responding element is respectively accelerated or delayed by 90°; an 
general the phase acceleration is equivalent, in our scheme, to a multip 
tion by a unit complex number e*. Some coefficients may be equal to 
and the corresponding elements of the array will be missing. In viey 
this possibility, we shall call ¢ the “apparent” separation between | 
elements; it is the greatest common measure of actual separations. W| 
the elements are equispaced the apparent separation is the actual separat 

Thus we have the fundamental 

Theorem I: Every linear array with commensurable separations betwee 
the elements can be represented by a polynomial and every polynomial can | 
interpreted as a linear array. 

The total length of the array is the product of the apparent separation: 
between the elements and the degree of the polynomial. The degree of th: 
polynomial is one less than the “apparent” number of elements. Thi 
actual number of elements is at most equal to the apparent number. 

The above analytical representation of arrays is accomplished with the 
aid of the following transformation 

$= ¢%, (4 
in which y = 8f cos 6 — # is a function of the angle 6 made by the line oi 
sources with a typical direction in space. Since y is always real, the ab 
solute value of z equals unity and gz itself is always on the circumference 
the unit circle (Fig. 2). As @ increases from 0° (which is in a direction of the 


line of sources) to 180° (which is in the opposite direction), y decreases an! 


2 For brevity’s sake, we shall call ./@ itself the “field strength.” 

3 Equation (3) could be derived directly from the physics of the situation in the sa 
manner as (1). The foregoing method of transition from (1) to (3) serves only the purpos 
of showing the relationship between a less familiar formula and a very well known one 

4If the separations are not commensurable the arrays are represented by an alge 
function with incommensurable exponents. 
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Z=1 





Fig. 2—A typical direction in space is represented by a complex variable which is repre 
sented in a complex plane by a point lying on the circumference of a circle of unit 
radius, having its center at the origin. As the angle @ made by a typical direc 
tion with the line of sources, increases from 0° to 180°, point z moves clockwise 





(A) (8) 


Fig. 3—(A) The active range of z, corresponding to # = 8f and one-quarter wave-length 
separation between the elements. (B) The active range of z, corresponding to 
3 = Bland ¢ = 2n. 


¢ moves in the clockwise direction. When 6 = 0, ~ = Bf — #8; and when 


6 = 180°,y = —Bf— 3. Hence the range y described by z is 
v = 286. (5) 
When the separation ¢ between the successive elements of the array is 


equal to one-half wavelength, the range of z = 2x and as @ varies from 0° 
to 180°, ¢ describes a complete cycle and returns to its original position. 
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In this case there is a one-to-one correspondence between the points o} 
circumference of the unit circle and conical surfaces coaxial with th 
of sources. Such conical surfaces, called radiation cones, are loci of d 
tions in which the radiation intensities are equal. If the separation bet, 
the elements < 4/2, the range of z is smaller than 27 and z describes o1 
portion of the unit circle (Fig. 3A). Finally, if ¢ > A/2, then the path 
overlaps itself (Fig. 3B). Such a path, winding upon itself, will be call: 
Riemann circle. In this instance, one and the same point on the circle 1 
correspond to several radiation cones; but if we regard different positior 
zc along its path as distinct points on the Riemann circle, then there wi 
a one-to-one correspondence between the points on the circle and 
radiation cones. 

Since the radiation intensity is a periodic function of y, the space fa 
of a given array will repeat itself if the separation between the elements 
is greater than one-half wavelength. 


COMPOSITION OF SPACE FACTORS 


Since the product of two polynomials is a polynomial, we obtain the f 
lowing corollary to Theorem I] 

Theorem II: There exists a linear array with a space factor equal to 
product of the space factors of any two linear arrays. 

In other words, there is a linear array such that its radiation intensit, 
in any given direction is the product of the radiation intensities in this dire 
tion of any two givenarrays. Thus we have 


VD, = | do + ays + aoe + ee + ay” 


JB, = | bo + bis + bes te) + Bm as |, 


\ 2; V ®, =|(dtaz+-:°:+ s"-!) (bp + Bis + e+ + Om z"") | 
= | dodo + (aod, + aybo)z + (aobo + ab; + aeby)sz + °-> 


The coefficients of the expanded product represent the amplitudes and the 


phases of the derived array. 
Naturally the process may be repeated and a linear array can be con 


structed with its space factor equal to the product of the space factors of 
any number of linear arrays or to ary power of the space factor of any array. 


For example, let us Start with a pair of equal sources, represented by 


V@=)|1+ 21, 


and construct a linear array with the space factor equal to the square o! 


(7). The field strength of the required array will be 


V@e=|1+2P=|14+224+2\|. 
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\is array consists of three elements with amplitudes proportional to 1, 
If the elements of the original couplet are one-quarter wavelength 

and 90° out of phase, the couplet is ‘‘unidirectional”’. 
factor of such a couplet is depicted by Curve A in Fig. 4. The space 
factor of the triplet represented by (8) is shown by Curve B. In the 


directions in which the couplet radiates half as much or a third as much 


The space 


power as in the principal direction, the triplet radiates correspondingly 
only a quarter or a ninth of the power radiated in the principal direction. 
[he above considerations suggest a simple method for suppressing 


subsidiary radiation lobes. It is well known that in a uniform linear array® 


a t= 2,99 = 


8 


~~ a a“ <a 


4 40° ——S=«G 780° 

a 

Fig. 4—Space Factors—Curve A is the space factor of a unidirectional couplet in which 
{=/4. Curve (B) represents the space factor of an array with amplitudes pro 
portional to 1, 2, 1. 


the difference in levels of the principal maximum of radiation and the first 


subsidiary is substantially independent of the number of elements, pro- 


vided this number is sufficiently large. Thus in the limit, the first sub- 
ordinate maximum is 13.5 decibels below the principal maximum. Con- 
sequently for the array with its space factor equal to the square of the 
space factor of the uniform array, the limiting difference in levels must 
be 26.9 decibels. 


Since the uniform array is represented by 
Vb@=l|14+24+ 24+ ---4 2%], (9) 


°A “uniform” array is an array made up of sources of equal strength with a uniform 
rogressive phase delay. 
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the other array is given by 
Ve = |14+24+ 

1+ 25 + 322° + +++ + mz + (n — 1)2” 


-+2 2e2n 84 ¢ e2n—2 


Thus the amplitudes of the individual sources are proportional to 1, 


n—1,n, n—1,-+--3, 2, 1. Figure 5 depicts the effect of 
‘triangular’ amplitude distribution. 


TA UNIFORM AMPLITUDE DISTRIBUTION | ly 
Mi f analy 
re |cos (2 cos @) | 7 
|? cos($ cos e) i Cons! 
TRI poate AMPLITUDE * eine ataaas f 
19 =| [SIN oo 
\4 cos (2 2 cos ail 


3-6-7-6-3-1 AMPLITUDE DISTRIBUTION | 
- |cos?(2F cos e)| 
Vo =| 
\ 
\|D: BINOMIAL AMPLITUDE DISTRIBUTION 
¥d= sin®(Z cose) 


27 cos*(¥ cos @) | 


Fig. 5—Space Factors—(A) is for a uniform array and (B) for an array with “triangular (Fig. 
amplitude distribution. Fig. 

If the 
A) ¢ 
prope 
Ea 


analy 


Evidently we could raise (9) to any given power 
VO =|Ltstet e-em |m 1 


This process does not change the number of separate radiation lobes. The 
so-called “binomial” distribution of amplitudes was first suggested by 
John Stone Stone.’ His scheme is a special case of (11) if we let » = . 
For the effect of the binomial amplitude distribution see Fig. 5. . == 

The relative merits of two forms for the radiation intensity as given !) throt 
(1) and (3) can now be appraised in the light of the foregoing examples p inth 
Using (1), we have for the instantaneous radiation intensity of the un F chan 
directional couplet 


®6U. S. Patents 1,643,323 and 1,715,433. 





OF LINEAR ARRAYS 


A 


MATHEMATICAL THEOR} 
~ cos 6 — 4 
) 


Vo; cos wt + cos( wt +- 
. T 
cos wt + sn( et + 5 cos s) 


just inspecting this equation, we find no evidence for existence of a 


By 
linear array with a space factor equal to the square of the space factor of 
the couplet. Still less obvious is the method of obtaining proper amplitude 


ratios. 

ARRAYS OF ARRAYS 
The foregoing method of composition of space factors is in reality an 
inalytical expression of geometric construction of “‘arrays of arrays’ 
a pair of equiphase sources of equal strengths 


- 
_ 





Consider, for instance, 
4 
jan 9 ——afes—— f 


- 
(8) 





sean 
(A) 
eg 


~- - , 


maya 
(Cc) 


Fig. 6 





Take two such pairs as elements of an array of the same type 
The middle sources add up to a single source of strength two. 


Fig. 6B). 
\) as elements of (B), then (C) is obtained; the amplitudes of (C) are 


(Fig. 6A). 

If the operation is repeated by taking (B) as elements of (A) or by taking 
fey Ae 

Each shift of a source to the right through distance f is represented 


proportional to 1, 3, 3, 
An algebraic identity 
(13) 


analytically as multiplication by z. 


(ao + ayz + aoz?)z = aos + ayz? + ao2? 


is an expression of an obvious fact that each element of an array is shifted 
Similarly a given change 


through the same distance as the entire array. 
(14) 


in the strength and the phase of the array is achieved by making the same 
change in all its elements; this fact is expressed by the identity 


b(ao + ayz + a2") = day + bayz + daz:?. 
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In general, if an array represented by 
f(z) = ado + ays + a2? + 
is taken as the element of an array given by 
F(z) = bo + diz + dos? +--+ + Oye”), 
then the resulting array of arrays is represented by 
f(s) F(z) = bof(z) + disf(z) + bee*f(z) + ++ + bmias™ f(z). 


DECOMPOSITION THEOREM 
Consider now a pair of non-directive sources with strengths proporti 
to 1, —¢; then 
/* | 
Ve =|z-#l. 
Geometrically, the complex number zs — ¢ is represented by a line draw: 
from point ¢ to point z (Fig. 7A). Accordingly, the radiation intensity 


(A) (8) 


Fig. 7—The radiation intensity of a linear array is represented by the square of the pro« 
of the lines joining the null points of ~/@ to a point z on the unit circle. 


of the pair of sources is represented by the distance between t and z. li yo 
vanishes for some particular direction in space, it vanishes for all dire: 
tions making the same angle with the line of sources; these directions forn 
a cone of silence of the radiation system. Obviously, a radiating couplet 


{ 


has a cone of silence if and only if the null point of ~/@ is in the range of =; 


in particular, there can be no cone of silence unless the null is on the u: 


circle. 

By the fundamental theorem of algebra a polynomial of degree (n — | 
has (w — 1) zeros (some of which may be multiple zeros) and can be fa: 
tored into (zn — 1) binomials; thus 


Vb = | (2 — t)(s — te) «++ (3 — tr-1) 
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Each binomial represents the directive pattern of a pair of elements sepa- 
rated by distance f. Hence 
lheorem III: The space factor of a linear array of n apparent elements 
s the product of the space factors of (n — 1) virtual couplets with their null 
ints at the seros of VU: ti, ta, + tho. 
\ccordingly the radiation intensity of an array is equal to the square 
the product of the distances from the null points of the array to that 
point s on the unit circle which corresponds to the chosen direction (Fig. 
To each null point lying in the range of z, there corresponds one 
ind only ue cone of silence provided each null point is counted as many 


times as s happens to pass it in describing the complete range. 


The null points of a uniform linear array and the point s = 1 representing the 
direction of the greatest radiation divide the unit circle into equal parts. The 
hollow circles represent the null points and the solid circles the points of maximum 
radiation. 


By summing the geometric progression (9) the radiation intensity of 
a uniform array can be represented as follows 


Ve =|-—_ |. (20) 


Hence the null points of such an array are the n-th roots of unity, ex- 
cluding z = 1. Since z is a unit complex number,’ any power of it is also 
iy 


a unit complex number. Moreover, each multiplication by z = e” rep- 
resents a displacement through an arc of y radians. Hence the n-th roots 


” A unit complex number is a complex number whose absolute value is equal to unity 
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of unity divide the circle into m equal parts (Fig. 8). Analytica 
have 
2mri 
z —-1=0, in me * , m= 1,2,3,°--n—1, 
2mr é 2mr ‘: 
Yn = — : COS On = — — : 
n Be nBt 
When z = 1, V@ is evidently a principal maximum. Other m: 


of smaller magnitude, the so-called subordinate or subsidiary maxin 


occur approximately half way between the null points. The g 
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Fig. 9—The field strength ./@ as a function of y for m = 5. The principal maximu 
reduced to unity. 


behavior of the field strength can readily be understood if we follow : 
around the unit circle. When plotted against y, ~/@ has the shape show: 
in Fig. 9. This is a universal radiation characteristic which can be inter 
preted for any particular spacing and phasing between the elements wit! 
the aid of the curve for y + 8 = Bf cos 6 (Fig. 10). 

It is easy to estimate the relative level of the first subordinate maximun 
For a fairly large number of elements, the difference in levels is deter 
mined largely by the distances of the maximum points from the near 
null points. The distances are approximately equal to the circular arcs 
joining the corresponding points. Since the arcs joining the first sul 


ordinate maximum with the nearest null points are nearly half as long as 








first 1 


and t 


maxi! 


li n 


ratio 


levels 





















A MATHEMATICAL THEORY OF LINEAR ARRAYS 91 


those for the principal maximum, the first subordinate maximum of the 
field strength is about one-quarter of the principal maximum. In other 
words, the subordinate maximum is approximately 12 decibels below the 
principal maximum. 











a: ER ee DRE Ae Sr, Se ae ve 





Fig. 10 


A more accurate value for this difference in levels can be obtained by 
first rewriting (20) in the form 


a 
> -= sin v 
s*7—sg ? 2 9 
Vv ® = 3} ae ot = ; y (22) 
sin 5 


ee ; 3m sll 
and then substituting successively y = 0 and y = —, one for the principal 
n 


maximum and the other for the first subordinate. Accordingly we obtain 


+t} / 

. VPRO) _ sect 3m 
39 seas: (23) 
n 

" ‘oe a , 3m : 

Z If m is large sin 7, 8 approximately equal to = and the field strength 

n 2n 
: . 3x ‘i ini : j 
ul ratio becomes = = 4.71. This ratio corresponds to the difference in 


~ 


levels equal to 13.5 decibels. 
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DIRECTIVITY OF ARRAYS 


, 


rr cc = 
he “decomposition theorem’ 


of the preceding section throws co 
able light on directive properties of arrays. The number of elements 
array is one greater than the number of virtual couplets. Hence to 
the greatest possible directivity with a given number of elements, the virt 
couplets must be properly combined. 


Q 
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The null points of several three-element arrays. The spacing betwe 
elements is \/4 and the progressive phase delay is 7/4 (T equals period 
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Fig. 12—Comparison of directivity of several three-element arrays. The spacing betw: 
the elements is \/4; the direction of principal radiation is @ = 0°. Curvs 
refers to the uniform array, (B) to an array with nulls at P and B (see Fig 
and Curve (C) refers to an array with its nulls at B and M. 


For example, the null points of a uniform array of three elements, ont 
quarter wavelength apart, are at Pand Q (Fig. 11). If = 2/2, the rang 
of z consists of the lower half of the unit circle and principal radiation take: 
place in the direction 6 = 0. Evidently, the virtual couplet with its nu 
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at ( is comparatively nondirective. Substituting for this couplet another 
couplet with a null at B should improve the directivity of the array. This 
is indeed the case: In Fig. 12, Curve A depicts directive properties of the 
uniform array and Curve B depicts those of an array with its nulls at P and 


The field strength of the second array is 





Vv b 


=~ = 9 - T 

itsvse +86 *|, == (cos 6 — 1); 

hence the amplitudes of the elements are proportional to 1, 1/3, 1 and the 
total progressive phase delay in the direction of maximum radiation is 


r , W 2r li S 
+ - = — radians. 
2 0 3 


The minor lobe of the second array is substantially smaller than that of 
the first array. The major lobes, however, are equally “wide’” although 
one lobe is somewhat sharper than the other. The width of the major lobe 
can be reduced at the expense of increasing the minor lobe by moving the 
null from P to M (Fig. 11). The effect of this change is shown by Curve 
C (Fig. 12). The corresponding field strength is” 





Ve=l(et+)e4+))=|% + 04+as+i 
.F ir (25) 
=|1-—il+is—if|= 14+ VJ2e *z+e 22; 
hence the amplitudes are proportional to 1, 1/2, 1 and the total progressive 
=e 7 3r 
phase delay is 2 “Beaar i 


For arrays of six elements, one-quarter wavelength apart and with # = 
r/2, we have Fig. 13. Curve A represents the directive characteristic of a 
uniform array, with its nulls as shown in Fig. 14A, and Curve B shows the 
directive properties of an array with its nulls equispaced on the lower half 
of the unit circle as shown in Fig. 14B. 

If the spacing between the elements is f = \/8 and if the phase delay 3 = 
r/2, then the effect of distribution of the null points is even more pronounced 

Figs. 15 and 16). This time z is confined to the fourth quadrant of the 
unit circle. In Fig. 15, 2 = 3; Curve A corresponds to an array with equal 
amplitudes in which case the nulls are equispaced on the complete unit 
circle (Fig. 17A) and Curve B corresponds to an array with its nulls equi- 
. 8 If the “width” of a lobe is measured by the angle of the cone of silence enclosing the 
ode, 

* When transforming the expressions for \/@, it is well to remember that the absolute 


value of a complex quantity does not change if this quantity is multiplied by a unit complex 
number. 
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spaced within the range of z (Fig. 17B). In Fig. 16, 2 = 6; Curve A: 
sents an array with nulls distributed evenly on the complete circl: 
Curve B represents an array with nulls evenly spaced in the range o 
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Fig. 13—Directive properties of 2 six-element arrays with f = \/4. Curve (A) refers t 
a uniform array and Curve (B) refers to an array with its nulls equispa: 
the range of z. 





(A) (B) 


Fig. 14—Disposition of null points for the arrays with directive characteristics as show: 
Fig. 13 


If the total length of an array is kept constant but the number of e|: 
ments is increased, the array may be made more directive; Figure 18 illus 
trates this point. This increase in directivity can be secured only if the nu 
points of the array are properly distributed within the range of z; in Fig. 18 
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Fig. 15—Directive properties of three-element linear arrays with = /8. Curve (A) 
refers to a uniform array and Curve (B) to an array with its nulls equispaced 
in the range of z. 
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Fig. 16—Directive properties of six-element linear arrays with f = /8. Curve (A 
refers to a uniform array and Curve (B) to an array with its nulls equispaced 
in the range of z. 
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ig. 17-—Disposition of nulls for the arrays whose directive properties are shown in Fig. 15. 
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the null points are evenly spaced in the range of 2, appropriate to each sep Ci 
aration between the elements. e the 
If the elements of the array are directive, the null points should be (is affec 
tributed with due reference to the directive pattern of the elements in order 
that a further increase in directivity could be secured. "* 
Let 
poly 
droj 
whe 
Con 
| p ‘ a ——- +s 
ae ~— Sa an i40° a ao 80° 
6 
Fig. 18—Directive properties of linear arrays with total length equal to \/4. (A), whe 
{ = d/4; (B), n = 3, € = d/8; (C), n = 5, € = 0/10. 
If 
MULTI-DIMENSIONAL ARRAYS bite 
The simplest method of constructing multi-dimensional arrays is to take a 
linear array as an element of another linear array. The axis of the second 
array may be chosen to make any angle with the axis of the first array. 
In this way only a special class of multi-dimensional arrays can be formed. yes 
Analytical expressions for the radiation intensities of more general arrays | 2 
can be formulated in terms of two or more complex variables. These | 
variables, however, will not be independent and a given direction in spac ds 
will be represented by a group of related points, one point on each circt bel 
representing the particular complex variable. At this time we shall not be ng 
concerned with any developments applicable to such general multi-dime: dh 
sional arrays. oth 
tha 
; ARRAYS WITH PRESCRIBED SPACE FACTORS tive 
If the minimum separation between the elements does not exceed X/2. ry 
it is theoretically possible to design a linear array with a space factor give 
by an arbitrary function f(y) or F(6) of direction of radiation. Naturally urb 
the number of required elements will be usually infinite; with a finite number = 
of elements the space factcr may only be approximate. 
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Consider an array with an odd number of elements n = 2m + 1. Since 
the modulus of z is unity the polynomial (3) can be divided by 2” without 


S afiecting  ; thus 
Vb= | age7™ aya + agg m2 4+. 
(26) 
T Om T Omyi2 TT °° * TT dem2™ |. 
Let us now assume that the coefficients equidistant from the ends of the 
polynomial are conjugate complex; then the polynomial is real and we can 
drop the bars. Thus setting 


Gm = Ao, Omak = Ap — 1B,y, k> 0, Gu-k = Omik, (27) 
where the 4’s and B’s are real; we have 
Ot ko" + Am_12 R= (A, — 1B,)e** + (4, + 1B, e *¥ 
(28) 
= 2A, cos ky + 28, sin ky. 
Consequently, (26) becomes 
7 V/b = D> «(Ax cos ky + By sin hy), (29) 


k=0 


where ¢, is the Neumann number.” 
If now we wish ~/@ to be a prescribed function f(y) of the variable y 
we need only expand this function in a Fourier series 


0 


Ve = f(y) = bu €x( px cos kW + ge sin ky), (30) 


k=O 


and approximate it with any desired accuracy by means of a finite series 


29). Once the .4’s and B’s are known, we calculate the a’s from (27 





It must be remembered that the real independent variable is not y but 
# and the directive pattern is to be assigned as a function of @. Besides 
being dependent on 6, y is a function of the distance f between the succes- 
he sive elements of the array. Since @ varies from 0° to 180°, the range of y 
is ¥ = 28¢. The function f(y) is prescribed within this range. On the 
other hand the period of the expressions (29) and (30) is 27. This means 
that if Y > 2x, that is if £> 2/2, it is impossible to obtain the desired direc- 
tive pattern with our scheme, because the pattern repeats itself automatically 
) as increases or decreases by 27. But if f</2, we have a considerable 
latitude in the design; outside the range of y, we can supplement f(W) by an 
irbitrary function of y. It is only when f= \/2 that there is a unique class 
: of linear arrays that will produce a directive pattern given by the first 


Me = 1, eq = 2 whenk # O. 
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(m + 1) terms of (30). Dr. T. C. Fry of these Laboratories has sugg 
that leaving ¢ undetermined and fixing the number of elements, an 
could be designed which would have the best fit to the prescribed pat 
In this connection, the “best fit’? means the least mean square devi 
of the approximating pattern from the given pattern. 

If ~/ @ is given as a function F(@) of 6, then by virtue of the definiti 
y we can write 


F(0) = F (cos WV 7 *) = f(y). 


Let us now consider a simple example for the sake of illustrating | 
method. Let f(y) be defined by 


f(y) = 9, O<y<r, 
= |, w<y<2rn. 


We shall assume that the separation between the elements is one-| 
wavelength. This makes the range of y equal to 27. It is also seen | 
regarded as a function of 6, f(W) retains its essential characteristic: b 
equal to zero over one-half of the range and to unity over the remai! 
half. 

Expanding (32) into a Fourier series we have 


ea) . 
> = 


k=1 


1 


b = ) 
fy) ==- ; Ly 


2 
kR—-1 


a) 


n ( 
2 


Consequently 


Ay = 3, Ax = Oif k ¥ 0; 


B, = 0, if k is even; (34 
; are 
B, = ——, if k is odd. 
kr 


Figure 19 shows several approximations to f(y) by means of a finite num 
ber of elements. The curve S,, corresponds to an approximation by th 
finite series (29). If Sy is deemed to be a sufficiently good approximatio: 
to the given directive pattern, then 


g+22 + 


~~ _1\1,1. 
YR * ~ laws 
T / 


1 
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on : : rn r 
The total length of this array is (m — 1) . = 2m==9). All ele 


except the three central ones are separated by one wavelength since th 
powers of z except 2° are missing. 
END-ON ARRAYS WITH EQUISPACED NULL PoINTs 


We now pass to a more detailed analysis of end-on arrays with null p 
equispaced on a given circular arc. 








=| 
z=t3 
Fig. 20 
For an end-on array 3 = Bf and 
z= ev. y = Bl(cos @ — 1). 36 


The range of z begins at z = 1 and extends clockwise to a point determined 


} 


by ¥ = —28f. Let m — 1 null points be equispaced on an arc of length ¥ 


as shown in Fig. 20; the field strength is then 


iv 


Vb = | (z —t)\(s —f)---(¢ — t" yt. t=e ™, (37 
This can be expressed as 
ae J 
@ = 2”"| sin — ) 
v 2 (v T n—1 
(38 


x sin} (+ oy )-+- sin} 
n—1 


“ 
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[he angle of the cone of silence enclosing the major radiation lobe is 
determined from 


Be(cos 5 — 1) = —- . (39) 


7 - Oy y 
— cos 6; = s me : 0 
: _— (n — 1)BC’ ”e 2(n — 1)B¢ (40) 


If the arc W is equal to the range of z, then (40) becomes 
| : 


2 . A 1 
1 — cos A = , sin — = — (41) 
n— 1 2 /n—1 
In this case, the size of the first cone of silence is determined solely by the 
number of elements. On the other hand, if Y = 27 —2z/n, the nulls are 
equispaced on the unit circle and we have an ordinary uniform array; then 


2 r . A r “ 
1 — cos; = a sank sin — = (42) 
np nf 2 4/ Int 


This time the size of the first cone of silence depends upon the total length 
L = (nm — 1)€of the array measured in wavelengths. 

When the number of elements in the first case and the total length of the 
array in the second are large, then we have approximately 


' 2 , d 
6, = , 6 =2 (43) 
wae nl : 


For a large m the ratio of the two cone angles is approximately. 


E ? 
* on 4/%. (44) 
6 N 


For example, if {= /8, the angle of the major lobe in the first case is one-half 
of that in the second case or one-quarter if we are to compare the solid angles. 

Equispacing the null points in the range of s not only makes the major 
lobe narrower but it also makes it sharper. Thus at the point lying halfway 
between the point of maximum radiation and the first null point, the field 
strength relative to the principal maximum is 


‘ : 3 P 5 . (Qn — 3 
sin v sin — v sin v -++ sin W 
4(m — 1) 4(n — 1) 4(n — 1) 4(n — 1) - 
x= = ae —s = « (4a) 
y , 2y ; 3y . (n — 1) 


sin — sin sin -++ sin 
2(m — 1) 2(n — 1) 2(n — 1) 2(n — 1) 
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For a quarter wavelength separation between the elements Y = 7 and 


ratio is equal to 


1 
x — 
V/2in — 1) 


so that the drop in the radiation intensity becomes [10 Logio (nm — 1) 
decibels. On the other hand, for a long uniform array the correspond 
drop is independent of ” and is equal to 4 decibels. 

Another consequence of equispacing the null points in the range o 
consists in substantial suppression of subsidiary radiation lobes. The tirst 
subordinate maximum is situated approximately halfway between the first 


3 : 
two null points where y = — : id 1) ; thus the field strength there, relativ: 
2(n — 
to the principal maximum, is 
a ee se . (2n — Sv 
sin — sin +++ sin ~— 
"“He—-1) 4 — , 4(n — 1) 4(n — 1) (a7 
x = = = = go (4 
, y ; 2y . 3y . (n — 1)p 
sin sin sin +++ sin 
~ Ste = 9) 2(n — 1) 2(n — 1) fig a 1) 


For a quarter wavelength separation this field strength becomes 


-—— . T 
. wen 4(n — 1) 
aie = t= oe 
feces ve e — 
V 2(n — 1) 2 sin ach 
(4s 
sin c 
= 4(n — 1) 
. 2n—3)r 
a(n — 
/2(n 1) sin din = 1) 
When » is sufficiently large, we have approximately 
en = 49 
x &n —1)V2(n — 1) ( 


and the-subsidiary maximum is [30 Logio (n — 1) + 5] decibels below the 
principal maximum. Each time the number of elements is doubled, the 
level of the subsidiary maximum is diminished by about 9 decibels. Thus 
an array of the type (37) with ~ = 26¢ has very sharp directive properties 
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in order to find the relative amplitudes and phase deviations of the ele- 
ents of the array represented by (37), we expand ¥ # into a single poly- 
nomial as follows" 


Jb = |(1 — f'2)(A — £2) --- (A — £""'s) | 


n—1l n+1 +2 n+k A +1) 
\1—-? )---A—¢ ) - k 
= /! (—) - = €£ 
+ 2 ERE ECE 
Ph. oes. n—l _ nk n—k 
n—1 2 t 2 ) ° (t 2 _ t 2 ) _kn 
=|1+ 2 (-))—71—3 Ea EE (50) 
ie (@?-—P#)(¢' -—t---@?-P) 
. (n — 1)p (n — k)yp 
= — Cs din — 1 kn 
=|1+2, (—)* ——~— — 2(n=1) ¢! 
x (-) . # ee We, 
sin -+* sin 
2(n — 1) 2(n — 1) 


Hence the progressive phase delay from one antenna to the next is equal to 


v=r— — ¥ and the amplitudes are in the ratio 
2(n — 1) 
. (n — 1) . (n— 1b... (n — 2)p 
ti 2(n — 1) ria 2(n — 1) — 2(m — 1) 
es «aa * ie 
sin 5 1) sin 55 1) sin 57 — 1) 
. (51) 
(n — Ie (n—2)yp. os (n — 3)p 
2(" — 1) 2(n—1) = 2(n — 1) 


, y me 2y . sy” 
sin 6 — 1) sin 2 — 1) sin 28 — 1) 


The amplitudes of the elements equidistant from the ends of the array are 
equal. In the special case of an end-on array with nulls equispaced in the 
range of s, ¥ = 28¢ and 8 = Bf; hence the progressive phase delay from one 
Be 
1—1- 
While (50) serves well for finding the amplitude and phase distribution in 
the individual elements of the array, another form is more general for cal- 
culating the directive properties. In order to obtain this form we set 


antenna to the next is  — 


4 Chrystal’s Algebra, Vol. 2, p. 340, (1926). 








where 








D = 2 cos 


(n 


Let D be the maximum value of /; then the gain of the array over 


single source is given by 


i 
6 «ies. a 


® dQ | ® sin 6 d6 
0 


where © is the solid angle and the integration is extended over a unit sphere 
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. (n— 1p . n— kp Since 
sin ‘++ sin 
2(n — 1) 2(n — 1) ter, | 
po = 1, Pn-1-k = Pr, Pr = 3 
sin v +++ Sin ai ® 
. ) , ee . ) 
2(n — 1) 2(m — 1) 
7 a 
v v 
P= + sai 
v Z 2(n — 1) . 
i(n—l)¢ 
divide the last expression in (50) bye ? and combine the terms « 
distant from the ends. Thus we obtain 
( 1 ( " Sinc 
n— 1) n— 3)y9 
V/V = 2 Cos ; + 2p; cos : : 
(n = 5)e 
. 
+ 2p2 cos ; + . 
. D up ° of ° 
where the last term is 2p»_, cos > if mis even and p»—1 if mis odd. we 


2p* 
= 10 Logy — decibels, (5 


be 
For an end-on array with nulls equispaced in the range of z, the maximum | 
radiation is in the direction Y = 0. Thus we shall have 


_ 1 Po ( om 3) D 
D ha + 2p: cos es ? es 
: (55 
-_—= o A 
+ 2p. cos ~~ ; Jeo ae 
- Su 
- un 
go = B+ ——_ +e. (56) 
= 1 9 


A corivenient expression for the radiation intensity can be obtained from 
(50) by taking its norm 


o = [Po + pie’? + pre?'¥ + eos + Pn—e'” 1 ¢] 


‘(Po 4. pre *? + aw + Pn ye7 tinh ¢]. 
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Since the set of coefficients po , ~1, p2 *** Pn-1 is symmetric about the cen- 
ter, we find 
® 2pi cos (n — 1)¢ + 2(popi + pipo) cos (n — 2g 
+ 2(pope + pipi + propo) cos (n — 3)¢ 
+ 2( pops + pipo + popi + pspo) cos (n — 4)o4+ --- <9 
+ 2(Popn-2 + Pipn-s + prpns+ + -+* + Pn—2po) cos ¢ 
+ (poPn1+ Pipn-2 + poPn-s + °° + Pn-ipo). 


Since 


| ® sin 0d0 = Ll [ @ dg, 
0 Bed, 


“ ¢1 


(59) 


y 


: T Qi = gvo- 28¢ 
2(n — 1) eh ; 


go = pl-8+—-+ 
we can write 


[o sin 0d0 = 1 [ 2p0 sin (» — 1)¢ 4 2(popi + pi po) sin (nm — 2)p 
ai pe n—1 n— 72 
(60) 


€o 


+ +++ + (popa-a + Pipaa2 + °°* + Parpoy - 


¥1 


For an end-on array with nulls equispaced in the range of z, (60) becomes 


rr 2 [2(-)" "pe .. — 1)8¢ 
® sin 6d9 = y? sin (n — 1)8f cos (n 1) 
“0 pt oa 9 ad 
oun n—2 5) a 2) 
” Am) a Pio) sin (n — 2)8¢ cos = eS eins 1 
Se n— 


+ (popna-1 + pipn-2 + +°* + pa-ip 3¢]. 


Substituting in (54), we shall obtain the gain of the array. 
Similar expressions can be obtained for an end-on array in which the 
amplitudes of the individual elements are equal. Thus we have 
1, n-ne ay , 
(n=1) i(n—2)¥ i 
$=—(er" +6 +e te" +1] 


n° 


7 ea oh So | 


Il 


“ [2 cos (n — 1)¥ + 4 cos (nm — 2)p 


+ 6cos (n — Z)\p + +++ + 2(n — 1) cosy + al, 
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Fig. 21—The directive gain in decibels of a pair of sources with equal amplitudes. (A thre 
the phase delay between the sources is 2xf/\; (B), the phase delay between t! rang 
sources is  — 2nf/X. \ 
end- 
tion 
no i 
10} | mit! 
| can! 
the 
‘ Im] 
a | ant 
VU 
3 | ttt tt larg 
Z5) + + +——_—_+__4 + + + + 4 + + + + + + — + fam” ] 
Z| dec 
0 t + + + + + -——+ > + 
] | 
| , 4 
| 
| | | | 
Bon 
. oS ae ae See ae i Se Eo 
0 005 0.10 O68 020 O25 
SEPARATION IN WAVELENGTHS 
Fig. 22—The gain as a function of separation in wavelengths: n is the number of elements 
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where y = BE (cos@— 1). In this case D = 1 and 


sin 2(mn — 1)Bf , 2 sin 2(m — 2)p¢ 
4+ 
n—1 n— 2 


ae 2 
| ® sin 6d0 = atl nl + 
(63) 


3 sin 2(n — 3)Bf 
sin 2(n — 3)80 


+ (n — 1) sin 26¢ 
n— 3 


When the separation between the elements is exactly an integral number 
of quarter wavelengths, (63) becomes 


al ? 
[ eS ee (64) 
0 


n 


and consequently the gain is 
G = 10 Logy n. (65) 


Figure 21 contrasts the directive gain of a pair of sources of equal strength 
with the phase delay 2rf/X (Curve A) with a directive gain of another pair 
of sources of equal strength but with the phase delay  — 2rf/X (Curve B). 
In one case the directive gain diminishes with separation between the ele- 
ments and in the other it increases. Figure 22 shows the directive gain of 
three-element and four-element end-on arrays with nulls equispaced in the 
range of z. 

As the separation between the elements decreases, the directive gain of an 
end-on array with nulls equispaced in the range of z increases but the radia- 
tion intensity per ampere-meter decreases. This circumstance would be of 
no importance if we had perfect conductors at our disposal to make trans- 
mitting and receiving antennas; but in reality parasitic losses in themselves 
cannot be removed and the efficiency of an array decreases, therefore, with 
the separation between the elements. This decrease in efficiency will 
impose an upper limit on the overall gain that can be obtained with small 
antenna arrays in spite of the fact that the directive gain could be made very 
large. 

Likewise the band width diminishes as the distance between the elements 
decreases. This imposes another limitation on arrays of this type. 














Memorial to the Classical Statistics 
By KARL K. DARROW 


NE of the most elusive and perplexing, hazy and confusing of the parts 

of theoretical physics is that which bears the name of “statist 
mechanics”’.* On the principle that a tree is to be judged by its fruits, | 
must be ranked as high as the tree which bore the golden apples of the 
Hesperides; for among its fruits are the Maxwell-Boltzmann distributi 


law, the black-body radiation law, the value of the chemical constant, the 
Fermi distribution-law for the electrons in metals, the alternating intensities 
in band-spectra—and indeed the tree might lay a valid claim to the whole oi 
quantum-theory. The singular thing is that such wonderful fruits should 
have grown from, or should have been grafted upon, so badly-rooted a tree. 
To change the metaphor, one frequently feels that the superstructure is 
sustaining the foundations, and the premises are flowing from the con 

sequences, rather than the other way about. Perhaps anyone who feels 
this way should be disqualified from writing about the subject; but on the 
present occasion, the attempt is going to be made. 

Statistical mechanics—hereinafter to be called “‘S.M.” at times for short 
did not of course arise from any desire to solve the problems suggested above, 
which came late. It seems to have sprung from attempts to answer older 
questions, of which the following may serve as an example. Consider a gas 
in a box, with an electric fan or something of the sort fitted inside to stir 
it up. The gas having been stirred up, the fan is stopped, leaving it in a 
state of surging and whirling about within the confines of the box. Very 
shortly, however, the surging and the whirling cease, the gas having passed 
of itself into a state of tranquillity and uniformity—uniform density, uni 
form pressure, uniform temperature. From this state it never departs, un 
less stirred up afresh. There is a tendency of the gas to go of itself from the 
state of surging into the state of uniformity, and no tendency at all for it to 
go from the state of uniformity back into the state of surging. This is very 
unlike the behavior of a pendulum, which having fallen from one end of the 
arc of its sweep to the middle thereof, moves on to the opposite end, re- 
traces its path and returns to its first situation. Why should the gas behav: 
that other way? 

*T acknowledge with gratitude the incentive given me by Smith College to explot 
this subject, by offering me the opportunity of giving a course on statistical and chen 
cal physics in the spring semester of 1942. 
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For an answer to this question and others of the kind, S.M. offers the 
following statement: 


Basic THEOREM OF STATISTICAL MECHANICS 


{ system is more likely to be found in a state of greater probability than 
in a State of lesser probability 


It may be that no reader of these lines has ever seen the basic theorem of 
S.M. set forth with such merciless candor, though in many a sober treatise 
there is an elaborate statement which when analyzed turns out to be just 
this and nothing more. Of course it is a tautological statement, and has no 
value except insofar as it may help to drive some contradictory notion out 


of the student’s mind or to prepare that mind for some meaning or other 
which is not yet in the statement but may be added to it later. Actually it 
can serve both these offices. 

l'o be expelled from the mind of the student is first of all the idea that 
S.M. is going to give him a description of the way in which the gas proceeds 
from the surging state to the uniform. From an astronomer he may 
earn the orbit of the moon from apogee to perigee. From an authority 
on ballistics he may find the trajectory of the bullet from the muzzle of the 
gun to the bull’s eye on the target. From a railroad office he may get a 
timetable showing the passage of the train from mile to mile over the rails 
from Boston to Chicago. All this sort of thing is out of the range of 
" statistical mechanics! If a railroad acted like a surging gas and its time- 
table were devised in the spirit of S.M., one would go to the office and be 
told that the trains were enormously more likely to be in Boston than in 
Chicago or anywhere in between. From this one would be expected to infer 
that at any moment chosen at random the chance of finding a train anywhere 
along the line except in Boston would be practically nil—unless indeed one 
got a train and put it on the rails at Springfield, and even this would be of 
little use for getting to Chicago, since at every subsequent instant the 
train would almost certainly be in Boston. Not a very useful timetable, 
and not a very useful railroad! 

S.M. thus starts off with a renunciation. It renounces the prospect of 
telling just how the gas proceeds from the surging state to the uniform state. 
lo that smooth unbroken sequence of times and places whereby the moon 
finds its way through the heavens and the bullet through the air and the 
train along the rails, there is no counterpart presented. 

This of course is a serious matter, for the smooth unbroken sequence is 
inseparably linked—or almost inseparably linked—with the notion of cause- 
and-effect, the notion of natural law, the notion of man as a being who can 
foretell the future. Mechanics harmonizes with these notions; for mech- 
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anics is the science which professes that, given the positions and m« 
and the forces in a system of particles at 10 A.M. sharp, it can pred 
positions and momenta of all the particles at 11 A.M. and every inst 
between and all through the endless future. Statistical mechanics, 
the implications of its name, is nowhere nearly so audacious. 

Suppose the electric fan, or whatever stirring-gadget was employe 
stopped an instant before 10 A.M. sharp. S.M. limits itself to afi 
that at 11 A.M. the most probable state for the gas in the container 
immensely most probable state, the almost-certain state—is the ur: 
state. It also says the same thing exactly for 10:15 A.M., and for 1 
A.M., and even for 10 A.M. sharp. If at 10 A.M. sharp the gas is in a 
of wild and furious surging, S.M. does not deny the fact, but sees no r 
for revising its own affirmation. If at 10:01 A.M. the gas is settling 
but has not yet quite reached the uniform state, that again does not deter 
S.M. from standing by its assertion. Whenever a freak of chance or a 
man may have produced one of the states which it calls improbable, S.M. 
just says ‘“‘wait, and you shall have the state which / am going to talk about 
To further questions it can only say “I know my limits”—and that is y 
its basic theorem says for it. 

If now the negative aspect of the basic theorem is sufficiently clear, we | 
address ourselves to the task of giving the theorem a positive meaning. For 
this there is but one way: the word “probability” must be replaced wit 
some word or phrase or mathematical expression which does have a meani 
After this is done we can of course restore the word “probability” as 
equivalent for that other word or expression. The basic theorem will 1 
be tautological upon the surface only, for actually it will have the meaning 
conferred upon it by the definition of its key-word. 

Various meanings have been offered for the key-word, by various peop 
who have been successful in getting useful results out of statistical mechanics 
Until 1924 the dominant meaning was that imposed by Gibbs and Bolt 
mann. From this meaning arises the form of S.M. which is called “‘the classi- 
cal statistics”. (The word “statistics”, by the way, is a bad but commo! 
abbreviation for “statistical mechanics”.) This is the topic of the present 
article. In 1924 there was proposed a novel meaning for the key-word, 
which led to results sometimes agreeing with, sometimes differing from, thos 
attained by the classical statistics. Where the results of the two agreed 
they agreed with experiment also; where the results of the two disagre« 
experiment sustained the new one. This event has left the classical statis 
tics in a strange situation, in which one cannot exclude the possibility t! 
all of its remarkable achievement is due to a happy but deceptive chanc: 


The classical statistics may indeed be only a past episode in the history o/ 


scientific thought, and it is for this reason that I have given to the article | 
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strange and sombre title “‘Memorial to the Classical Statistics’. Yet even 
1s a past episode, it is worthy of remembrance; its didactic value may yet be 
reat; and perhaps the human mind may some day stretch its powers to the 
point of conceiving the classical and the new statistics as aspects of a single 

hole, as it has lately stretched itself to the extent of uniting the wave- 
picture and the corpuscular picture of matter and of light. 


THE MAXWELL STATISTICS 


Since the main concern of S.M. is with the ‘most probable state’’, one 
sees that its principal content must be made up of assertions about that most 
probable state. Maxwell made such an assertion. He wrote down a for- 
mula for the distribution-in-velocity of the molecules of a gas. It is the 
formula now called ‘‘the Maxwell-Boltzmann distribution-law”’, which is so 
well known to the readers of this journal that I will not bother to write it down 
intil there is actual need for having it on the page. Maxwell might have 
said bluntly: “This is the distribution which I will assume for the most 
probable state”; and having said so, left it at that. He did not leave it at 
that, and presumably he would have been dissatisfied so to leave it, as most 
of us would be. Instead, he postulated a pair of attributes for the most 
probable state, and showed that if these are the attributes, then the distri- 
bution is according to that formula. 

rhe attributes which Maxwell postulated are “isotropy”’ and ‘“independ- 
ence’. 

The former is easy enough. One assumes that in the most probable state, 
the distribution of velocities of the molecules is isotropic. Nothing can 
usefully be added to this simple statement. 

rhe latter is a little harder to grasp. Perhaps it can best be exhibited by 
describing a couple of imagined cases for which it would not be valid. Sup- 
pose for instance that all of the molecules have the same speed-—the same 
magnitude of velocity, though their velocity-vectors be pointed in all 
directions. Let this common value of speed be denoted by V, and let any 
direction chosen at random be made the axis of x in an ordinary coordinate- 
frame. If a molecule happened to be travelling with such a velocity that 
the component thereof along the x-axis, v, let us call it, was just equal to V, 
then it would be a certainty that v, and v,, the y and z components of the 
velocity, were both of them zero. If a molecule happened to be travelling 
in such a way that v, was zero, then either 2, or v, or both of them would have 
to be different from zero, and the square root of the sum of the squares of 
v, and v,; would have to be equal to V. There would consequently be a 
correlation between the values of the three components, and the probable 
nay even the possible—values of any one of them would be affected by those 


of the other two. If the molecules had a uniform distribution of speeds up 
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to a maximum value of V, there would still be a correlation of a 

sort, though not so marked a one: the higher the velocity-component 
x-direction, the lower would the y and the s components be likely 
One could imagine distributions for which the higher the velocity-com; 
in the x-direction, the higher would the y and the s components bx 
to be. 

The “assumption of independence” is, that in the most probabl 
there is no correlation at all. Whether the x-component of the velocit 
molecule is high or low is a detail which has no influence whatever « 
possible or the probable values of the y and the s components. Low \ 
of v, go just as well and just as abundantly with low values of 2, as with 
and reversely. 

The Maxwell-Boltzmann law, as I said, is the distribution-law 
conforms to both the assumption of isotropy and the assumption of i) 
pendence. So the question arises: do those two assumptions havi 
quality of plausibility and of convincingness, which make the average pe: 
son say “‘Surely these must be the attributes of the most probable state of 
gas!” I do not know what result a referendum on this question would give, 
but it is my guess that most physicists would feel more satisfied with thes 
than they would with the Maxwell-Boltzmann distribution-law if it wer 
tossed out to them with the bare affirmation “This is assumed to be 
attribute of the most probable state”. Clearly this is how Maxwell felt, ai 
there is no better guide than the intuition of a Maxwell. 

The foregoing question is something else than the question whether +! 
assumptions, and the Maxwell-Boltzmann distribution-law which follows 
from them, are truly the attributes of the most probable state. It is a strang 
historical fact that not for many years after the promulgation of this famous 
law, and not till after both of its sponsors were dead, was there any proper 
test of it. The derivations of the law were exercises in abstract and u 
renumerated thought. Nevertheless experiment—applied to thermioni 
electrons, to molecules of ordinary gases, to thermal neutrons—came at long 
last to justify Maxwell. To any who may feel that the assumption 
independence is in itself too reasonable to require any proof, I disclose that 
in other forms of statistics this assumption is declared to be false, except as 
an approximation. 

The “Maxwell statistics” therefore consists in the main of the statement: 

The most probable state of a gas is that in which isotropy and independenc: 
prevail among the velocity-vectors of the molecules. 

We now require some terminology and some notation. 

I take for granted an understanding of the terms “‘velocity-vector” and 
“‘distribution-in-velocity”, these being learned by physicists out of kinetii 
theory if not out of S.M. A velocity-vector may be replaced by a point 
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which serves just as well for all of its purposes andevenbetterforsome. Let 
the velocity-components 2, , v, ,and 2, be laid out along the axes of a Cartesian 
coordinate-frame, and the vector for any molecule be drawn from the origin: 
the point at its tip is the point in question. Point and coordinate-frame are 
said to be “in velocity-space’’. Statistical mechanics prefers as a rule to 
deal with the momenta of the molecules rather than their velocities. This 
is for valid and powerful reasons, one of which is that the transition to the 
case of photons becomes much easier.'! In the case of material gases it 
makes no practical difference, since the momentum of a molecule is its 
velocity-vector multiplied by the mass of the molecule which is practically a 
constant, and every statement about the distribution-in-velocity can with 
the utmost ease be translated into a statement about the distribution-in- 
momentum and vice versa. The momentum-vector may be replaced by the 
point at its tip, having coordinates p,, p, and p, in a coordinate-frame in 
“momentum-space’”’. If we consider together with these the three co- 
ordinates x, y, 3 of the molecule in ordinary space, we may say that we are 
locating the molecule in six-dimensional space. I have yet to meet someone 
who claims that he can visualize a six-dimensional space, and yet there is no 
doubt that the phrase fulfills a psychological need and has a practical value. 
The six-dimensional space of these particular six variables is called “the 
u-space”’. 

It seems odd to bring in the u-space before considering by itself the three- 
dimensional “‘ordinary”’ or ‘“‘coordinate-space” in which the gas is located. 
Is there nothing to be said about the most probable distribution of the 
molecules in the coordinate-space? Well, “every schoolboy knows” that 
the state to which a gas tends and in which it remains is a state of uniform 
density. Maxwell, I think, accepted this as one of the facts behind which 
one cannot, or does not, go. For a complete statement of the Maxwell 
statistics I therefore offer the following: 

A gas is very much more likely to be in its “most probable state” than in any 
other. The most probable state is that in which isotropy and independence pre- 
vail among the momentum-vectors, while the distribution in coordinate-s pace is 
uniform, 

So in the Maxwell statistics the distribution-in-momentum of the mole- 
cules is derived from assumptions ostensibly more basic, while the distri- 
bution-in-ordinary-space is simply affirmed. If a theory could be devised in 
which both were derived from assumptions apparently more basic, one would 
be likely to feel that something had been gained. Now this is a char- 
acteristic, and one of the principal virtues, of Boltzmann’s theory known as 
the “Boltzmann Statistics” or as the ‘Classical Statistics’. 


‘Another reason has to do with “Liouville’s theorem,” for which unfortunately I 
cannot make room without overloading this article. 
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THe BOLTZMANN STATISTICS 

Boltzmann invented a way of appraising the probability of any im 
state of a gas, which has the following very remarkable features: 

(a) It gives so sharp a definition to the key-word “probability”, th 
only can the state of maximum probability be identified, but the ra: 
the probabilities of any two states can be computed. 

(6) For the distribution-in-momentum of the molecules in the mos 
probable state, it derives a formula identical with that which springs 
the Maxwell statistics. This of course is why the formula is known as t} 
Maxwell-Boltzmann law. 

(c) For the distribution-in-space of the molecules in the most probal 
state, it derives the uniform distribution. 

All this does not entail that the Boltzmann statistics is necessarily righ 
It does, however, lead to consequences, which it is the privilege and th 
affair of experimental physics to verify or to reject. 

I can now write down a phrase into which the Boltzmann statistics, 
equally well those which came later, can be fitted: 

The probability of a state is the number of different ways in which the stat: 
be realized. 

This is another of those oracular sayings which acquire a meaning 01 
after some meaning is given to the key-word, which is this case is ways. | 
could now rewrite the basi¢ theorem without the word “probability”, 
so can the reader; but the only effect would be to transfer the mystery out oi 
the word “probability” and into the word “way”. Boltzmann, however 
assigned a meaning to the latter word. It is this meaning which we nov 
must strive to realize. 

For this purpose I propose a game of which the outfit consists of a sack 
an enormous number NV of balls, and a smaller number M of baskets. Th 
game is played by reaching into the sack, drawing out the balls one after 
another, and tossing them into the baskets. All of the balls feel precise!) 
alike to the hand, so that there is never the least inclination to put one asic 
and pick up another as one’s hand gropes around in the sack. Nevertheless 
when one looks at the balls after they have fallen into the baskets, one sees 
that they are nicely adorned with the integer numbers running from 1 to .\ 
Incidentally the baskets also are numbered. It is this numbering whi: 


gives point to the game. 


Someone or other—someone who might be designated as the caller, after 


the man who calls the figures of a square-dance—has prescribed a sequenc 
of M numbers \, and N»2 and N; and soon to Ny, all of them positive 


integers and totalling up to NV. A single inning of the game consists in 


drawing all of the balls out of the sack one after another, and dropping th 
first V, which come out into the basket I, the next V2 which emerge into the 
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basket II, and so on until every one of the balls is reposing in one or another 
basket. Now along comes the umpire with pencil and pen, and he writes 


own on one sheet of his pad of paper the numbers of all the balls which are 
n basket I, and on a second sheet the numbers of all which are in II, and so 
on until he has got an inventory of the contents of all of the baskets. The 
noventory does not state the order in which the balls in any basket were dropped 
into that basket. ‘That order is blotted out and forgotten. The inventory 
tates which balls are in which baskets, and lets it go at that. 

This does not seem a very entertaining game, but entertainment is not 
what it is for. The present question is: how many different inventories can 
there be, consistent with that sequence of figures V,, V2,.N3, «+: Vy which 
the caller prescribed at the start? 

The answer is obtained in what must seem, to anyone meeting for the 
first time such a question, a strangely devious way. 

First we evaluate the whole number of different orders in which the balls 
can be drawn from the sack. This is \-factorial or V!; for the first ball to 


emerge may be any one of the .V, and the next may then be any one of the 


\Y—1) remaining, and the next may then be any one of the (VW —2) remain- 
ing, and so on to the end. 

If each order corresponded to a different inventory, V! would be our 
answer. Clearly this is so, if and only if there are as many baskets as balls 
and one ball in every basket. In all other cases V! is larger, and often 
colossally larger, than the number which we seek. It is necessary now to see 
that this great multitude of V! different orders falls into groups composed of 
Y orders apiece, all of those in a single group corresponding to a single 
inventory—necessary to see this, and to calculate X; whereupon we shall 
find that X, the ‘‘number of orders per inventory”’, is the same for all of the 
inventories—so that the number which we seek is V! divided by thiscommon 
value of X,. 

It seems to be helpful to think of some one inventory, and of some one 
order which leads to that inventory. By a certain amount of mental effort, 
which varies from person to person, it can be seen that this particular order 
is but one among N,! No! N3! +--+ Nyy! different orders all leading to the 
very same inventory. For think of the \V, numbered balis which lie in the 
lirst of the baskets: there are \,! different orders in which they could have 


/come out of the sack, and every one of these corresponds to the very same 


inventory. Think next of the V2 numbered balls which rest in the second 
basket: they might have come out of the sack in V,! different orders, without 
changing the inventory. Think now of the contents of both of these baskets 
at once. Each of the N2! orders in which the second basketful may come 


} out of the sack may follow on any one of the \,! orders in which it is possible 
for the first basketful to emerge. 


The product .V;! V2! is therefore the total 
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number of ways in which the first (NV, + Ne) of the balls might have 
out of the sack without changing the inventory. 

The process of proof need not be carried further. X has been eval 
It bears no earmark of whatever particular inventory the student may hay 
chosen to adopt at the beginning. It depends only upon the seque: 
numbers V,, Ne, ---, Vy fixed by the caller, which sequence I will hereafter 
term a “‘distribution”’. 

The number of inventories—or ‘“‘complexions”’, to use a commoner w 


for the distribution V,, ---, Vy is therefore given by the formula, 
W = N!/N,! No! «++ Nu! = N!/ILN;! 1 


The theorem to which we are advancing affirms that this number has its 
maximum value for the uniform distribution—the distribution in which th 
caller assigns the same number of balls, V/M, to each of the baskets. 

The usual argument for this statement may be put as follows: Let 
assume the uniform distribution, with A = (V/M) balls in each basket, and 
compare its value of W with that of one of the neighboring distributions 
such as the one in which there are (A + 1) balls in the first of the baskets 
(A — 1) in the second and A in each of the rest. It is not even necessary 
to get out a pencil and paper to see that W for the latter is less than W’ for 
the former, being in fact just 4/(A + 1) times as great. The same is 
evidently true for disarrangements of the uniform distribution which involv: 
more than two baskets and more than one ball per basket. The conclusio 
is clinched by the obvious fact that when all of the balls are in any one bas- 
ket, W has its least possible value, viz. unity. (To unite this formally wit! 
the previous statements, one must follow the mathematicians’ practice of 
using a symbol 0! or “‘zero-factorial”’ and giving it the value unity). 

We shall have to play this not so very entertaining game on several oc- 
casions in S.M., altering the meaning of the balls and the meaning of the 
baskets from one occasion to the next. The reader has probably guessed 
that the balls stand for the molecules. The guess is right in the classica 
statistics, wrong in the newer forms. To get at the meaning of the baskets 
suppose the gas contained in a box of volume V, the interior of which is 
divided up by impalpable coordinate-planes into compartments or cells a! 
of the same volume V). The baskets stand for the cells. 

Now we have the theorem that W is greatest for the uniform distributio: 
of the balls in the baskets, and the assertion that the most probable state 
of a gas is the state of uniform density, all ready to be fitted together. The 
process of fitting-together is of the simplest. W is christened the “‘prob- 
ability” of the state described by the “distribution” N,, Ne, «++ Nw, the 
quantities V,; now standing for the numbers of molecules in the various cells. 
Not only is the state of uniformity the most probable one by this definition, 


, 








fact 
prot 
was 
not 
inkl; 
e 
take 
dens 
infin 
OC 
at tl 
corre 
him 


sphe 








MEMORIAL TO CLASSICAL 





STATISTICS 117 


t so long as the number of molecules N is many times as great as the number of 
ompartments M—a condition easy to realize—those distributions which are 

rkedly far from uniform have probabilities which are fantastically smaller 
than the value of W for the uniform state. 

[he Boltzmann statistics manages thus to derive the assertion aforesaid 
he assertion that the uniform distribution in ordinary space is of all the 
most probable—from a principle which (at least in appearance) is more fun- 
damental. It has indeed a couple of bothersome points—more than a couple 
perhaps, but there are two in particular which the newer statistics will at- 
tempt to assuage. One of these is the size to be assigned to the cells Vo ; but 
we are borrowing trouble to think too much of that now, since whatever 
choice be made so long as N/M be large will not affect the achievement just 
cited. The other is, that one would much rather think of the molecules of a 
gas (of a single chemical kind) as being alike absolutely, than as being dis- 
tinguished one from another by a mysterious something-or-other represented 
in this theory by numbers painted on balls. In the Boltzmann statistics, 
however, the numbers must stay on the balls. 

We go over into the momentum-space, setting up a coordinate-frame and 
representing the molecules by dots, the coordinates of which are the momen- 
tum-components p,, py, p: of the molecules in question. To each position of 
a dot corresponds an energy-value, equal to (1/2m) (p: + p, + p:); we will 
callit £. E vanishes at the origin, and has a constant value over any spheri- 
cal surface centered at the origin. To any distribution of the dots will cor- 
respond a specific value for the total energy of the gas. For this we need a 
symbol different from £; and as we shall have a good deal to do with thermo- 
dynamics later on, I choose the thermodynamical symbol U. The average 
energy of the molecules of the gas will then be U’/N, to be denoted by U. 

lhe entry of E and U into the situation is of the first importance. It is in 
fact all that will save us from the highly unwanted conclusion that the most 
probable distribution in the momentum-space is the uniform one, just as it 
was in the coordinate-space. To see why it makes so great a difference is 
not altogether easy. I think that the reflections which follow may give an 
inkling of the reason. 

The momentum-space must be taken either as infinite or as finite. If we 
take it as infinite and demand a distribution of uniform density, then the 
density goes to zero and at the same time the energies of the molecules go to 
infinity, producing an impossible situation. Let us then take it as finite, 
blocking off all of the parts of it which lie beyond a certain sphere centered 
at the origin. Assume a uniform distribution within the sphere. This will 
correspond to a certain value of U. (The student may suppose, if it makes 
him happier, that the U-value was preassigned and the radius of the ‘certain 
sphere” chosen accordingly.) The W-value of this distribution will surely 
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be greater than that for any non-uniform distribution, whether of the sam: 
U’ or of a different U’, confined within the sphere. However, by blockin 

the whole of the momentum-space beyond the sphere, we have barr: 
whole lot of distributions corresponding to the same U and having so1 
their dots beyond the sphere. By no means have we proved that the || 
value for the uniform distribution within the sphere is greater than that for 
any and all of the barred distributions. Now if we can agree that the b|: 
ing-off of part of the momentum-space is a silly thing to do and unacceptahile 
to Nature, the argument for the uniform distribution is spoiled, and we hav 
to look for a new idea. 

At this point it seems best to go through the mathematical process | 
finding the distribution of greatest W in the coordinate-space and 
momentum-space, just as that process is presented in the textbooks. 

We return to equation (1) and make it a manageable one by having re- 
course to that godsend of statistical mechanics, the ‘Stirling approxima 


tion”, which may be written thus: 
InN! = NInN — N+ In V24N 2 


This is valid only for large values of VV, though writers on S.M. never seem 
to remember how large the values must be. For still larger values of .V we 
can drop off the last two terms, arriving at a sort of super-Stirling approxima 
tion which however itself is commonly called the Stirling approximation: 


InN! = NInN 3 
Putting (3) into (1), we find: 
nW=NinN—-ZN;InN; } 
Defining some quantities w, by the equations V; = Nw; , we make this over 
into: 
nW = —N Zw, In», 5 


having availed ourselves of the obvious fact that Sw; is equal to unity. 

We might now convert this into an equation for W, but this would bea 
waste of time and energy, since whenever W has a maximum so also will 
InW. With In W, therefore, we operate from now on. Making small varia- 
tions in the quantities V; , and making therefore small variations—call them 
éw;—in the quantities w;, we find in first approximation for the ensuing 


change in In W, 


dlnW = —NI(1+ Inw) dw, 6) 
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Now we are restricting ourselves to variations in the quantities V; which 
leave unchanged the total number of molecules in the cells, or of balls in the 
haskets—to variations, therefore, for which 


> N; = N = constant, ~b6w; = 0 7) 


This restriction being introduced into (6), 6 In W proceeds to vanish if and only 
if w; has the same value for all of the cells. Now, the vanishment of 6 /n W is 
a necessary condition for having a maximum of W’ at the situation in ques- 
tion. I do not refer to it as a sufficient condition, because it admits ofa 
minimum or of what is technically known as a “stationary” value of W in the 
situation in question. However it has already been shown, without the aid 
of the Stirling approximation, that the expression to which we are approxi- 
nating is greater for the uniform distribution than for the neighboring non- 
uniform ones. It may therefore be accepted that here we have a maximum 
of W for the uniform distribution, and have reached the old result in a new 
way; an achievement nearly useless, were it not a prelude to the performance 
in momentum-space. 

I continue to use the symbols W and NV and JV; and w; , but now with ref- 
erence to the distribution of the representative dots in momentum-space. 
\ new symbol, £; , shall signify the energy of a molecule in the ith cell of the 
momentum-space. We wish at all costs to avoid the conclusion that the 
stable distribution in the momentum-space is the uniform one. Boltzmann 
managed to avoid it, and his was the following way: 

Let us write, for the number of molecules in the 7th cell, the expression: 


Nw; = NA exp (—BE;) (8) 
and insert it into (6). We shall find: 
6lInW = — NY(1 + InA — BE,) Gw;. (9) 


Of the three terms on the right, two vanish for all variations in which the 
total number of molecules remains the same. The third does not—but it 
will vanish for a restricted class of these variations, to wit, those and those 
only for which the total energy of all the molecules remains the same; for 
\xw £; is precisely that total energy. 

Some writers at this point ask the student to imagine a. gas in acontainer 
being completely cut off from energy-interchange with the container-walls 
and with the whole of the outside world, and therefore being limited to the 
particular U-value with which it started out. Others import the word ‘‘tem- 
perature” which I am desperately (and vainly) trying to keep out until I am 
ready to bring it formally into the discourse, and aver that the gas is nearly 
or quite so limited if the walls of the container have the same temperature as 
the gas itself.. The student may take his choice, but must suppose that 
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under such conditions Nature rejects that distribution which so to sp: 
“stable against” every conceivable variation, and elects that peculiar 
bution which is stable not against any conceivable variation but only ag 
the possible ones. Perhaps this is because the uniform distribution wo 
entail the consequences mentioned on page 117, or perhaps there is no 
in saying that it is ‘‘because”’ of this or “because” of that. Anyhoy 
peculiar distribution is the one which the data sustain. 

However I have not really defined the peculiar distribution as yet, h 
merely thrown the symbols 4 and B into equation (8) as though they 
for completely disposable constants. It can readily be seen that a 
most there can be but one disposable constant, for A and B interlinke 
the obvious equation: 


= w= AL exp (—BE;) = 1 1 





But even B is not disposable, if the total energy U and the average energy 


per molecule U are preassigned; for there is another obvious equation: 


U => Ew; = A Exp (—BE,) 11 


What with equations (10) and (11), there is no longer anything disposab) 


about the constants A and B. The peculiar distribution in the momentum- 


space is completely defined. It is the Maxwell-Boltzmann distribution-| 


obtained from the Maxwell statistics, and sometimes known as the ‘‘cano1 


cal’”’ distribution. 


To summarize now the Boltzmann statistics as on page 113 the Maxw: 


statistics was summarized: 
A gas is more likely to be found in its most probable state than in any ot! 


é? 


The probability of a state is found by imagining it as a distribution of numbered 
molecules among cells, in the coordinate-space and in the momentum-s pac: 


That of any distribution is measured by the number of inventories com patil 


therewith. By this criterion the most probable distribution in coordinate-s pac 
is the uniform one, and by this criterion carefully hedged about, the most probabil 


distribution in momentum-space is the Maxwell-Boltzmann or canonical 
It is necessary to liken molecules of a single kind to numbered balls, differin 
no way except the numbering. 


This point was reached by statistical mechanics about fifty years ay 


Had it not been for Planck’s wish and tenacious will to explain the bla 
body. radiation-law, it might have been the stopping-point. 


A HELPFUL AND TROUBLESOME COINCIDENCE BETWEEN Two 
DIFFERENT QUANTITIES 


Let us return to the game with the sack, the balls and the basket, pla; 


in the manner which led to good results when applied to the molecules in | 


coordinate-space. 


fi 
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[he most probable distribution is the one evoked-by the caller, when he 

s for an equal number of balls in every basket. If there are N balls and 

V. baskets, this means V/M balls to each basket, and a maximum number of 

inventories which I will call Wy... Looking back to equation (1), we see 

Wnax iS a fraction the numerator of which is \-factorial, while the de- 

nominator is (V/M)-factorial raised to the power M@. Taking logarithms 
| using the super-Stirling approximation, we find: 


ng In Wax = Nin M (12) 


[he logarithm of the probability of the most probable distribution (of num- 
bered balls in numbered baskets, or molecules in equal cells of coordinate- 
space) is equal to the logarithm of the number of baskets (cells), multiplied 
by the number of balls (molecules). 

Next suppose the caller, in a fit of uncontrollable zest for the game, calling 
rg in succession every one of the conceivable distributions. What is the total 
number of inventories compatible with all of them together? To sum over 
every conceivable expression of the type of (1) seems a hopeless assignment, 
but there is a short-cut to the result. 
abl Fix a particular order for the drawing of the balls from the sacks—it may 
— as well be the very order of their numbering. The first of the balls to be 
drawn may be tossed into any one of the baskets, giving M distinct “possibil- 
1oni- ities”. The second may be tossed into any one of the baskets, the same or 
another, giving in conjunction with the fate of the first M@° different possibil- 
we ities. The third may be tossed—but we leap to the conclusion. ‘There 
ire M¥% possibilities altogether, and these are the inventories. Thus the total 
thes number of inventories consistent with all of the distributions, which I will 
hered call Woe, is a number whereof the logarithm is, 


In Wee = Vin M (13) 


pace But this is the same as the expression for Jn Winax! 
hat The meaning of this strange coincidence can only be, that when N and 
1 \/M are both so great that the super-Stirling approximation is a good one, 
40 41 then the logarithm of the number of inventories belonging to the most prob- 
ible distribution is nearly as great as the /ogarithm of the total number of 
inventories belonging to all of the distributions put together—so nearly as 
lack- great, that either Jogarithm is a good approximation to the other. 

In the foregoing very important paragraph, I have italicized the word 

garithm” because if it were left out the statement would become a false 
ne. The statement is not true if applied to the numbers themselves. W tor 
is manyfold greater than Wyax, and the ratio between the two actually in- 
creases with rising V. So does the difference between In Wo and In Wiyax 
rease with rising V, but not so fast as either by itself; wherefore the truth 
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of the statement. The student may convince himself of this by app 


the second-degree Stirling approximation (equation 2)?. 

I have called this both a helpful and a troublesome coincidence. It 
be deemed a helpful one, because the expression for the total number « 
ventories is easier to derive and easier to remember than the expression {or 
the number of inventories belonging to the most probable distribution. [i 
therefore one has good ground for believing (as here is the case) that th 


logarithms of the two are approximately equal, one may serenely remem)er 
and use In Wor instead of In Winax . The troublesome feature is, that some 
expositors speak of In Wor throughout and never allude to In Wysx , thus 
confusing the student to an extent which (if my experience is typical) may 


well be serious. I shall later dwell on the fact that In W for any distribution 
is regarded as a measure of the entropy of that distribution, and In |i 
therefore as a measure of the entropy of the most probable distributio: 
Some people imply that In W’,,, is the true measure of the entropy of the gas, 
instead of being an approximation to it. They commit no numerical error 
in so doing, but they blot out the most remarkable quality of the Boltzman: 
statistics, to wit, the clear distinction which it makes between the most prob- 
able distribution and those of lesser probability. This mistake is mor 
commonly made in treating the newer statistics. Here I am not so sure 
that it is a mistake, but I think so. 


MEANINGS OF THE WorD “STATE” 


The word ‘‘state”’, which turns up continually in this essay, is one of those 
words of which a proper definition is hardly less than a full description of th« 
theory which employs it. When the theory changes so also does the meaning 
of the word. In the welter of statistical theories, the word “‘state’’ has 
several different meanings. In thermodynamics also it has more than one 
meaning, but one is preeminent. 

Thermodynamics usually concerns itself with gases (not to speak of |i- 
quids and solids) which are in what I earlier called a “uniform” state: uniform 
density, uniform pressure, uniform temperature. For a gas of a single kind 
(‘‘kind”’ being a word which it is the business of chemistry to define) it is a 
fact of experience that any two of these three variables suffice to define the 
third and also all of the other variables which thermodynamics cares about 
Of these others there are two in particular which I mention at this point, the 


2 Actually if one goes from the “most probable state’ V; = const. = N/M to t 
“next most probable” in which one ball is taken out of one of the baskets and put into 
another, the change in W is in the ratio of (N/M) to (N/M) + 1, which is practical! 
no change at all when V/M is so high as is commonly taken. This shows that the stat: 
ment could not be true if it were made about the numbers Wax and Wot rather that 
about the logarithms thereof. It certainly looks as though the statement could not bi 
true even when made of the logarithms, but this is evidently one of the cases wher 
“intuition” is a fallible guide. 
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energy and the entropy. This makes five altogether, and any two of the 
five suffice to determine the state-—THE STATE, the uniform state, the 
only one about which thermodynamics really knows or cares. When asked 
about what I earlier called ‘“‘a surging state’, thermodynamics mutters 
something to the effect that the entropy of such a state is smaller than that 
of THE STATE, and then puts an end to the conversation by refusing to 
commit itself further. Thermodynamics takes no cognizance of the molecu- 
lar structure of matter. A gas might be a continuum, for all that it knows 
or cares. 

Statistical mechanics talks about a mental image of the gas, in the form of 
a flock of dots in the coordinate-space and another flock of dots in the mo- 
mentum-space, or one may call them a single flock of dots in the u-space. In 
Boltzmann statistics, the ‘‘state”’ of this image is what I have been calling 
the ‘distribution’. The most probable state of the image—to wit, the one 
with the greatest number of inventories or complexions—is identified with 
THE STATE of thermodynamics. All of the rest belong to the category of 
which thermodynamics would say, that the entropy is smaller than it is for 
THE STATE. But since according to S.M. they belong to a category for 
which the probability is smaller than it is for THE STATE, one sees a con- 
nection between entropy of the gas and probability of the image beginning 
to take shape. 

Now it is time to make a formal introduction of the concepts of entropy 
and temperature—the latter word having already sneaked into this article 
two or three times in spite of all my efforts to keep it out. 


FORMAL ENTRANCE OF ENTROPY AND TEMPERATURE 
For a substance, meaning now a gas, of a single kind, entropy and tempera- 
ture are defined by the equation, 


dU = TdS — PdV (14) 


P stands for pressure, V for volume, and S for entropy. For energy I use 
the symbol U’ already employed in that sense—but notice that formerly it 
stood for the kinetic energy of the molecules! To use the same symbol in 
both senses implies that the energy of the gas is entirely the kinetic energy 
of translatory motion) of the molecules. This identification turns out to be 
valid for the ‘‘monatomic” gases, which are luckily numerous and well- 
studied. To these we confine ourselves throughout this article. 7 stands 
for the temperature called absolute; this being the only kind of temperature 
which will ever figure in this article, the adjective henceforth is discarded. 
Density was the fifth variable in my list given above, but volume is usually 
preferred to it. To make them equally useful, the quantity of gas must be 
stated; here it will be taken as one gramme-molecule. 
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It is evident that the equation is a comparison between two states 


not go astray by supposing that these are like two of the states whi 
have been considering, having the same l’ and V and differing in the nu 
of inventories! These on the contrary are two examples of THE STA’! 
of the thermodynamic state, of the most probable state—of a gas, difi 
in the values of some at least among the five variables. The quantitic 
and dS and dV are the differences between the U-values and the.S, alues 
the V-values of the two states, while P and 7 may be taken as reterri 
either, the smallness of the difference between the two states—implie: 
the differential notation—permitting of this. 

It is also evident from my wording that the one equation is being | 
to define the two quantities S and 7. This is unluckily no verbal 
nor is it a temporary shortcut to be replaced by a royal road as the 
gument proceeds. The meanings of entropy and temperature are 
coiled up together in thermodynamics, that it is impossible to take t 
apart unmutilated. One cannot seize either by storm and then invest | 
other, at least not without the aid of statistical theory: one has to surroi 
them both ina single campaign. As F, jington has vigorously written, this 
isa common thing in physics. Elect iorce is defined as that which acts 0; 
electric charge, electric charge as tl vhich is acted upon by electric force 


and so on.... Common as it mags.be, it is probably nowhere else so har 
assing as in thermodynamics. Thre are three ways of intruding upon th 
vicious circle. 

First, to apprehend both concepts in a single mental act. This is th 
counsel of perfection. 

Second, to use a temporary definition of temperature, with the promise oi 
confirming or correcting it later. The ideal-gas thermometer is the devic 
used for this purpose in thermodynamics. Anyone trained in this way is 
likely to think for the rest of his life of temperature as the primary concept, 
entropy as a derived one—as indeed was the case, when thermodynamics 
started. 

Third, to produce a theory which makes a pronouncement as to the nature 
of entropy. ; 

This last is the major office of statistical mechanics. To those who accept 
it, entropy becomes the primary concept and temperature the derived one, 
and both are visualized by the aid of the key-word “‘probability” of the basi 
theorem, interpreted in some particular way. 


OLD STATISTICAL THEORY OF ENTROPY 


In the classical statistics, the entropy of a distribution is considered to be t/ 
logarithm of the number of inventories or complexions compatible with that dis- 
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tribution, multiplied by a constant (always denoted by k) which is adjusted to 


bring about agreement with experiment: 


entropy S = kin W. (15) 


To illustrate this doctrine and to evaluate k, I now take the student back 
to the coor linate-space, where a box of volume V populated with V mole- 
cules is “ led mentally into M equal cells of volume V,,and the most 
probable distribution is characterized by the value \V /n M for the logarithm 
of the probability. The entropy—or no, not the entire entropy of the gas, 
but merely what I will call “the contribution of the volume of the entropy” 
and denote by S,—is then supposed to be kN Jn M, or: 


Se = kN In V — kN In Vo (16) 


Reverting to the equation (14) in which the definitions of entropy and 
temperature were tangled up together, and rearranging it, we get: 


TdS = dU + PdV (17) 


Now, an “ideal gas” is defined by tw. ttributes. First, there exists between 
its pressure and its volume and its perature the relation P = aT/V, 
wherein a stands for a constant. Sec , its energy U depends upon the 
temperature only, and not upon any ot. - variable, in particular not upon 
the volume. Therefore we may write: 


TdS = C,dT + (aT/V)dV (18) 


C, here standing for something of which we need only know that it is a 
function of T alone. Integrating, we find: 


S = R In V + (function of temperature) + constant (19) 


and lo! it is seen that the dependence of entropy on volume is precisely of 
the sort which the theory is fitted to explain. 

The next step is to adjust the value of the constant k. The constant a 
aforesaid is proportional to the amount of gas in the box, proportional there- 
fore to NV: it is the constant ratio of a to V to which k must be equated. For 
the amount of gas let us choose one gramme-molecule. Then a assumes the 
value always symbolized by R and called the ‘‘gas-constant’’, and V assumes 
the value usually symbolized by No and called the ““Avogadro number’’. 
Both of these are known from experiment, and & is fixed by the equation 


k = R/No (20) 
The constant & is named in Boltzmann’s honor, though in his time its value 


was not known because the value of No was only vaguely apprehended. 
Now we have settled what I called “‘the contribution of volume to en- 
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tropy”’. It remains to interpret the rest of the right-hand member 0! (19 Com] 
which I will call ‘‘the contribution of temperature to entropy”. To do this state 
we must re-enter the momentum-space. hed \ 
From (15) and (5) and (8) we get, for the entropy S,, of the flock cons! 
in the momentum-space: Uji 
proc 
Sn = —kInW = — kN zu; Inw; ae 
dS. 
, : . BE : 
= —kNAS(InA — BE,e °"' 2 I 
Refreshing our memory from (10), we see that the first term of this « 
sion reduces to —kN In A. Refreshing our memory from (11), we see tha: rhis 
the second term reduces to + RNBU or RBU. Referring now ols 
aiu 
gramme-molecule of gas, I put R for Vk, and find: lie 
Sm = —RInA + kBU ) tinu 
= ose ; ae H 
Sm is hereby given as a function of U’, but a more complicated functio: 
appears on the surface, since A depends upon B (equation 10) and B upoi 
(equation 11). Yet when we differentiate S,, with respect to U’, and 
doing take account of these complications, it turns out that we might as w: the 
have been oblivious of them! for the result is the same as though A a Thi 
Al 
were constants: divi 
ii 
dS,/dU = kB 23 con 
Now the temperature, which has so often slipped into this argument 
ways more or less surreptitious, is about to make its formal and ceremo: ai 
. . . . . . - 0 
entry into the statistical picture. We turn back to equation (17), r 
- 101 
deduce: : 
qul 
dS/dU = 1/T 24 the 
. : : : . . one 
rhe derivative here standing on the left is the derivative of entropy wit 1 
ae ; 
respect to energy under the condition of constant volume: a thermodynan 
‘ i. it sme knee . P P ee. Va 
cist would write it (0S/dU’), . It is therefore properly to be identified wit ) 
ss sen ; ik ‘ ap 
the derivative in (23), and we make the two identical by putting: o 
i 
B = 1/kT 25 va 
. : ‘ : ‘ . ter 
Now taking the entropy S to be the sum of S, and S,, , we find: 
S=S.+ Sn = —RINA+ U/T+ RIinV — Rin Vo 2 
and this is to be compared with (19), the thermodynamic expression for sO 
entropy, which I repeat to make the comparison easier: 
S = | (C,/T)dT + Rin V + constant 
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Comparing these, we see first of all that R/n V appears in both, as was already 
stated. It also seems at first glance that (—R/n A + U'/T) is to be identi- 
fied with the integral in (27), and that —R/n V9 is to be identified with the 
constant in (27). This however is not necessarily the case, for (—Rin A 4+ 
|'/ 7) may prove to include constant terms. Indeed they do; and we must 
proceed to evaluate both A and U in terms of 7 in order to round off the 
task. 


I recall equation (10) and write it thus: 
1/A = 2 exp (—E;/kT) 28) 


[his is a summation, to which each cell contributes one term having the 
value of E appropriate to that cell—£; for the ith cell. Of the volumes of 
these cells I have thus far said nothing, except that all are equal. I con- 
tinue to say nothing further, but I give to their common volume the symbol 
IT). Let us now form the integral: 


, Pf 


[| exp (—E/kT) dp.dp,dp., E = (1/2m) (p: + p, + p:) (29 


the range of integration extending over the whole of momentum-space. 
This integral may be described as follows. Let the momentum-space be 
divided into cells of unit volume. Each of these cells of unit volume makes a 
contribution 


exp (—E/kT) 


to the integral, & standing now for the average value of E in the cell in ques- 
tion. The integral is the sum of all of these contributions. Now let us in- 
quire how much of a contribution is made by this same cell of unit volume to 
the summation (28). This second contribution is made up of 1/H terms, 
one for each of the cells of volume Hy which occupy the cell of unit volume. 
The values E; corresponding to these cells will not be exactly equal to the 
value E corresponding to the entire cell of unit volume; but to the degree of 
approximation which is now being used, the difference may be neglected. 
The summation (28) is then equal to 1/H times the integral (25). Now the 
value of the integral (29) is given in all tables of definite integrals, and in 
terms of our symbols it amounts to 


(2amkT)*” 
so we come to the conclusion: 
In A = — In (2rmkT)*” + In Ho 


> 


= —3InT — In (24mk)*” + In Ho (30) 
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Now we have attended to every term in (26) except the term U/T. Near 
every reader will remember that the average kinetic energy of an atom of. 
monatomic gas at temperature T is kT. I therefore leave out the ci 
tion of this result, except for showing the student how to begin on it 
first step is to go back to equation (11) where an expression was given f{ 
and in that expression to replace the summation 2 E; exp (—BE 


PrP fp Pp 


(1/Ho) times the integral | | E exp (—BE)dp, dp, dp,. It follows : 


U/T is (3/2) Nk, which for one gramme-molecule of gas is (3/2)R, which | 
write as R In e°* 

The picture of entropy for a monatomic gas limned by the Boltzmann sta- 
tistics, is now completed. Entropy is the function which follows: 


3 (2xmke)"” 
S=-RinT + Rin V + Rin ——— 
2 Vo Ho 
The dependence on volume is correct, 7.e., just the same as in the thermody- 
namic formula. The dependence on temperature is correct, for (3/2)R is 
the value of the specific heat at constant volume per gramme-molecule of a 
gas, the quantity C, of equation (18). The additive constant, as to the 
value of which thermodynamics says nothing, is fixed when the volumes 
Vand Ho of the elementary cells in the ordinary space and the momentum- 
space are fixed. 
MIXTURES OF GASES 
Now we will go through the mental operation which is called ‘considering 
a mixture” of two different monatomic gases, V’ atoms of the one and .\ 
atoms of the other, in the same box and (necessarily) in the same momentun 
space. Let me denote by U’ and U”’, respectively, the energies of thes 
two gases; and by V; = N’w; and NV, = Nw;, respectively, the numbers of 


’ 


atoms of the two kinds in the ith cell of momentum-space. 
If we seek the most probable distribution of the first gas in the momentum 

space, making the stipulation that we will admit only such variations of the 
quantities w; as leave V’ and U’ unchanged—well, of course, we get the sam 
result as before, the distribution (8), with N’ in place of N and (let me say 

A’ in place of A and B’ in place of B. A’ will depend upon B’ and B’ will 
depend upon U’/N’. If we do the like with the second gas, we get anew to 
the distribution (8) with V’’, A” and B” in place of NV’, A’ and B’. A” will 
not be the same as 4’ nor will B” be the same as B’, unless it happens that 
U”/N" is equal to U’/N’. There is no cause for surprise in this. In acting 
this way we are only treating each gas by itself, and have as yet done nothing 
which can be regarded as ‘‘considering a mixture”’. 

Let us however seek the most probable distribution of the two gases, mak- 
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ing the stipulation that we will admit only such variations of the quantities 
and w; aS leave N’ and N” and the sum of the energies U’ and U’"—not how- 
ever the individual energies U’ and U’’—unchanged. In acting this way 
we are doing something which may be regarded as “‘considering a mixture’”’, 
since we are allowing for the possibility that energy may pass from the one 
gas to the other and the other to the one. Equally well are we considering 
the case of two gases separated by a partition through which energy may 
pass, but not the atoms. Since in such a case we really ought to take into 
account the atoms and the energy of the partition also, we must appease the 
critics by providing that the partition shall be very thin. 

Choose any set of values of the quantities N;, which is to say, any particu- 
lar distribution of the first gas; and choose any set of values of the quantities 
Ni 
to equation (1) and put primes on all the symbols V, ¥;, V2, --- on the 


’; , which is to say, any particular distribution of the second gas. Go back 
right-hand side of that equation. The resulting expression gives the total 
number of inventories or complexions of the first gas. Take off the primes 
and affix double primes to each of these symbols. The resulting expression 
gives the total number of inventories or complexions of the second gas. 
Every complexion of either may coexist with any complexion of the other. 
Therefore the total number of complexions of the pair of gases is the product 
of the two expressions. It is this product which is W for the pair of gases, be 
they mixed or side-by-side. 

With use of the Stirling approximation, the logarithm of W for the pair 
is the sum of two such expressions as we have seen in (5): 


vl 


, tin ? , ris ” 2 
InW = —N’Sw; In w; — NS; In w; 32) 
and its variation is: 
dln W = —N’Z(1 + In w,)dw; — N’Z(1 + In w” bw, 33) 
Let us now give a trial to the tentative distribution, 
, , open 5 / > , 2 
w; = A’ exp (—B’E,), w’’ = A” exp (—B” E,) (34) 


On substituting this into (33) we find that if B’ is unequal to B’’, the dis- 
tribution has a stationary value of W with respect only to such variations 
as leave the energies of the two gases separately unchanged—the result 
which we had before. If however B’ and B” are the same, then W is 
stationary with respect to variations which leave the sum of the energies 
unchanged, either being allowed to gain or lose so long as the other loses or 
gains by an equal amount. Since each B is controlled by the corresponding 
U/N, the distribution (33) has a stationary value of W for variations of the 
type in question if and only if the average energy of the atoms of each gas 
is the same. Since each B controls the corresponding A, this condition of 
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equal average energy makes the distributions of the two gases jus 
same. 

We have already seen that &B is the reciprocal of the temperature 
it is the reciprocal of (0U’/AS), in our statistical picture, and the defir 
of absolute temperature 7 is precisely that 7 is this derivative. The si 
ment to which we have come is, that the most probable state of the mixtw 
the one in which T is the same for both components. It is often expresse: 
this way: classical statistics shows that for two (or more) gases in equilibri 
with each other, the temperature must be the same. It is indeed a fact 
experience, and a most important one, that when two systems (be they gases 
or be they not) are in thermal equilibrium, their temperatures are the sam: 
This has not hitherto been mentioned, and yet we seem to have derived 
(Juite a rabbit for the magician of the classical statistics to have pulled | 
of the hat! 

However, skeptical people who see a rabbit pulled out of a hat are inclined 
suspect that either the rabbit was in the hat beforehand, or else there is 1 
rabbit. Let us inquire into the contents of the hat and see whether we cai 
find the rabbit there. 

The first (and the last) question to be asked is: what is the difference } 
tween ‘‘different” kinds of gas in the statistical picture? 

To the physicist or the chemist, different kinds of gas will be (for exampli 
mercury and helium. These differ in their spectra, boiling-points, chemica! 
properties, and quantities of other features. None of these features however 
appears in the theory, and therefore none of them can contribute to th 
result. The atoms also differ in mass, and for a moment this seems to be 
difference of which the statistical picture takes account, since the letter ™ 
appears in some of our equations. However, it appears only in the ultimate 
equations, those such as (29) in which the distribution-in-momentum is 
expressed. It does not appear in the original form of the Maxwell-Bolt: 
mann distribution-in-energy, the form shown in equation (8). It appears in 
particular in the last term of equation (31), but not elsewhere. Apart fron 
this it may be said that in the classical statistics, all gases are the same gas 

This is a paradox, but only one of two. The other paradox is, that in 
the classical statistics two parts of the same gas are different gases. This second 
paradox arises from the numbering of the molecules, which is an essentia 
feature of the classical statistics. 

Therefore in the statistical picture a mixture of NV’ atoms of mercury and 
N” atoms of helium is distinguished by the fact that the mercury atoms bear 
one set of integer numbers (say those from 1 to NV’) and the helium atoms 
another set (say those from V’ + 1 to VN’ + N”). But if the atoms were al! 
helium atoms or all mercury atoms, they would also be divisible in man) 
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lifferent ways into a set of NV’ atoms bearing one set of numbers and a set 
of NV” atoms bearing another set of numbers. Each set would obviously 
have to have the same distribution, with the same 1 and the same BA, as any 
other set or as the totality of all the atoms. This conclusion, which is self- 
evident in the case in which all the atoms are called ‘‘mercury”’, remains true 
when some of the atoms are called ‘‘mercury” and others are called ‘thelium”’. 
We have done nothing but change the names of some of the atoms: we have 
not imported into our theory anything which differentiates one kind of 
atom from another kind. No wonder we have arrived at the conclusion that 
all kinds have the same distribution-in-energy, the same 1, the same B and 
the same temperature! The rabbit was indeed in the hat, but it does not 
look like so much of a rabbit. 

The classical statistics therefore doesn’t recognize any of the real dif- 
ferences between atoms of different kinds, except for alterations in the last 
term of (31); but it does make an artificial difference which creates the 
astonishing result, that any two samples of the same gas are different gases! 
At this point we may begin to wonder whether this peculiarity, which has 
led to so apparently brilliant a result in respect of the equality of tempera- 
tures in thermal equilibrium, might elsewhere lead us astray. It does; and 
here appears the rift in the lute of classical statistics. 


THE RIFT IN THE LUTE 


Let us imagine two boxes of equal size separated by a common partition, 
each containing a gas consisting of V atoms, both gases at the same tempera- 
ture. We will baptize one gas “mercury” and the other gas ‘“helium”’. 
Let an opening be made through the partition. It is known that in such a 
situation in Nature, the two gases diffuse into one another, the final and 
permanent condition being that in which the mercury and the helium are 
equally distributed between the two boxes. The process of diffusion is an 
example of what in thermodynamics is called an “irreversible” process. The 
state of uniform mixing ought to correspond to the most probable state in 
the statistical picture. But what does the statistical theory say? 

The statistical theory says nothing about diffusion and nothing about 
mixing. ‘The statistical theory takes account of nothing but the facts that 
the mercury had at its disposal the volume V before and the volume 2V 
after the breaking of the partition, and ditto for the helium. The value V 
contains M cells (M = V/Vo) and the volume 2V contains 2M cells. The 
(approximate) probabilities of the uniform distribution are M* before and 
(2M) after. The latter is greater than the former; the entropy goes up by 
Nk \n 2 for each gas, by 2 Nk In 2 for the two of them, when the private pre- 
serve of each is thrown open to the other. This gain is what is called the 
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b] 


though as we have seen it is really the “entro; 


“entropy of mixing’ 
expansion”. It is the alteration in the second term of the righthand mer 
of (31). 

But now suppose both of the boxes hold helium. One may indeed 
tinue to suppose that when the partition is opened each one of the 
samples of helium undergoes an expansion, doubling its volume. 
entropy would then go up by 2 VkIn 2. However this looks so silly a t! 
to say that no one, I feel almost secure in affirming, has ever said it. 
natural thing to say is, that the 2V atoms of helium distributed thro 
the two boxes at uniform temperature and uniform pressure have just | 
same entropy-value whether or not the partition is broken. 

What does the classical statistics say about this situation? Its answer ca: 
be foretold. Since the two samples of helium are different by virtue of the 
different numberings of the two sets of atoms, the classical statistics insists 
that the entropy increase by 2NVk/n 2 when the partition is broken, even 
though the gases are the same. This is indeed, if I may pervert the poem, 
“the little rift within the lute, which makes the classical statistics mute.” 
The achievement of predicting the uniform distribution in ordinary space, 
the achievement of predicting the Maxwell-Boltzmann distribution-law in 
momentum-space, the achievement of providing the proper relation between 
temperature and mean kinetic energy—all of these are unsettled by this 
calamity. 

Were I writing a strictly logical article I should quit at this point. No- 
thing further can apparently be done, except to tamper with the classical! 
statistics in an effort to remove the unwanted result which has sprung forth 
to plague us. To violate the logic of the classical statistics in order to 
banish the undesired while keeping the desired results is a very questionable 
act. In theoretical physics, it is not admissible that the end justifies an) 
and all means. Nevertheless so successful a feat of tampering has been 
done, that I cannot refrain from mentioning it as I close. 

Let me first express in a slightly different way the nature of the “rift”’. 
Compare two samples of the same gas at the same temperature, one con- 
sisting of N atoms in a volume V, the other consisting of xN atoms in a 
volume xV. That which is called entropy in thermodynamics—and there- 
fore that which is entropy, since it is the privilege of thermodynamics to 
give the definition of entropy—is x times as great for the latter as for the 
former. But that which the classical statistics calls entropy—or, as we must 
admit, miscalls entropy—is not » times as great for the latter as for the 
former. It would be, if there were x times as many atoms but just the same 
number of cells. However, there are « times as many atoms but also 
times as many cells into which to put them. The number of complexions is 
approximately M* in the former case and («M)*" in the latter, M standing 
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the number of cells in the former box. The thing miscalled entropy is 
iM in the first case and (kxNInM + kxN]nx) in the second case. It is 
term kxNInx which is the rift. 

Clearly we could abolish this term by allowing the volume of the cells to 
swell in the ratio +:1 when going from the former case to the latter. This 
; the same as making Hp proportional to the number of atoms in the sample 

f gas which happens to be under study. Since in equation (31) the volumes 
1, and Ho (of the elementary cells in ordinary space and in momentum- 
space) are indissolubly bound together in the product Vof/» , this is the same 
as making Vof79 equal to some constant multiplied by the number of atoms 
under study. 

Such, if I interpret correctly, was the idea proposed by Sackur in 1912. 
While it does the task required, it is an “‘ad hoc” assumption of the most 
barefaced character. If the gas under study is at first divided into two 
parts by a partition and the partition is then abolished, the cells must be 
supposed to swell up at the moment when the partition vanishes. 

We can also abolish the fatal term by going back to equation (1) for the 
number of complexions, and removing the factor V! in the numerator and 
replacing it by unity. We then have unity divided by the original de- 
nominator, which in the (most probable) case of the uniform distribution is 
N/M)! raised to the power M, as I remarked on page 121. Using the 
super-Stirling approximation, we find that the logarithm of one fraction is 
NiInM — NinN). The factor NV! which we formerly had in the numerator 
killed off the term (— NIn.V), but now that we have taken it out, this term 
survives. If now we say that & times the logarithm of W/N! shall be the 
picture of entropy in the classical statistics, then the term (—kNInJN) 
comes over into the right-hand member of (31). It may be amalgamated 
with the last term already standing there; and when this is done, we find 
Vol7y multiplied by N exactly as Sackur put it there, and with the same 
wished-for result. 

This, if I interpret correctly, is the idea proposed in 1913 by Tetrode. It 
does the task required of it, but its drawback is that the removal of the 
factor N! from the right-hand member of (1), a drastic piece of surgery as 
it were, violates the system of the classical statistics.’ 

I was not, however, thinking merely of this achievement when on Page 
132 I spoke of ‘‘a remarkably successful feat of tampering.” To show the 

’ This may seem too strong a statement. We are, after all, only asked to accept & | 
W/N!) as our picture of entropy, instead of kInW; why be reluctant? But in effect, as I 
see it, we are asked first to accept klnWf as our picture of entropy, f being an arbitrary 
function of N; and then we are asked so to choose f, that the dependence of k In Wf on A 
shall conform to the actual behavior of entropy. This is different from and much less 

npressive than our original procedure, which consisted in first realizing that W is the 


number of complexions, and then discovering that k In W depends on volume and on 
temperature in just the right ways for entropy. 
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magnitude of the achievement, I will rewrite equation (31) with two 
tions. The first consists in replacing R with NV, so that the expressio 

refer not to a gramme-molecule of gas but to any number N of atoms. 1 
second consists in following Tetrode by affirming that the entropy 

k ln W, but & times the logarithm of W/N! I follow him still furt! 
using, not the super-Stirling approximation in which N In N is written for 
InNV!, but the better approximation in which (V InN — N) or (NInN — 

is written for In NV! The result is: 


S = (3/2)NkRInT + NkRIn V — NkIn N 


+ Nk In [(2rmk?e??/V oH] 35 


’ 


This quantity newly chosen as the picture of ‘‘entropy”’ depends on vo 
and on temperature in the right way, as did the other. The dependenc: 

N the number of atoms is now correct, and no wonder, for the new quantity 
was chosen with that purpose. There is a fourth term in the right-ha: 
member which is proportional to NV, and its value is completely determined 
if the value of VofZp is fixed. The value which it takes when N is made equa 
to Nomay be called ‘‘the chemical constant”’; but this name has been spoile 
through being used with several different meanings, and should probably }y 
abandoned. 

When to Vo/7y, the volume of the elementary cell in six-dimensional space, 
there is given the value /*—the cube of Planck’s constant—the resulting 
value of the fourth term is excellently confirmed by experiments on all of th: 
noble gases, and (with less precision) by experiments on many of the mona- 
tomic vapors of metallic elements. This is the achievement known as “‘the 
verification of the Sackur-Tetrode formula” and it is indeed a grand on 

Anyone versed in thermodynamics will probably regard this not as a grand 
result, but as an incomprehensible one! Are we not taught in thermody 
namics that nothing is ever measured about entropy except the differenc: 
between its values under different conditions, so an additive constant lik« 
the one in question must drop out of every verifiable equation, and its 
value can never be found? How then can it make sense to speak of co. 
firming the value of the fourth term on the right-hand side of (35)? 

Well, actually it 7s a difference which is measured: the difference betwee: 
the entropy of the gas at any convenient temperature and volume and the 
entropy of its solidified crystalline form at the absolute zero. This dif 
ference is found to be such, that if for the entropy of the gas one puts thi 
value (35) with /* substituted for Vof», then for the entropy of the crystal- 
line solid at the absolute zero one finds the value: sero. This result—this 
conclusion that the entropy of a crystal is zero at the absolute zero—is i 
itself so desirable and welcome that it is taken as the confirmation of the 
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Sackur-Tetrode formula. By “desirable and welcome’’ I mean that it is 
harmonious with the idea that entropy is a measure of disorder, an idea 
plausible in itself and fruitful in its applications. A chemical element per- 
fer tly crystallized at the absolute zero is supposed to be the exemplar of 
supreme order, and therefore its entropy ought to be nil. But this is an 
enormous subject requiring at least one other article, and I am glad that my 
attempt at writing such an article stands already in print in the June( 1942 
issue of this Journal. 

Here then is the astonishing history of the Classical Statistics. By a 
strangely artificial device, the numbering of atoms deemed identical, it 
arrived at the proper distributions—that is, the distributions ratified by 
experiment—in ordinary space and in momentum-space. It then proposed 
a picture of entropy partially right, yet wrong in its dependence on the 
number of atoms, and therefore fatally wrong. With another artificial and 
dubious device, it corrected itself by adopting a new picture of entropy, this 
time depending in the right way upon the number of atoms. With a third 
artificial device (the introduction of Planck’s constant in a peculiar way) it 
completed the formula for entropy in a manner leading to the consequence 
that the entropies of solidified crystallized elements are zero at absolute zero. 
All of these feats and more were subsequently achieved by the New Statistics, 
in a manner which I hope to explore on a later occasion. 

























Abstracts of Technical Articles by Bell System Author: 


Poles and Pole Treatment! Recinatp H. Cotiey. Studies ma 
pole use and drainage on the southern pine forest have brought out t! 
the demand for poles 35-feet and longer were to continue at the present 
a situation would soon develop in which these poles would be at a prem 
It would seem wise to use as many circumference classes as possible ai 
broaden the use of poles of other species whenever it is practicable to d: 
Successful full length treatment of red pine, lodgepole pine, western « 
and Douglas fir poles will help to broaden this use. 

There has been a definite trend toward greater mechanization in 
production. Machine shaving smooths the pole surface and accele: 
drying. Poles with square cut roofs and slab gains are all-purpose poles 
one design taking the place of four. New types of preservative treatment 
greensalt, creosote-petroleum-pentachlorphenol and _ salt-creosote—are 0! 
promise and must be considered, when current restrictions are removed, 
those cases where clean poles are mandatory. 

Interesting breaking test data are reported which show that the modulus 
of rupture of pole top sections average 90 per cent of the modulus of 1! 
poles as a whole, and that pole tops are sufficiently strong to meet thei! 
specified class breaking loads. This is of considerable importance wher 
poles are guyed. 

Current groundline treatment methods, it is pointed out, should be applic« 
to untreated poles in line that are worth saving at the time of inspection, as 
a part of the regular inspection procedure. 


Hearing, the Determining Factor for High-Fidelity Transmission”, Harve 
FLETCHER. This paper gives the requirements for ideal systems for th: 
transmission of speech and music. These requirements are based on: | 
Measurements of the threshold and frequency limits of the hearing of mor 
than 500,000 people at the New York and San Francisco World’s Fairs 
2. measurements of the discomfort level of sound; 3. measurements of roon 
noise in a wide variety of locations; and 4. measurements of the frequency) 
limits and the maximum and minimum levels of speech, orchestral music, 
and various instruments of the orchestra. 

From this information and from judgment tests it is concluded that 
substantially complete fidelity in the transmission of orchestral music ‘s 


1 Elec’l. Engg., Transactions Section, September 1942. 
2 Proc. I. R. E., June 1942. 
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ained by use of a system having a volume range of 65 decibels and a 
frequency range from 60 to 8000 cycles per second. Substantially complete 
lity for the transmission of speech is obtained by a system having a 


frequency range from 100 to 7000 cycles per second and a volume range of 
1) decibels. 

Preliminary experiments comparing a single-channel system and a two- 
channel stereophonic (auditory perspective) system showed that stereo- 
phonic transmission with an upper frequency limit of 5000 cycles per second 
was preferred to single-channel transmission with an upper limit of 15,000 
cycles per second. A definite improvement was obtained in the stereo- 
phonic system by using three channels instead of two. 

A New Direct Crystal-Controlled Oscillator for Ultra-Short-Wave Fre- 
quencies.» W.P. Mason and I. E. Farr. An ultra-high-frequency crystal 
oscillator is described which utilizes a mechanical harmonic of an AT or BT 
crystal. With the oscillator frequencies as high as 197 megacycles, hormon- 
ics as high as the 23rd have been excited. Taking the second electrical 
harmonic of the oscillator, frequencies as high as 300 megacycles, or 1 meter 
have been obtained. Since a mechanical harmonic is used, the crystal can be 
of a practical size to handle and adjust. The harmonic vibration of the AT 
and BT crystals have as low a temperature coefficient as the fundamental 
mode, and temperature coefficients of less than two parts per million per 
degree centigrade are easily obtained. Stability curves for this type of 
oscillator are shown and the results indicate that at 120 megacycles stabilities 
in the same order of magnitude as for ordinary crystal oscillators can be 
obtained. Without temperature or voltage control it appears likely that 
the frequency should remain constant to + 0.0025 per cent. 

Some measurements have been made of the properties of harmonic crystals 
at high frequencies. It was found that the Q of a crystal is independent of 
the frequency but in general increases with harmonic order. The ratio of 
capacitances r of a crystal increases as the square of the harmonic order. It 
is shown that in order to obtain a positive reactance in the crystal Q > 2r. 
Chis relation will only be satisfied for harmonics of AT crystals less than the 
7th. Asa result oscillator circuits such as the Pierce circuit cannot be used 
to drive crystals at high harmonic frequencies. A discussion of oscillator 
circuits is given and it is shown that a capacitance-bridge oscillator circuit 
with the crystal in one arm is the best type to use for high-frequency har- 
monic crystals. 


War Activities of the Bell Telephone System KeitH S. McHucu. The 
scope of the Bell System’s nation-wide service is, even in peace, difficult to 


° Proc. 1. R. E., October 1942. 
‘ Bell Telephone Magazine, November 1942. 
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visualize in its entirety. In war, when practically every phase of t 
tional effort to overthrow the Axis aggressors depends in some part 0 
communication, both the extent and the importance of the System 
tributions to the winning of the conflict are beyond summarizing. | 
past two years, numerous articles in the Bell Telephone Magazine (list 
the end of this article), and in the employee publications of the Associat; 
Companies, have described many aspects of the System’s cooperatio: 
the armed forces, with industry, and with the civilian population. No 
nearly a year after Pearl Harbor, it seems appropriate to review bot 
System’s preparations for the national emergency and the steps which 
taken since war became no longer a threat but a fact. To the extent 

it is possible in limited space, this article rounds out the previous frag 
mentary parts of the whole picture. 


The Number of Two-Terminal Series-Parallel Networks.’ Youn Rutos 
and C. E. SHANNON. This paper is concerned with the number of ways 
n abstract (electrical) elements may be connected in series-parallel arrang: 
ments and in particular with the way the number behaves for m large. After 
a proof of a generating identity for the numbers given without proof by I 
A. MacMahon in 1892, the paper gives recurrences and schemes of computa 
tion by means of which MacMahon’s table for the numbers is extended from 
n = 10 ton = 30. The behaviour for » large is shown to be of the form 


AN n-3?2 


> 


with A a fixed constant and \ a real number between 2 + 1/2 = 3.414 
and 4 and closer to the former than the latter; indeed an approximating 
function for which A is about 3.56 agrees with the numbers within 3°% over 
the range 7 to 20. These results are used to show that almost all switching 
functions of variables require at least 

n 


(i — «) — e>0 
logon 


switching elements (make or break contacts) in series-parallel realization. 


The Electrical Oscillations of a Perfectly Conducting Prolate Spheroid 
RoBert M. Ryper. The forced oscillations of a perfectly-conducting pro- 
late spheroid of eccentricity nearly unity are shown to be decomposable into 
“harmonics” corresponding to different modes of vibration, each harmoni: 
being quantitatively connected with a certain portion of the impressed elec- 
tric field which drives the antenna. The harmonics contribute additivels 


to the current and field of the spheroid; each offers a characteristic imped- 


5 Jour. Mathematics and Physics, August 1942, 
6 Jour. Applied Physics, May 1942. 
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ance to the driving field, and the properties of the antenna are a composite 
depending upon the proportions of the various harmonics present. The 


) behavior of the harmonics with frequency is discussed qualitatively; analyti- 
}cal expressions obtained are useful chiefly at the resonant frequencies of the 
‘antenna, where the most important harmonic becomes sinusoidal in char- 


t 


p acter. 


On Radiation from Antennas.” S. A. ScueLKuNorr and C. B. FELDMAN. 
[his paper presents some theoretical remarks and experimental data relating 
to applications of the transmission-line theory to antennas. It is empha- 
sized that the voltage, the current, and the charge are affected by radiation 
in different ways, a fact which should be considered in any adaptation of 
line equations to antennas. 

It is shown experimentally and theoretically that in an antenna of length 
equal to an integral number of half wave-lengths, which is energized at a 
current antinode, the effect of radiation on the current and the charge 
but not on the voltage) can roughly be represented by adding to the resist- 


ance of the wires another fairly simple term. 


The Use of Secondary Electron Emission to Obtain Trigger or Relay Action.® 
\. M. SkeLtetr. The use of secondary electrons to obtain trigger action 
similar to that of a thyratron is described. An experimental tube and the 
necessary circuits by which this action is achieved are discussed. This 
combination gives the features of a triode with a relay or on and off feature, 
resulting in an amplifier, oscillator, modulator, or other vacuum tube device 
which may be turned on or off abruptly at high or low frequencies. In 
addition, it can be used to replace thyratrons in many of their circuits where 
very low impedance is not necessary and is capable of much greater speeds 
of operation in such applications. 


{ New Frequency-Modulation Broadcasting Transmitter.2 A, A. SKENI 
and N. C. O_msteap. A new frequency-modulation transmitter is de- 
scribed which uses a novel amplifier circuit permitting an unusually simple 
mechanical design and an economical vacuum-tube complement. 

The choice and design of circuit components, governed by both mechani- 


cal and electrical considerations, are discussed in detail. 


A Secondary Frequency Standard Using Regenerative Frequency-Dividing 
Circuits. F, R. STANSEL. A secondary frequency standard is described 


7 Proc. I. R. E., November 1942 

’ Jour. Applied Physics, August 1942. 
’ Proc. I. R. E., July 1942. 

10 Proc. I. R. E., April 1942. 
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in which standard frequencies are derived from a 5-megacycle oscill: 
a series of frequency dividers. The advantage of obtaining standa 
quencies by frequency division rather than by frequency multiplica 
pointed out and the characteristics of the regenerative frequency di 
used are discussed. 


Some Mechanical Aspects of Telephone Apparatus J. D. TEBo aid | 
G. MenHLHouse. Pari J. It is seldom realized that the vastness : 


Bell System requires such an enormous amount of equipment and wires {o: 
handling the 100,000,000 calls per day made by the people in the | 
States. A total of 44,000 kinds of apparatus involving 170,000 difi 
parts are required. The crossbar switch, described in this article, is | 
specific telephone switching mechanism used in the latest dial tele 
system. 

Of particular interest in the crossbar switch is the design of contact spri: 
These springs are essentially thin, metal beams of a rectangular cross sect 
but of varying cross sectional area along their length. They are clamp 
at one end and are subject to bending as compound cantilevers. Si 
the clamping is necessarily between relatively soft insulating mater 
the determination of the effective length of the spring in determining its 
true deflection curve required the use of unique methods. Likewise, si: 
the springs are punched out of sheet stock at an angle to the grain directio 
of the material, the modulus of elasticity does not remain the same for equ 
cross sections of the same material. Consideration of these points wa 
necessary in determining the strength of the magnets for operating the co: 
tact springs, as well as to insure that the stresses introduced in the springs 
would not be excessive. 

To study the motions of the various parts of the switch, both high spe: 
rapid record” oscillograph were used. The os 


se 


motion pictures and the 
lograph was provided with means for obtaining ‘“‘ shadowgrams”’ of the actua 
movement of parts simultaneously with the changes in the electrical chara 
teristics of the magnets and contacts. The use of both the camera a1 
oscillograph provided valuable data for making improvements, both 
design and operating characteristics. 

Part IT. The manufacture of crossbar apparatus is accomplished 01 
product basis; that is, the entire range of manufacturing operations is segre 
gated into one division, and practically all operations from raw mate! 
to the completed product are performed in this division. More than 15) 
kinds of parts totaling an annual demand of 200,000,000 individual pieces 
are required to produce the crossbar switches. 

To produce these parts requires a number of special machines, tools, a 


1 Mech. Engg., May 1942 and June 1942. 
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operations. Of particular interest are the use of special progressive punch 
and die tools, coil winding, contact welding and conveyorized assembly and 
{justment. 

lhe use of heavy presses with large progressive punches and dies was 


necessitated by the degree of accuracy required for this grade of equipment. 
For example, the vertical unit base, weighing only ? lb., is produced by a 
75,000 Ib. press, using a 3,500 lb. tool. Again, since 40,000,000 contact 
springs are required annually, each of which must be attached to insulators, 
automatic presses, conveyor belts and handling devices are required. Weld- 
ing two contacts on each spring is accomplished by special welding presses 
using rolls of contact metal tape, each contact being cut off just prior to 
the welding process. A quality of less than one defective contact out of 
20,000 is maintained. 

The coils are wound in special machines in ‘“‘sticks” of 5 to 7 coils with 
0007 inch thick cellulose acetate between each layer. A wire of a gage size 
halfway between *37 and #38 B &S is used on most of the coils because 
of critical capability conditions—the annual amount of wire of this size 
being 125,000 Ibs. 

Assembly and adjustment of the switches proceeds on conveyor belts 
from one end of a large room to the other—the procedure being set up in 
such a way as to create a continuous flow of completed parts for wiring into 
the large frames ready for installation in telephone exchanges. 


Regulated Rectifiers in Telephone Offices." D. E. Trucksess. For many 
years rectifiers of the garage type were used in converting alternating cur- 
rent to direct current for charging batteries used for communication pur- 
poses. These batteries furnish power for relay operation, for talking, and 
filament and plate supplies for repeaters. The rectifiers were of the manual- 
control type where the operator selected the charging current by means of 
tap switches or rheostats. 

With the development of the thyratron type of tube, a rectifying means 
was made available in which the grid of the rectifier tube could be used to 
control its own output current by an electronic circuit. Rectifier circuits 
were designed to maintain a constant output voltage. If a regulated recti- 
fer is connected to a battery and the constant rectifier voltage is 2.15 volts 
per cell, the load current will automatically come from the rectifier and not 
from the battery. Also the battery will draw from the rectifier sufficient 
additional current to maintain its charge. If the circuit voltage is held 
within limits of less than plus or minus one per cent, the maintenance of 
the battery is reduced and its life is extended. 

The thyratron tube differs from the vacuum tube in that the grid does 


2 Flec’l. Engg., Transactions Section, August 1942. 
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not have a continuous control of the plate current. When a positive | 


tial is applied to the plate, current does not flow until the magnitude 
negative grid voltage is reduced to the critical value, at which time the 
current flows, and the magnitude of the plate current depends upo 
voltages and impedances in the circuit. The grid has no further co 
and plate current flows until it is stopped by reducing the plate v: 
to zero. 

Thyratron tubes use various gases and mixtures of gases. The earlies 
type used mercury vapor, but this type of tube is quite sensitive to temper 
ture changes. The grid characteristics are shifted materially by cha 
in the room temperature in which it is operated, and in low temperatures 
it is almost a vacuum tube. Thyratron tubes using argon gas are 
affected by temperature changes, but high-pressure argon tubes have a 
inverse voltage which limits their application to low-voltage rectifiers 
Tubes using low-pressure argon have a higher inverse voltage, but are 
companied by a high are drop which makes their efficiency low. A mixtur 
of mercury vapor and argon has been found which provides the temperatu: 
stable grid characteristic of the argon tube and the low arc drop of the mer 
cury-vapor tube. This type of tube has been very successful with certai 
regulating circuits, particularly at voltages less than 60 volts. 

Five kinds of regulating circuits are used in telephone offices to hold th 
output voltage of rectifiers constant. The selection of the circuit to be used 
depends upon the magnitude of the current, d-c voltage, and type of recti- 
fying means to be used. Two forms of regulating circuits using thyratron 
tubes and one using two-element high-pressure tubes were developed. A 
fourth circuit using all vacuum tubes was adapted for telephone use. Th« 
fifth kind uses a negative resistance. 

In this paper a table shows the voltage and current output, type of contro 
of the rectifiers, and the rectifying means that have found widespread use 
in the Bell System. 

The regulated rectifier finds its applications in telephone offices where 
constant voltage, independent of load and a-c line-voltage variations, is 
required to supply filament grid bias and plate voltage to telephone re- 
peaters. Certain measuring circuits require a regulated rectifier to supp!) 
a stabilized voltage. Regulated rectifiers also find applications where con 
stant voltage is of secondary importance but an automatic power plant is 
desired for maintaining storage batteries in a fully charged condition to be 
ready to supply the power for telephone offices if the a-c power fails. A 
further compensation of regulating the voltage is the increase in life ob 
tained from storage batteries if they are not continually being charged and 
discharged but are fully floated. 
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that there was no April 1943 issue of the Bell System 
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It is also to be noted that there was only one issue of 
Volume 21 of the Bell System Technical Journal. It 
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